Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 14,124 T→A 5.7% intergenic (+45/‑44) dnaK → / → dnaJ chaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091314,1240TA 5.7% 255.1 / 2.8 91intergenic (+45/‑44)dnaK/dnaJchaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone
Reads supporting (aligned to +/- strand):  ref base T (41/44);  new base A (3/3);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

AACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGA  >  NC_000913/13976‑14271
                                                                                                                                                    |                                                                                                                                                   
aaCATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAt                                                                                                                                                     <  1:900213/149‑1 (MQ=255)
aaCATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAt                                                                                                                                                     <  1:284722/149‑1 (MQ=255)
 ccATGCCCAGCAGCAGACTGCCGGTGCTGCTGCTTCTGCAAACAACGCGAAAGATGACGATGTTGGCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATTAACGGGTAATTATACTGACACGGGCGAAGGGGAAtt                                                                                                                                                    >  2:1041722/2‑149 (MQ=255)
 aCATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAtt                                                                                                                                                    >  2:127057/1‑149 (MQ=255)
  cATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAttt                                                                                                                                                   >  2:363653/1‑149 (MQ=255)
   aTGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTc                                                                                                                                                  >  2:611222/1‑149 (MQ=255)
   aTGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGGCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTc                                                                                                                                                  <  1:35386/149‑1 (MQ=255)
     gCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCct                                                                                                                                                >  2:419303/1‑149 (MQ=255)
           cagctgACTGCCTGTGCTGATGCTTCTGCAAACAACGCGATAGATGACGATGTTGGCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGc                                                                                                                                          <  1:548992/149‑1 (MQ=255)
            agcagACTGCCGGTGCTTATGCTTCAGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAAAAATCGCCCTATAAACGGGTAAATATACTGACACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                                                                         <  2:1066817/149‑1 (MQ=255)
             ccagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGAGGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGccc                                                                                                                                        <  1:1088900/148‑1 (MQ=255)
              cagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGCCGAAGGGGAATTTCCTCTCCGCCCg                                                                                                                                       >  2:120117/1‑149 (MQ=255)
               agACTGCCGTTGCCGATGCTTCTGCAAACAACGCGAAAGATGACGATGGCGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTCATTATCCTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                                                                      <  1:1177769/149‑1 (MQ=255)
                 aCTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTATTTATACTGACACGGGCGAAGGGGAATTTCATCTCCTCCCGtcc                                                                                                                                    >  1:791268/1‑147 (MQ=255)
                   tGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCAt                                                                                                                                  >  1:755644/1‑149 (MQ=255)
                           cTGATGCTTCTGCAAACAACGCTAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATATCCTCTCCGCCCGTGCATTCATCTAg                                                                                                                          <  1:910689/149‑1 (MQ=255)
                              aTGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAgggg                                                                                                                       >  2:1001044/1‑149 (MQ=255)
                                  ttCTGCAAACAACGCGAAAGCTGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAt                                                                                                                   >  1:176426/1‑149 (MQ=255)
                                  ttCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAt                                                                                                                   >  1:1117416/1‑149 (MQ=255)
                                      gCAAACAACGCGAATGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAACGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGCGAATTTCCTCTCCGCCCTTGCATTCATCTAGGGGCTATTTaa                                                                                                               >  1:781892/1‑149 (MQ=255)
                                        aaaCAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaa                                                                                                             <  1:609637/149‑1 (MQ=255)
                                        aaaCAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaa                                                                                                             <  1:1058821/149‑1 (MQ=255)
                                          acaacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaaa                                                                                                           >  1:303851/1‑149 (MQ=255)
                                          acaacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGCCGAAGGGGAATTTCGTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaaa                                                                                                           >  2:363032/1‑149 (MQ=255)
                                            aacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                                                                         <  2:723094/149‑1 (MQ=255)
                                            aacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                                                                         >  1:873460/1‑149 (MQ=255)
                                             acGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                                                                        <  1:362852/149‑1 (MQ=255)
                                                 gAAATAAGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGCCAATAAATATTGGCCCTATAAACGGTTTATTATACTGACAAGGGCGAAGGGGAATTTCCTGTCCGCCCGTGCATTCATCTAGGGGCAATTTATAAATGATGGCt                                                                                                    <  1:240475/149‑1 (MQ=255)
                                                 gAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCt                                                                                                    >  1:942262/1‑149 (MQ=255)
                                                 gAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCATCTCCGCCCGTGCATTCATCTAGGGGCAATATAAAAAAGATGGCt                                                                                                    >  1:307302/1‑149 (MQ=255)
                                                   aaGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTGTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTaa                                                                                                  <  1:788177/149‑1 (MQ=255)
                                                    aGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAg                                                                                                 >  2:778669/1‑149 (MQ=255)
                                                     gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTGTAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGc                                                                                                <  1:372139/149‑1 (MQ=255)
                                                     gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGc                                                                                                <  1:127057/149‑1 (MQ=255)
                                                        gACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAg                                                                                             >  2:204730/1‑149 (MQ=255)
                                                          cGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCTTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                                                                           >  2:603976/1‑149 (MQ=255)
                                                           gATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCAGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatt                                                                                          >  2:1065607/1‑149 (MQ=255)
                                                              gttgtCGACGCTGTATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAAGGGGTAATTATAGTGACTCGGGAGTAGGTGAATTTCCTCTCCGCCCGTGCATTCATCTATTGGCAATTTAAAAAAGATGGCTAAGCAAGattatt                                                                                       <  2:526371/149‑1 (MQ=255)
                                                              gttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGattatt                                                                                       <  1:1127813/149‑1 (MQ=255)
                                                               ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACTGGTAATTATACTGACCCGGGCGAAGGGGAATTTCCTCTCCGCCGGTGCATTCTTCTAGGGGCAATTTAAAAAAGATGGGTAAGCAAGAttatta                                                                                      <  2:995438/149‑1 (MQ=255)
                                                                tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTAc                                                                                     <  1:224451/149‑1 (MQ=255)
                                                                 gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTGTCCGCCCGTGCATTAATCTAGGGCAAATTTAAAAAAGATGGCTAAGCATGTTTATTAcc                                                                                    >  2:660934/1‑148 (MQ=255)
                                                                 gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACg                                                                                    >  1:21593/1‑149 (MQ=255)
                                                                    gACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACgaga                                                                                 >  1:1177044/1‑149 (MQ=255)
                                                                     aCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAt                                                                                >  1:261301/1‑149 (MQ=255)
                                                                        atcaaTTTCCAGAAGTCAATGAAAAAAAATAATCGCCCTAAAAACGGCTAATTATACTTACAAGTGCGAAGGGGAATTTCCTCTACGCCCGTGCATTCATCTAGTTGAAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAtttt                                                                             <  2:768308/146‑1 (MQ=255)
                                                                         tGAATTTGAAGAAGTCAAAGCCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCAGTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                                                            >  2:526896/1‑149 (MQ=255)
                                                                           aaTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                          <  2:596181/149‑1 (MQ=255)
                                                                           aaTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                          <  1:452322/149‑1 (MQ=255)
                                                                             tttGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGTgtca                                                                                                                                  <  1:1132508/91‑4 (MQ=25)
                                                                             tttGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGTgtca                                                                                                                                  >  2:1132508/1‑88 (MQ=25)
                                                                               tGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGtt                                                                      >  2:350661/1‑149 (MQ=255)
                                                                                gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGttt                                                                     <  1:947493/149‑1 (MQ=255)
                                                                                gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGttt                                                                     >  1:494028/1‑149 (MQ=255)
                                                                                gaagaaGTCAAAGACAAAAAATAATCGCACTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCtttt                                                                     >  2:302551/1‑149 (MQ=255)
                                                                                 aagaagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTc                                                                    <  1:172104/149‑1 (MQ=255)
                                                                                    aagTCAAACACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaa                                                                 <  2:568609/149‑1 (MQ=255)
                                                                                      gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa                                                               <  1:290868/149‑1 (MQ=255)
                                                                                      gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa                                                               >  1:75820/1‑149 (MQ=255)
                                                                                        cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACa                                                             >  1:389424/1‑149 (MQ=255)
                                                                                        cAAAGACAAAAAATAAGCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAGGGCTACGCCAGATTATTACGAGATTTTAGGCGTTTCCAAAACa                                                             <  1:206645/149‑1 (MQ=255)
                                                                                         aaaGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAg                                                            <  2:656595/149‑1 (MQ=255)
                                                                                         aaaGACAAAAAATAATCGCCCTATAAACGGGTAATTATACGGACACCGGCGAAGGGGAATTTCCTCTCCGCCAGTGCATACTTCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACTAGAGTTTAGGATTTTCCAAAACat                                                            >  1:683795/1‑148 (MQ=255)
                                                                                           aGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCg                                                          <  1:1109858/149‑1 (MQ=255)
                                                                                            gACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCgg                                                         <  1:1001044/149‑1 (MQ=255)
                                                                                             aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGTCGTTTCCACAACAGCGGa                                                        >  1:178601/1‑149 (MQ=255)
                                                                                              cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCAGTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTATCCAAAACAGCGGaa                                                       >  1:839428/1‑149 (MQ=255)
                                                                                                aaaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAaga                                                     <  2:923239/149‑1 (MQ=255)
                                                                                                 aaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAagac                                                    >  2:418532/1‑148 (MQ=255)
                                                                                                 aaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAGTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAagag                                                    >  2:889523/1‑149 (MQ=255)
                                                                                                 aaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAATATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAagac                                                    >  1:1172851/1‑148 (MQ=255)
                                                                                                    ataatCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGt                                                 <  1:500939/149‑1 (MQ=255)
                                                                                                       atCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGaa                                              <  1:204730/149‑1 (MQ=255)
                                                                                                         cGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAAt                                            <  1:363653/149‑1 (MQ=255)
                                                                                                          tcccTATAAACGGGTAATTATATTGACACGGGCGAAGGTGAATTTCCTCTCAGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATc                                           <  2:247476/148‑1 (MQ=255)
                                                                                                              tataAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTACTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCTAGATTATTACTAGATTTTAGGAGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaa                                       >  1:59239/1‑149 (MQ=255)
                                                                                                               ataAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaaa                                      >  1:1177249/1‑149 (MQ=255)
                                                                                                                   aCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTTTTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGTcc                                  >  1:1117218/1‑149 (MQ=255)
                                                                                                                    cGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCt                                 <  1:724384/149‑1 (MQ=255)
                                                                                                                     gggTAATTATACTGACACGGGCGAAGGGGAATTTCGTCTCCGCCCGTGCATTCATCTAGGGGCAATATAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTa                                <  2:1163392/149‑1 (MQ=255)
                                                                                                                     gggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTa                                <  1:937262/149‑1 (MQ=255)
                                                                                                                     gggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTa                                >  1:961842/1‑149 (MQ=255)
                                                                                                                      ggTAACTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTAc                               <  1:550521/149‑1 (MQ=255)
                                                                                                                       gTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACa                              >  1:38847/1‑149 (MQ=255)
                                                                                                                            tataCTGACACTGGCCTAGGGGAATTTCCTCTCCGCGCGTGCATTCATCTAGGCGCACTTTAAAAAAGCTGGCTTAGCAAGATTTTTACGAGATTTTCGGCGTTTCCAACACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGc                         <  2:898136/149‑1 (MQ=255)
                                                                                                                            tataCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGc                         <  1:412355/149‑1 (MQ=255)
                                                                                                                                     acGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg                >  2:507573/1‑149 (MQ=255)
                                                                                                                                     acGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg                >  2:503152/1‑149 (MQ=255)
                                                                                                                                      cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGa               <  2:746640/149‑1 (MQ=255)
                                                                                                                                       gggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGaa              <  2:791268/149‑1 (MQ=255)
                                                                                                                                         gCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAAt            >  2:447900/1‑149 (MQ=255)
                                                                                                                                          cGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATa           <  1:611222/149‑1 (MQ=255)
                                                                                                                                           gAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATac          <  1:350661/149‑1 (MQ=255)
                                                                                                                                                gtaaTTTCCTCACAGCCAGGTCATTCATCTAGGGGCAATTTAAAAAACATGGCTAAGCAAGATTATTAGGATATCTTAGGCGTTTGCAAAACAGCGTAAGAGCGTGAAATCAGAATGGCCTTCAAACGCCTGGCCATGAAATAACAccc     <  2:188329/147‑1 (MQ=255)
                                                                                                                                                 gAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCg    >  1:1053400/1‑149 (MQ=255)
                                                                                                                                                   aTTTCCTCTCCGCCCGTGCATTAATCTAGGTGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGa  <  2:839428/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGA  >  NC_000913/13976‑14271

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: