Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,154 T→A 8.2% intergenic (‑124/+22) dsbB ← / ← nhaB protein thiol:quinone oxidoreductase DsbB/Na(+):H(+) antiporter NhaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,233,1540TA 8.2% 289.7 / 10.9 96intergenic (‑124/+22)dsbB/nhaBprotein thiol:quinone oxidoreductase DsbB/Na(+):H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base T (39/49);  new base A (4/4);  total (43/53)
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.34e-01

GAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGC  >  NC_000913/1233005‑1233303
                                                                                                                                                     |                                                                                                                                                     
gAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCAtt                                                                                                                                                       >  2:1372962/1‑150 (MQ=255)
   cATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCTCCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtat                                                                                                                                                    >  2:699703/1‑150 (MQ=255)
    aTTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATTCATAGAGTGCGAAGTTTAAACCATTTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGAAGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                                                                   <  2:528856/150‑1 (MQ=255)
    aTTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTtata                                                                                                                                                   <  1:135676/150‑1 (MQ=255)
        gTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                                                                               >  1:856959/1‑150 (MQ=255)
        gTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGcc                                                                                                                                               >  2:1391407/1‑150 (MQ=255)
         ttCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGccc                                                                                                                                              <  2:1429897/150‑1 (MQ=255)
          tCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGTACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCg                                                                                                                                             <  2:1930638/150‑1 (MQ=255)
           cAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCgg                                                                                                                                            <  1:1731643/150‑1 (MQ=255)
            aaaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGTCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGt                                                                                                                                           >  2:1014974/1‑150 (MQ=255)
            aaaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGt                                                                                                                                           <  2:236716/150‑1 (MQ=255)
             aaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                                                                          <  2:992347/150‑1 (MQ=255)
             aaaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTa                                                                                                                                          <  1:4310/150‑1 (MQ=255)
              aaaTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTaa                                                                                                                                         >  1:866974/1‑150 (MQ=255)
              aaaTCGCAACATAATAATTTCCCTTCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGttt                                                                                                                                         <  1:1847353/150‑3 (MQ=255)
                aTCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgc                                                                                                                                       <  2:478830/150‑5 (MQ=255)
                  cGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATATAAATGCCCGGtttgcct                                                                                                                                     <  1:1338968/150‑7 (MQ=255)
                   gCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAgttgt                                                                                                                                    >  1:1052510/1‑150 (MQ=255)
                     aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCTCCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                                                                  <  1:699703/150‑1 (MQ=255)
                     aaCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTa                                                                                                                                  <  1:1367516/150‑1 (MQ=255)
                      aCATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTAt                                                                                                                                 <  2:52435/150‑1 (MQ=255)
                       cATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATc                                                                                                                                <  2:1463007/150‑1 (MQ=255)
                        aTAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCa                                                                                                                               <  1:346957/150‑1 (MQ=255)
                          aaTCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGccc                                                                                                                                              >  2:31266/1‑133 (MQ=255)
                          aaTCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGccc                                                                                                                                              <  1:31266/133‑1 (MQ=255)
                           aTCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGTCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAg                                                                                                                            <  2:1761339/150‑1 (MQ=255)
                                 tCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGc                                                                                                                      <  2:89864/150‑1 (MQ=255)
                                    cTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTAt                                                                                                                   >  1:820470/1‑150 (MQ=255)
                                      gCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                                                                                 <  2:253749/150‑1 (MQ=255)
                                      gCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                                                                                 >  1:1079862/1‑150 (MQ=255)
                                      gCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                                                                                 >  1:755405/1‑150 (MQ=255)
                                            tGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCa                                                                                                           >  1:553928/1‑150 (MQ=255)
                                             gCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                                                                          >  2:1480967/1‑150 (MQ=255)
                                             gCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCAt                                                                                                          >  2:1512449/1‑150 (MQ=255)
                                                  gagTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCa                                                                                                     >  2:1724616/1‑150 (MQ=255)
                                                    gTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATa                                                                                                   <  1:1387426/150‑1 (MQ=255)
                                                      gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                                                                                 >  2:1837827/1‑150 (MQ=255)
                                                      gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                                                                                 >  1:1125674/1‑150 (MQ=255)
                                                          aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                                                                       <  2:699480/140‑1 (MQ=255)
                                                          aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                                                                       >  1:699480/1‑140 (MQ=255)
                                                           gTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCAtt                                                                                            >  2:727713/1‑150 (MQ=255)
                                                             ttAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCg                                                                                          <  2:1079862/150‑1 (MQ=255)
                                                               aaaCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGt                                                                                        >  2:1120811/1‑150 (MQ=255)
                                                                aaCCAATTCATTAGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCGGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGtt                                                                                       >  1:1655063/1‑150 (MQ=255)
                                                                 aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAg                                                                                                               <  2:128756/125‑1 (MQ=255)
                                                                 aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAg                                                                                                               >  1:128756/1‑125 (MQ=255)
                                                                 aCCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTa                                                                                      <  2:1431924/150‑1 (MQ=255)
                                                                  ccAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTAc                                                                                     >  1:1227691/1‑150 (MQ=255)
                                                                   cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACa                                                                                    >  1:1922419/1‑150 (MQ=255)
                                                                   cAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACa                                                                                    <  1:1187504/150‑1 (MQ=255)
                                                                     attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgg                                                                                  <  1:1261155/150‑1 (MQ=255)
                                                                     attcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgg                                                                                  <  2:1052510/150‑1 (MQ=255)
                                                                      ttcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAggg                                                                                 <  1:962928/150‑1 (MQ=255)
                                                                       tcattcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAgggg                                                                                <  2:153640/150‑1 (MQ=255)
                                                                          ttcGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCaa                                                                             >  1:1723526/1‑150 (MQ=255)
                                                                              cTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGt                                                                         >  1:496194/1‑150 (MQ=255)
                                                                                gcgcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTaa                                                                       <  1:1816339/150‑1 (MQ=255)
                                                                                  gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                                     <  1:1541324/150‑1 (MQ=255)
                                                                                  gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                                     >  1:1248331/1‑150 (MQ=255)
                                                                                  gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                                     <  1:1391407/150‑1 (MQ=255)
                                                                                  gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                                     >  1:593836/1‑150 (MQ=255)
                                                                                   ccaccaGGAGGCATCAAAAATAGGACTCAATTGGACAATTATCGGGCGAAAATGCAAACCGCGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAAGCATTCGGTTACAGGGGCAAGCGTAAACt                                                                    <  1:334306/150‑1 (MQ=255)
                                                                                   ccaccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACt                                                                    <  1:334304/150‑1 (MQ=255)
                                                                                    caccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTc                                                                   <  2:254454/150‑1 (MQ=255)
                                                                                      ccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa                                                                 >  1:992629/1‑150 (MQ=255)
                                                                                      ccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa                                                                 >  1:420637/1‑150 (MQ=255)
                                                                                         ggaggCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGc                                                              <  1:1480967/150‑1 (MQ=255)
                                                                                             gcagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCagag                                                          <  2:554750/150‑1 (MQ=255)
                                                                                             gcagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCagag                                                          <  1:716422/150‑1 (MQ=255)
                                                                                              cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGc                                                         >  1:716279/1‑150 (MQ=255)
                                                                                              cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGc                                                         <  2:553928/150‑1 (MQ=255)
                                                                                              cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGc                                                         >  2:644426/1‑150 (MQ=255)
                                                                                                gcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCaa                                                       >  1:1302725/1‑150 (MQ=255)
                                                                                                  aaaaaTAGGACGCAATTGGACAATTATCTGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGc                                                     <  2:1319032/150‑1 (MQ=255)
                                                                                                  aaaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGc                                                     <  2:1186316/150‑1 (MQ=255)
                                                                                                   aaaaTAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGcc                                                    <  1:232385/150‑1 (MQ=255)
                                                                                                         ggACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAg                                              >  2:584716/1‑150 (MQ=255)
                                                                                                           aCGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAgtg                                            <  2:716279/150‑1 (MQ=255)
                                                                                                           aCGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAgtg                                            >  1:1690590/1‑150 (MQ=255)
                                                                                                           aCGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAgtg                                            <  1:1512449/150‑1 (MQ=255)
                                                                                                              cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCa                                         <  1:1014974/150‑1 (MQ=255)
                                                                                                              cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCa                                         >  2:215785/1‑150 (MQ=255)
                                                                                                              cAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGATATCAAAGCGTTGCTAACCAGAACATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAATCGACGCAGAGCAAGCCGACGAGTGTCa                                         >  1:1291826/1‑150 (MQ=255)
                                                                                                               aaTTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAg                                        <  2:1285347/150‑1 (MQ=255)
                                                                                                                      cAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAgg                                 >  1:44235/1‑150 (MQ=255)
                                                                                                                         ttATCGGTCTAAAAGGCAAACCGGGCATTTATATTCCCGGTAAGTTTTTATCATAGCGTTTCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGAAAGCGTATTCTCGACGGAGAGCTAGACTACGAGTGTCAGGAGGAGGttg                              <  2:1346542/150‑1 (MQ=255)
                                                                                                                                   aaaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAggg                    <  2:866974/150‑1 (MQ=255)
                                                                                                                                   aaaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATAAAAGCGTTGCTATACAGACCATTTGCATAAACCATTCGGTTACAGGCGCAAGCCTAAAAACGACGCAGAGCAAGCCGACGAGTGTCAGGACCAGGGACAAAGGCAggc                    >  2:1368410/1‑149 (MQ=255)
                                                                                                                                    aaaGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGc                   <  2:1227691/150‑1 (MQ=255)
                                                                                                                                        gCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATc               <  2:820470/150‑1 (MQ=255)
                                                                                                                                        acttaCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATc               >  2:1604398/5‑150 (MQ=255)
                                                                                                                                         cAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATcc              <  1:1368410/150‑1 (MQ=255)
                                                                                                                                           taCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCcac            >  1:1886975/2‑150 (MQ=255)
                                                                                                                                           taCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCcac            >  2:219652/2‑150 (MQ=255)
                                                                                                                                            aCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCcaca           >  2:1539780/1‑150 (MQ=255)
                                                                                                                                              cGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACAcc         >  1:434396/1‑150 (MQ=255)
                                                                                                                                              cGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACAcc         <  1:1372962/150‑1 (MQ=255)
                                                                                                                                                     ttATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGc  <  2:125874/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGC  >  NC_000913/1233005‑1233303

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