Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,233,157 T→A 8.5% intergenic (‑127/+19) 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,1570TA 8.5% 281.3 / 12.8 94intergenic (‑127/+19)dsbB/nhaBprotein thiol:quinone oxidoreductase DsbB/Na(+):H(+) antiporter NhaB
Reads supporting (aligned to +/- strand):  ref base T (38/48);  new base A (4/4);  total (42/52)
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.31e-01

CATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCAT  >  NC_000913/1233008‑1233306
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                                    >  1:755405/1‑150 (MQ=255)
                                   gCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                                                                                    >  1:1079862/1‑150 (MQ=255)
                                   gCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATcc                                                                                                                    <  2:253749/150‑1 (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                                                                                                    >  1:1125674/1‑150 (MQ=255)
                                                   gcgcAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATaaa                                                                                                    >  2:1837827/1‑150 (MQ=255)
                                                       aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                                                                          >  1:699480/1‑140 (MQ=255)
                                                       aGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTg                                                                                                          <  2:699480/140‑1 (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:593836/1‑150 (MQ=255)
                                                                               gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                                        <  1:1391407/150‑1 (MQ=255)
                                                                               gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                                        >  1:1248331/1‑150 (MQ=255)
                                                                               gcCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAAc                                                                        <  1:1541324/150‑1 (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:420637/1‑150 (MQ=255)
                                                                                   ccaGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGa                                                                    >  1:992629/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                                                            <  2:553928/150‑1 (MQ=255)
                                                                                           cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGc                                                            >  2:644426/1‑150 (MQ=255)
                                                                                           cagcaAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGc                                                            >  1:716279/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               >  2:219652/2‑150 (MQ=255)
                                                                                                                                        taCCGGGCATATAAATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCcac               >  1:1886975/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)
                                                                                                                                                     tatGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCAt  >  2:102681/1‑150 (MQ=255)
                                                                                                                                                     tatGCCCGGTAAGTTGTTATCAAAGCGTTGCTATACAGCCCATTTGCATAAACCATTCGGTTACAGTGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGATTGTCAGGACGAGGGTGTAAGTCAGGGCCATCCACACCATGCGGACaa  >  2:1838309/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATTGGTTCAAAAATCGCAACATAATCATTTCCCTGCAATATGCATAGAGTGCGCAGTTTAAACCAATTCATTCGCTGCGCCACCAGGAGGCAGCAAAAATAGGACGCAATTGGACAATTATCGGGCGAAAAGGCAAACCGGGCATTTATATGCCCGGTAAGTTGTTATCAAAGCGTTGCTATCCAGCCCATTTGCATAAACCATTCGGTTACAGGGGCAAGCGTAAACTCGACGCAGAGCAAGCCGACGAGTGTCAGGACGAGGGTGTAAGGCAGGGCCATCCACACCATGCGGCCAT  >  NC_000913/1233008‑1233306

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: