Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,935,841 T→G 5.2% Q192P (CAG→CCG)  flhF ← flagellar biosynthesis protein FlhF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,935,8410TG 5.2% 694.1 / 4.3 248Q192P (CAG→CCG) flhFflagellar biosynthesis protein FlhF
Reads supporting (aligned to +/- strand):  ref base T (122/106);  new base G (2/17);  total (125/123)
Fisher's exact test for biased strand distribution p-value = 2.36e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

ACGTAGCGCGCGGCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGAGTTGCACT  >  NC_002947/4935694‑4935990
                                                                                                                                                   |                                                                                                                                                     
aCGTAGCGCGCGGCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTgg                                                                                                                                                     >  2:273766/1‑150 (MQ=255)
  gTAGCGCGCGGCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGgcg                                                                                                                                                   >  2:338331/1‑150 (MQ=255)
    aGCGCGCGGCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGcgcg                                                                                                                                                 <  1:1622072/150‑1 (MQ=255)
       gcgcgGCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGgcgcat                                                                                                                                                >  2:2863870/1‑146 (MQ=255)
           ggCCAGCTTGGCCTGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCTTGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCACGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTc                                                                                                                                          <  1:100857/150‑1 (MQ=255)
           ggCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTc                                                                                                                                          <  2:3603241/150‑1 (MQ=255)
           ggCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTc                                                                                                                                          >  1:4379925/1‑150 (MQ=255)
           ggCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCATCTTCGATCGGCTCTATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTc                                                                                                                                          <  1:100794/150‑1 (MQ=255)
            gCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCg                                                                                                                                         >  1:4205199/1‑150 (MQ=255)
              cAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCTTGCGCCAGGCCTGGCGCGGCTCTTCGAt                                                                                                                                       >  2:314033/1‑150 (MQ=255)
              cAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGAt                                                                                                                                       <  2:1422796/150‑1 (MQ=255)
               aGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGAtc                                                                                                                                      >  2:3255697/1‑150 (MQ=255)
               aGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGAtc                                                                                                                                      >  1:2643669/1‑150 (MQ=255)
                gCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGAtct                                                                                                                                     >  2:4005819/1‑150 (MQ=255)
                   tGGCCTGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGAATGCGCCAGGCCTGGCGCGGCTCTTCGATCTcgg                                                                                                                                  <  2:1515644/150‑1 (MQ=255)
                         gggTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCa                                                                                                                            <  2:2671806/150‑1 (MQ=255)
                           gtggtggtCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCTCCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAgg                                                                                                                          <  2:4046868/150‑1 (MQ=255)
                               tggtCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCca                                                                                                                      >  1:966905/1‑150 (MQ=255)
                               tggtCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCca                                                                                                                      >  1:77610/1‑150 (MQ=255)
                                ggtCTGGCCCTGGCCCGCCGGGCCGACCTGGGCGAACCCCCCACCCTCTTCGATGGGCTCTATTTCCGGCACTTCGATCATCCGTGCCAAGTGCGCCAGCAGCATGCGCCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCcag                                                                                                                     <  1:2001082/150‑1 (MQ=255)
                                 gtCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCcagc                                                                                                                    >  2:2107026/1‑150 (MQ=255)
                                  tCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCcagca                                                                                                                   >  2:2377276/1‑150 (MQ=255)
                                   cTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCcagcag                                                                                                                  >  2:3825989/1‑150 (MQ=255)
                                        ccATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTcgc                                                                                                             >  2:4398744/1‑150 (MQ=255)
                                        ccATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTcgc                                                                                                             >  2:680658/1‑150 (MQ=255)
                                          aTGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCgg                                                                                                           >  2:556430/1‑150 (MQ=255)
                                             cccGCCGGGCCGACCATGGCGATCACCCCACCCTCTGCGATCGGCTCGATTGCAGGCACTTCGATCATCCGTGCCAAGTGAGCCCGCAGCATGCGCCAGGCCGGGCGCGGCGCTTCGAGCTCGGCGGGCAGGTCCAGCAGCTCGCGGGCg                                                                                                        >  1:3012810/1‑150 (MQ=255)
                                              ccgccgGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGa                                                                                                       >  1:2399980/1‑150 (MQ=255)
                                              ccgccgGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGa                                                                                                       <  2:2025553/150‑1 (MQ=255)
                                              ccgccgGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGa                                                                                                       >  2:3771300/1‑150 (MQ=255)
                                                gccgGGCCTCCCTTGGCTATCACCCCACCCTCTTCGATCGGATCGATTACAGGCACTTCGAACATCCGTCCCAAGTGCGCCAGCAGCTGGCGCCCGGCCGGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATc                                                                                                     <  1:1233635/150‑1 (MQ=255)
                                                gccgGGCCGACCAGGGCGATCACCCCACGCTGTTCGATGGGCTCTATTTCCGGCTTTTGTATCATCCGTGCCTTGTGACCCAGCAGCATGCGCCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATc                                                                                                     <  1:4226931/150‑1 (MQ=255)
                                                    ggCCGACCCTGGCGACCACCCCACCCGCTTCGCTCGGCTCGATTTCCGGCACTTCGATCATCCGTGCAATGTGAGCCTGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGcc                                                                                                 <  1:805570/150‑1 (MQ=255)
                                                        gACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGa                                                                                             >  1:4061184/1‑150 (MQ=255)
                                                           cATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAg                                                                                          >  1:363197/1‑150 (MQ=255)
                                                           cATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAg                                                                                          >  1:1867926/1‑150 (MQ=255)
                                                             tGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCCGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGc                                                                                        >  2:3614990/1‑150 (MQ=255)
                                                             tGGCGACCCCCCCACCCACTTCTTTCGGCTTTTTTTCCGCCAATTCTATCATCGTGGCCACGTGAGCCCGCAGCCTGCGCCAGGCCGGCCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGc                                                                                        <  1:1130107/150‑1 (MQ=255)
                                                                 gATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCATGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGAt                                                                                    <  2:1219546/150‑1 (MQ=255)
                                                                    caccccacccTCTTCGATCGGCTCGATTTCCGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCGGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATgcg                                                                                 <  1:4452611/150‑1 (MQ=255)
                                                                      ccccacccTCTTCGATCGGCTCGTTTTCCGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCt                                                                               <  1:4004565/150‑1 (MQ=255)
                                                                      ccccacccTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGACCGGACAGGCCGATGCGCt                                                                               >  1:794494/1‑150 (MQ=255)
                                                                         cacccTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCACGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCa                                                                            >  1:4417630/1‑150 (MQ=255)
                                                                          acccTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAg                                                                           >  2:6336/1‑150 (MQ=255)
                                                                          acccTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAg                                                                           >  1:3753209/1‑150 (MQ=255)
                                                                            ccTCTTCGATCGGCTCGATTTCAGGCACTTCGATCACCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAggc                                                                         <  2:772339/150‑1 (MQ=255)
                                                                                ttCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGgcggcg                                                                     >  2:2644555/1‑150 (MQ=255)
                                                                                ttCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGgcggcg                                                                     >  1:3241603/1‑150 (MQ=255)
                                                                                ttCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGgcggcg                                                                     >  2:1672632/1‑150 (MQ=255)
                                                                                     tCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCaaa                                                                >  2:2605497/1‑150 (MQ=255)
                                                                                     gcGGCTCGTTTTCCGGCACTTCGACCCCCCGGGCCACGTGCGCCCGCCGCAGGCGCCCGGCCGGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCaaa                                                                <  1:1292297/149‑1 (MQ=255)
                                                                                         ctctttttCAGGAACTTCTATCTTCCGTGCAATGTGAGCCAGGAGCTTGCGCCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGtt                                                            <  1:2463684/145‑1 (MQ=255)
                                                                                         cTCGATTTCCGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGtt                                                            <  1:2982884/150‑1 (MQ=255)
                                                                                         cTCGATTTCAGCCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGtt                                                            <  2:3571691/150‑1 (MQ=255)
                                                                                          tCGATTTCCGCCACTTCGATCATCCGGGCCAAGTGAGCCAGCCGCATGCGCCCGGCCGGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTg                                                           <  1:600738/150‑1 (MQ=255)
                                                                                          tCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTg                                                           >  1:990396/1‑150 (MQ=255)
                                                                                          tCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTg                                                           >  1:2779011/1‑150 (MQ=255)
                                                                                          tCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTACAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTg                                                           >  1:4478868/1‑150 (MQ=255)
                                                                                           ctaTTTCAGGCACTTCGTTCATCCGTGCCAAGTGAGCCCGCAGCTTGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTgg                                                          <  2:1892472/148‑1 (MQ=255)
                                                                                            gATTCCTGGCACTTCGATCATCCGTGCCATGTGACCCAGCAGCAGGCGCCCGGCCGGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGc                                                         <  1:1345161/150‑1 (MQ=255)
                                                                                                 ccggCTCTTCTATCATCCGTGACATGTGCGCCAGCAGCTGGCGCGCGGCCGGCCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCctgct                                                    <  1:893102/148‑1 (MQ=255)
                                                                                                  aGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCctgctg                                                   >  1:3528412/1‑150 (MQ=255)
                                                                                                  aGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCctgctg                                                   >  1:3711961/1‑150 (MQ=255)
                                                                                                    gCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTAGATCTCGGCGGTCAGGTCCAGCAGCACGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCCAAGGTTGGCCTgcggcg                                                 >  2:3447195/1‑150 (MQ=255)
                                                                                                    gCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGACTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCg                                                 <  1:2831340/150‑1 (MQ=255)
                                                                                                        ttCGATCATCCGTGCCATGTGAGCCAGCAGCATGCGCCCGGCCGGCCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCtt                                             <  1:2562725/150‑1 (MQ=255)
                                                                                                        ttCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCtt                                             >  1:220005/1‑150 (MQ=255)
                                                                                                        ttCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCtt                                             >  2:1312335/1‑150 (MQ=255)
                                                                                                         tCGATCATCCGTGCCAAGTGAGCCAGCCGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTtg                                            <  2:2703666/150‑1 (MQ=255)
                                                                                                         tCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTtg                                            <  2:853439/150‑1 (MQ=255)
                                                                                                              catcCGTGCCAAGTGAACCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGcc                                       <  2:1231393/150‑1 (MQ=255)
                                                                                                                tcCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCt                                     >  2:2360362/1‑150 (MQ=255)
                                                                                                                   gTGCCATGTGCCCCACCAGCCGGCGCCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCa                                  <  1:1921380/150‑1 (MQ=255)
                                                                                                                   gTGCCATGTGAGCCAGCAGCATGCGCCTGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCa                                  <  2:2436764/150‑1 (MQ=255)
                                                                                                                     gCCAAGTGAGCCAGCCGCATGCGCCCGCCCGGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGc                                <  1:2283578/150‑1 (MQ=255)
                                                                                                                        atgTGAGCCAGCAGCAGGCGCCCGGCCGGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGa                             <  1:4174173/148‑1 (MQ=255)
                                                                                                                        aaGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGa                             <  2:3174205/150‑1 (MQ=255)
                                                                                                                        aaGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGa                             <  1:680658/150‑1 (MQ=255)
                                                                                                                        aaGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTACTGCGGCTTGCTGCCCTGCAGCTGa                             <  2:3554786/150‑1 (MQ=255)
                                                                                                                         aGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGAt                            >  2:1050221/1‑150 (MQ=255)
                                                                                                                          gTGCGCCAGCAGCAGGCGCCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGAtt                           <  1:1711425/150‑1 (MQ=255)
                                                                                                                           tGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTc                          <  2:4205199/150‑1 (MQ=255)
                                                                                                                           tGAGCCAGCAGCAGGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTc                          <  2:3880564/150‑1 (MQ=255)
                                                                                                                             cgCCAGCAGCATGCGCCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCa                        <  2:1855556/149‑1 (MQ=255)
                                                                                                                                cagcagCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCCGGCGGCGCCACAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGc                     >  2:893472/1‑150 (MQ=255)
                                                                                                                                cagcagCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGc                     <  2:3317211/150‑1 (MQ=255)
                                                                                                                                  gcatcaTGCGCCAGCCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCa                   <  2:4478868/150‑1 (MQ=255)
                                                                                                                                  gcagcaTGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCa                   >  1:340854/1‑150 (MQ=255)
                                                                                                                                  gcagcaGGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCa                   <  1:3665225/150‑1 (MQ=255)
                                                                                                                                  ccagccTGCGCCAGGCCTGCCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCa                   <  1:2176964/149‑1 (MQ=255)
                                                                                                                                    agccggcgcCCGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTACAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGc                 <  1:2321020/145‑1 (MQ=255)
                                                                                                                                       ccgcgcAAGGCCTGGCGCGGCTCTTCGATCTCGGAGGTCAGGTACAGCATCTCGCGGGCGATCGGCCCGGACCGGCCGATGCGCTGCAGGCGGCGCCCAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCcgcgcc              >  2:1095705/3‑148 (MQ=255)
                                                                                                                                       aTGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCGGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCgca              >  1:1601743/1‑148 (MQ=255)
                                                                                                                                         gcgcCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAAcc            <  2:4385186/150‑1 (MQ=255)
                                                                                                                                           gcCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCCAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGa          >  2:12602/1‑150 (MQ=255)
                                                                                                                                           gcCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGa          >  2:1753646/1‑150 (MQ=255)
                                                                                                                                               ggCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGAGTTg      <  2:751582/150‑1 (MQ=255)
                                                                                                                                               gcccTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGAGTTg      <  1:4291926/148‑1 (MQ=255)
                                                                                                                                                  cTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGAGTTGCAc   >  1:3316598/1‑150 (MQ=255)
                                                                                                                                                   tGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGAGTTGCACt  >  2:2967837/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACGTAGCGCGCGGCCAGCTTGGCCAGGGTGGTGGTCTTGCCCATGCCCGCCGGGCCGACCATGGCGATCACCCCACCCTCTTCGATCGGCTCGATTTCAGGCACTTCGATCATCCGTGCCAAGTGAGCCAGCAGCATGCGCCAGGCCTGGCGCGGCTCTTCGATCTCGGCGGTCAGGTCCAGCAGCTCGCGGGCGATCGGCCCGGACAGGCCGATGCGCTGCAGGCGGCGCCAAAGGTTGGCCTGCTGCGGCTTGCTGCCCTGCAGCTGATTCCAGGCCAGCGAACCGAGTTGCACT  >  NC_002947/4935694‑4935990

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