Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,139,787 C→G 5.1% D470E (GAC→GAG flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,139,7870CG 5.1% 365.7 / 9.4 117D470E (GAC→GAGflgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (49/62);  new base G (3/3);  total (52/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 120 total aligned reads displayed.

CAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCA  >  NC_000913/1139640‑1139936
                                                                                                                                                   |                                                                                                                                                     
cAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGc                                                                                                                                                     >  1:732456/1‑150 (MQ=255)
 aaCAGGGATGTAGTAATCACCGACGAGGCGGAAGTGGGGGTGGCGAGCGAAGAAGAGGCGGGTGATAGCGATAACCGCAGCGGTCAGGCCCTGCTGGATCTGCAGGGCAACGGTGGTGCGGTGGGCGGTGCGAAGTCCTTTAACGACGcg                                                                                                                                                    >  2:1293838/1‑149 (MQ=255)
   cATTGATGTATTAAACACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTATAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTa                                                                                                                                                  <  1:1047319/150‑1 (MQ=255)
     tGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATg                                                                                                                                                >  1:446313/1‑150 (MQ=255)
     tGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATg                                                                                                                                                <  2:1878313/150‑1 (MQ=255)
     tGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATg                                                                                                                                                <  1:1503330/150‑1 (MQ=255)
     tGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGACCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATg                                                                                                                                                >  1:1198482/1‑150 (MQ=255)
            aTTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATATGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTACCGACGCTTATGCCTCGtt                                                                                                                                         >  2:667699/1‑150 (MQ=255)
              tAATCACCGACGAAGCGAAAATAGCGATGGCTAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTTGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTag                                                                                                                                       <  2:1289924/150‑1 (MQ=255)
                aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                     <  2:1381870/150‑1 (MQ=255)
                aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                     <  1:1470669/150‑1 (MQ=255)
                aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                     >  2:1520821/1‑150 (MQ=255)
                 tCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtga                                                                                                                                    >  1:1792004/1‑150 (MQ=255)
                 tCACCGACGAAGCGAAAATAGAGATGGCTAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGAGGCATAAGCGTCGTTaaagg                                                                                                                                    <  1:1542993/150‑5 (MQ=255)
                 tCAACGACTAAGCGAAAATAGCGATTGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGGGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtga                                                                                                                                    <  1:172160/150‑1 (MQ=255)
                  ctccGACGAAGCGAAAATAGCGATGGCGAGCGACGCAGATGCGGGGGATAGCGATACCCGCAACGGTCAGGCCCTGCTGGATCTGCAACGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGAGGCATAAGCGTCGTTaaagga                                                                                                                                   <  2:302879/148‑6 (MQ=255)
                  cACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgag                                                                                                                                   >  1:555278/1‑150 (MQ=255)
                  cACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgag                                                                                                                                   <  2:712618/150‑1 (MQ=255)
                  cACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGGCAGGCCCTGCTGGAGCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtgag                                                                                                                                   <  1:1649213/150‑1 (MQ=255)
                   aCCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGAGGCATAAGCGTCGTTaaaggat                                                                                                                                  <  2:1654124/150‑7 (MQ=255)
                        cgaAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGatat                                                                                                                             >  1:1292524/1‑150 (MQ=255)
                           aGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGAAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCGTCGTTAGTGAGTGATATCgg                                                                                                                          >  1:1820896/1‑150 (MQ=255)
                            gCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGt                                                                                                                         <  2:394/150‑1 (MQ=255)
                             cGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGta                                                                                                                        <  2:1885385/150‑1 (MQ=255)
                             cGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACAGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGta                                                                                                                        <  2:657243/150‑1 (MQ=255)
                              gAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaa                                                                                                                       >  2:449850/1‑150 (MQ=255)
                              gAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGAAAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACACTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaa                                                                                                                       <  1:1172954/150‑1 (MQ=255)
                              gAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGAGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTAGGCCGCGTTAGTGAGTGAGATCGGtaa                                                                                                                       >  1:852789/1‑150 (MQ=255)
                               aaaaTAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaat                                                                                                                      <  1:244735/150‑1 (MQ=255)
                                 aaTAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaataa                                                                                                                    <  2:1037808/150‑1 (MQ=255)
                                  aTAGCGATGTCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATaaa                                                                                                                   <  1:1204685/150‑1 (MQ=255)
                                  aTAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGCCGCTTATGCCTCGTTAGGGAGTGCTATCGGTAATaaa                                                                                                                   >  2:277158/1‑150 (MQ=255)
                                        aTGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACcgcg                                                                                                             <  2:196931/150‑1 (MQ=255)
                                        aTGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACcgcg                                                                                                             <  1:1640184/150‑1 (MQ=255)
                                          ggCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAATGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGAc                                                                                                           <  1:1432947/150‑1 (MQ=255)
                                            cgagcgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGt                                                                                                         >  1:1803663/1‑150 (MQ=255)
                                            cgagcgaAGAAGATGCGGGTGATAGCGATAAACGAAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGt                                                                                                         >  1:1228210/1‑150 (MQ=255)
                                              agcgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTTCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTg                                                                                                       <  2:1060945/150‑1 (MQ=255)
                                               gcgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGa                                                                                                      >  1:7253/1‑150 (MQ=255)
                                                cgaAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaa                                                                                                     >  1:1439856/1‑150 (MQ=255)
                                                 gaagaaGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaa                                                                                                    <  2:1758956/150‑1 (MQ=255)
                                                  aagaagATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTTCTGGATCTGCAAATCAACAGTAAAAAGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGAAATCGGTAATAAAACCGCGACGTTGaaaa                                                                                                   <  1:35078/150‑1 (MQ=255)
                                                    gaagaTGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAAcc                                                                                                 <  2:1052134/150‑1 (MQ=255)
                                                      agaTGCGGGTGATAGCGATAACCGGAACGGTCAGGACCTGCTGGATCTGCAAAGCAGCAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAg                                                                                               <  1:1795215/150‑1 (MQ=255)
                                                      agaTGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAg                                                                                               <  1:123290/150‑1 (MQ=255)
                                                         tGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAg                                                                                            >  1:776780/1‑150 (MQ=255)
                                                          gCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAgc                                                                                           <  2:1649556/150‑1 (MQ=255)
                                                           cGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTTCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAgcg                                                                                          <  1:1308497/150‑1 (MQ=255)
                                                           cGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAACCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGATAACCAGTAgcg                                                                                          >  2:1796206/1‑150 (MQ=255)
                                                           cGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAgcg                                                                                          >  2:1280750/1‑150 (MQ=255)
                                                                     agATACCCGCAACGGTCAGGCCCTGCTGGACCGGCAAAGCAACAGTAAAACGGTGGGCGGGGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATAGGTAATTAATCCCCGACGTTGAAAACCAGTATCGCCACGCAAgg                                                                                <  1:319867/149‑1 (MQ=255)
                                                                          aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                           <  2:1271973/150‑1 (MQ=255)
                                                                          aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                           >  1:723723/1‑150 (MQ=255)
                                                                           ccGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGGGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTTATAAAACCGCGACGTTGAAAACCAGTAGCGCCCCGCAAGGTTAtgt                                                                          >  1:308641/1‑150 (MQ=255)
                                                                              cAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATgtggt                                                                       <  2:1891621/150‑1 (MQ=255)
                                                                                 cGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGAc                                                                    >  1:1578583/1‑150 (MQ=255)
                                                                                 cGGTCAGGCCCTGCAGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCATCGTTAGTGAGTGATATCGGTAATAAAACAGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGAc                                                                    >  1:1695844/1‑150 (MQ=255)
                                                                                  ggTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACg                                                                   <  1:423481/150‑1 (MQ=255)
                                                                                     cAGGCCCTGCTGGATCTGCAAAGCATCAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGCCGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAg                                                                >  2:866459/1‑150 (MQ=255)
                                                                                           ctgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAACCCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTcc                                                          <  2:894180/150‑1 (MQ=255)
                                                                                           ctgctgGATCTGCAAAGAAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTcc                                                          <  1:960911/150‑1 (MQ=255)
                                                                                            tgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCa                                                         >  2:414882/1‑150 (MQ=255)
                                                                                            tgctgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTCAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCa                                                         >  2:1479683/1‑150 (MQ=255)
                                                                                               tgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATc                                                      <  1:483015/150‑1 (MQ=255)
                                                                                                ggATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGCTATCGGTTATACAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcc                                                     >  1:1730141/1‑149 (MQ=255)
                                                                                                  aTCTGCAATGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATAACGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagc                                                   <  1:1265635/150‑1 (MQ=255)
                                                                                                     tGCAAAGCAACAGTAACACGGTGGGCGGTGCGAAATCCTTACACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTACTAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGGGACGCAGCTTTCCAATcagcagc                                                <  2:1616430/150‑1 (MQ=255)
                                                                                                      gCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGGTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcagcc                                               >  2:941990/1‑149 (MQ=255)
                                                                                                       cAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTTACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcagcag                                              >  1:119537/1‑150 (MQ=255)
                                                                                                       cAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATcagcagcag                                              >  1:398992/1‑150 (MQ=255)
                                                                                                        aaaGCAACAGTAAAACGGTGGGCGGTGCGAGATCCTTTAACGACGCTTAGGCCTCGTTAGTGAGTGATATCGGTTATAAAACCGCGACGTTGAAGCCCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAAGccgcagcagg                                             >  2:3411/1‑149 (MQ=255)
                                                                                                         aaGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTc                                            >  2:118463/1‑150 (MQ=255)
                                                                                                          aGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCg                                           <  1:875456/150‑1 (MQ=255)
                                                                                                            cAACAGTAAAACGGTGGGCGTTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGAt                                         <  2:852789/150‑1 (MQ=255)
                                                                                                            cAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGAt                                         >  2:654392/1‑150 (MQ=255)
                                                                                                                 gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCg                                    <  1:1590859/150‑1 (MQ=255)
                                                                                                                 gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCg                                    >  1:640886/1‑150 (MQ=255)
                                                                                                                    aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                 <  2:1748931/150‑1 (MQ=255)
                                                                                                                    aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                 <  2:7096/150‑1 (MQ=255)
                                                                                                                         gTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAAt                            >  1:492214/1‑150 (MQ=255)
                                                                                                                            ggCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAAtctc                         <  1:1520821/150‑1 (MQ=255)
                                                                                                                            ggCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAAtctc                         <  1:454366/150‑1 (MQ=255)
                                                                                                                             gCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATCTCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCg                        <  2:963187/150‑1 (MQ=255)
                                                                                                                              cGGTGCGAAATCCTTTAACGACGCTTATGCCTCGGTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGa                       <  2:385400/150‑1 (MQ=255)
                                                                                                                                gTGCGAAATCCTTTAACGCCGCTTATGTCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATg                     <  1:754494/150‑1 (MQ=255)
                                                                                                                                gTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATg                     <  1:1737224/150‑1 (MQ=255)
                                                                                                                                   cGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgagg                  <  2:1093046/150‑1 (MQ=255)
                                                                                                                                      aaTCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGt               >  1:841061/1‑150 (MQ=255)
                                                                                                                                       aTCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTa              >  1:942487/1‑150 (MQ=255)
                                                                                                                                         ccTTTAACGACGCTTATGCCTCGTTAGTGAGTGGTATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACg            >  1:1745875/1‑150 (MQ=255)
                                                                                                                                          cTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACgg           <  2:203670/150‑1 (MQ=255)
                                                                                                                                          atttAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACgg           >  1:7945/2‑150 (MQ=255)
                                                                                                                                           tttAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGa          >  2:1560887/1‑150 (MQ=255)
                                                                                                                                            ctAACGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGaa         >  1:604407/2‑150 (MQ=255)
                                                                                                                                             tAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGaaa        >  2:147680/1‑150 (MQ=255)
                                                                                                                                              aaCGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAAt       >  2:1582566/1‑150 (MQ=255)
                                                                                                                                              aaCGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAAt       >  1:86399/1‑150 (MQ=255)
                                                                                                                                                 gacgCTTATGCCTCGTTAGTGAGGGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTg    <  1:942784/150‑1 (MQ=255)
                                                                                                                                                  acgCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGc   >  2:1035420/1‑150 (MQ=255)
                                                                                                                                                   cgCTTATGCCTCGTTAGTGAGTGATAGCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCa  <  2:1468556/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCA  >  NC_000913/1139640‑1139936

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: