Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,139,796 C→G 5.8% A473A (GCC→GCG 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,7960CG 5.8% 358.9 / 10.8 119A473A (GCC→GCGflgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base C (48/63);  new base G (4/3);  total (53/66)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.03e-01

Only 100 of 121 total aligned reads displayed.

TATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGTTTTCA  >  NC_000913/1139651‑1139945
                                                                                                                                                 |                                                                                                                                                     
tATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGt                                                                                                                                                   >  1:1748931/1‑150 (MQ=255)
   tAATCACCGACGAAGCGAAAATAGCGATGGCTAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTTGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTag                                                                                                                                                <  2:1289924/150‑1 (MQ=255)
     aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                              <  2:1381870/150‑1 (MQ=255)
     aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                              >  2:1520821/1‑150 (MQ=255)
     aTCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTagtg                                                                                                                                              <  1:1470669/150‑1 (MQ=255)
      tCACCGACGAAGCGAAAATAGAGATGGCTAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGAGGCATAAGCGTCGTTaaagg                                                                                                                                             <  1:1542993/150‑5 (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)
                 gCGAATATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCTACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGt                                                                                                                                  <  1:1223328/150‑1 (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)
                    aaaaTAGCGATGGCGAGCGATTAAGATTCGGGTTATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtaat                                                                                                                               <  1:1693383/150‑1 (MQ=255)
                    aaaaTAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGtata                                                                                                                               <  2:1887318/150‑3 (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)
                           cGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGAAAAGCAACAGTAAATCTGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACcg                                                                                                                        <  1:245611/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                                                                                                              <  2:1690203/150‑1 (MQ=255)
                                      gaagaaGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGaaa                                                                                                             <  2:1758956/150‑1 (MQ=255)
                                       aagaagATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTTCTGGATCTGCAAATCAACAGTAAAAAGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGAAATCGGTAATAAAACCGCGACGTTGaaaa                                                                                                            <  1:35078/150‑1 (MQ=255)
                                        agaagaTGCGGGTGATAGCGATAACCGCAACGGTCATGCCCTGCTGGATCTGCAAAGCAACAGTAACACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATAGCGGTACTAAACCCGCGACGGTGAacac                                                                                                           <  2:1370787/150‑1 (MQ=255)
                                         gaagaTGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAAcc                                                                                                          <  2:1052134/150‑1 (MQ=255)
                                           agaTGCGGGTGATAGCGATAACCGGAACGGTCAGGACCTGCTGGATCTGCAAAGCAGCAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAg                                                                                                        <  1:1795215/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)
                                                cGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAgcg                                                                                                   >  2:1280750/1‑150 (MQ=255)
                                                   gTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCa                                                                                                <  2:1555383/150‑1 (MQ=255)
                                                     gATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACg                                                                                              <  1:579772/150‑1 (MQ=255)
                                                          agATACCCGCAACGGTCAGGCCCTGCTGGACCGGCAAAGCAACAGTAAAACGGTGGGCGGGGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATAGGTAATTAATCCCCGACGTTGAAAACCAGTATCGCCACGCAAgg                                                                                         <  1:319867/149‑1 (MQ=255)
                                                               aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                    >  1:723723/1‑150 (MQ=255)
                                                               aCCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtg                                                                                    <  2:1271973/150‑1 (MQ=255)
                                                                ccGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAAtgt                                                                                   >  1:1422512/1‑150 (MQ=255)
                                                                ccGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGGGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTTATAAAACCGCGACGTTGAAAACCAGTAGCGCCCCGCAAGGTTAtgt                                                                                   >  1:308641/1‑150 (MQ=255)
                                                                  gCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATgtgg                                                                                 <  2:1104327/150‑1 (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)
                                                                                ctgctgGATCTGCAAAGCTACAGTAAAACGGTGGGCGGTGCTAAATCCTTTAACGACGCTAATGCCTCGTTAGTGAGTGCAATCGGTAAAAAAACCGCTAAGTTGAAAACCAGTACCGACACGCCATCTAATGATTTGAAGAAGCTTTcc                                                                   >  2:336151/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)
                                                                                    tgGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATc                                                               >  1:981970/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)
                                                                                               aGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCg                                                    <  1:875456/150‑1 (MQ=255)
                                                                                                 cAACAGTAAAACGGTGGGCGTTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGAt                                                  <  2:852789/150‑1 (MQ=255)
                                                                                                 cAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGAt                                                  >  2:654392/1‑150 (MQ=255)
                                                                                                     aGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTcc                                              <  1:1528553/150‑1 (MQ=255)
                                                                                                      gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCg                                             <  1:1590859/150‑1 (MQ=255)
                                                                                                      gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCg                                             >  2:790598/1‑150 (MQ=255)
                                                                                                      gTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCg                                             >  1:640886/1‑150 (MQ=255)
                                                                                                         aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                          <  2:7096/150‑1 (MQ=255)
                                                                                                         aaaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtg                                          <  2:1748931/150‑1 (MQ=255)
                                                                                                          aaCGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGgtgt                                         >  1:862564/1‑150 (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)
                                                                                                                     gTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATg                              <  1:1737224/150‑1 (MQ=255)
                                                                                                                        cGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgagg                           <  2:1093046/150‑1 (MQ=255)
                                                                                                                          aaaTCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATgaggag                         <  1:746993/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)
                                                                                                                                  tAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGaaa                 >  2:147680/1‑150 (MQ=255)
                                                                                                                                   aaCGAGGCATAAGCGTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAAt                >  2:1582566/1‑150 (MQ=255)
                                                                                                                                      gacgCTTATGCCTCGTTAGTGAGGGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTg             <  1:942784/150‑1 (MQ=255)
                                                                                                                                        cgCTTATGCCTCGTTAGTGAGTGATAGCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCa           <  2:1468556/150‑1 (MQ=255)
                                                                                                                                          cTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAAc         <  2:1198482/150‑1 (MQ=255)
                                                                                                                                          cTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAAc         >  2:914366/1‑150 (MQ=255)
                                                                                                                                          cTTATGCCTCGTTAGGGAGTGAGATCGGTATTAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCTAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGCTGTCAATATCGATGAGGAGGACGGAAATCTGCacc         >  2:129312/1‑148 (MQ=255)
                                                                                                                                              tGCCTCGTTAGTGAGAGATATCGGTAATAAACCCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCATCAGCTGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGttt     <  2:180030/150‑1 (MQ=255)
                                                                                                                                                 cTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGTTTTca  >  2:812218/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCTGCTGGATCTGCAAAGCAACAGTAAAACGGTGGGCGGTGCGAAATCCTTTAACGACGCTTATGCCTCGTTAGTGAGTGATATCGGTAATAAAACCGCGACGTTGAAAACCAGTAGCGCCACGCAAGGTAATGTGGTGACGCAGCTTTCCAATCAGCAGCAGTCGATTTCCGGTGTCAATCTCGATGAGGAGTACGGAAATCTGCAACGTTTTCA  >  NC_000913/1139651‑1139945

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: