Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,400,840 C→A 5.1% L268L (CTC→CTA ybgF → TolQRA transport system periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,400,8400CA 5.1% 257.2 / 2.4 78L268L (CTC→CTAybgFTolQRA transport system periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (35/39);  new base A (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GCAAGGTGCCAGATTCGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTT  >  NC_002947/1400692‑1400988
                                                                                                                                                    |                                                                                                                                                    
gCAAGGTGCCAGATTCGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGctct                                                                                                                                                     <  1:655266/150‑1 (MQ=255)
    ggTGCCAGATTCGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGc                                                                                                                                                 <  1:734078/150‑1 (MQ=255)
    ggTGCCAGATTCGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGc                                                                                                                                                 <  2:358493/150‑1 (MQ=255)
            aTTCGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCCGTGACCGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGc                                                                                                                                         <  1:1059777/150‑1 (MQ=255)
             ttCGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCt                                                                                                                                        <  1:758842/150‑1 (MQ=255)
             tgcGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCt                                                                                                                                        <  2:238919/148‑1 (MQ=255)
               cGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGt                                                                                                                                      <  2:293740/150‑1 (MQ=255)
                   gTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCg                                                                                                                                  >  1:970736/1‑150 (MQ=255)
                          cTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCgaaa                                                                                                                           <  1:1302895/150‑1 (MQ=255)
                          cTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCgaaa                                                                                                                           <  2:367738/150‑1 (MQ=255)
                          cTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCgaaa                                                                                                                           <  2:568565/150‑1 (MQ=255)
                               tGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCTCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAAc                                                                                                                      >  1:921474/1‑150 (MQ=255)
                               tGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAAc                                                                                                                      <  1:316566/150‑1 (MQ=255)
                                   gTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCcgc                                                                                                                  >  2:886977/1‑150 (MQ=255)
                                              tATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCgg                                                                                                       >  2:779341/1‑150 (MQ=255)
                                                    ccATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGttttt                                                                                                 >  1:492597/1‑150 (MQ=255)
                                                    ccATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTATAGGCGTTGTCg                                                                                                                                          <  1:1800402/109‑1 (MQ=255)
                                                    ccATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTATAGGCGTTGTCg                                                                                                                                          >  2:1800402/1‑109 (MQ=255)
                                                    ccATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTATAGGCGTTGTCg                                                                                                                                          <  1:1800923/109‑1 (MQ=255)
                                                    ccATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTATAGGCGTTGTCg                                                                                                                                          >  2:1800923/1‑109 (MQ=255)
                                                      aTACCGACAACGTCAACGGCATCCTTGAGCAGGTCGTGACCCAGTCCCCCGGCACCTCTGCCGCTCAACTGGCCCCGCGTGGCCTGCAGAAGGTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTc                                                                                               <  1:858123/150‑1 (MQ=255)
                                                       tACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGACTGGCGCGGGtttctttc                                                                                              <  2:1614752/150‑3 (MQ=255)
                                                          cGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTa                                                                                           >  2:1908415/1‑150 (MQ=255)
                                                          cGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTa                                                                                           >  1:879841/1‑150 (MQ=255)
                                                           gACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAg                                                                                          >  2:757518/1‑150 (MQ=255)
                                                             cAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGaa                                                                                        <  1:63424/150‑1 (MQ=255)
                                                                ggTAAAAGGCATCCTTCAGCAGGTCGTTACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAACGAAACCCGCGCCATTCGCGGGTTTTTTCGTTAGAATCa                                                                                     >  1:397911/1‑150 (MQ=255)
                                                                ggTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCTCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCa                                                                                     <  1:38058/150‑1 (MQ=255)
                                                                ggTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCa                                                                                     >  1:1310574/1‑150 (MQ=255)
                                                                 gTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGt                                                                                             >  2:1141527/1‑141 (MQ=255)
                                                                 gTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGt                                                                                             <  1:1141527/141‑1 (MQ=255)
                                                                  tAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTAGTTTTGCCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACt                                                                                   <  2:1910699/150‑1 (MQ=255)
                                                                    aaaGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAGCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGc                                                                                 <  2:1177445/150‑1 (MQ=255)
                                                                          aTCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCtttt                                                                           >  1:1738779/1‑150 (MQ=255)
                                                                           tCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCATCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCTGTCGCGGGTTTTTTCGTTAGAATCACTGCCCttttt                                                                          <  2:392258/150‑1 (MQ=255)
                                                                           tCCGTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCAAGAATGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCttttt                                                                          >  2:1206528/1‑150 (MQ=255)
                                                                            ccTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAAACACTGCCCtttttt                                                                         <  1:1908415/150‑1 (MQ=255)
                                                                             tttCAGCAGGACGTGACCAAGTTCAACGGCACCTCTGCCGATCAACTGGCCCAGCGTGAACTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTc                                                                        <  2:397911/149‑1 (MQ=255)
                                                                             cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTc                                                                        <  2:970736/150‑1 (MQ=255)
                                                                             cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTc                                                                        >  1:226785/1‑150 (MQ=255)
                                                                             cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTc                                                                        >  1:1765702/1‑150 (MQ=255)
                                                                             cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTc                                                                        <  1:1816741/150‑1 (MQ=255)
                                                                             cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAAAACCCGCGACTGGCGCGGGtttctttc                                                                                              >  2:253552/1‑122 (MQ=38)
                                                                             cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAAAACCCGCGACTGGCGCGGGtttctttc                                                                                              <  1:253552/127‑6 (MQ=38)
                                                                                  gcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTaaa                                                                   <  1:322927/150‑1 (MQ=255)
                                                                                      gTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCAAGGATGGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACg                                                               >  1:987632/1‑150 (MQ=255)
                                                                                        cGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACgct                                                             <  2:1765702/150‑1 (MQ=255)
                                                                                        cGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACgct                                                             <  1:1576885/150‑1 (MQ=255)
                                                                                             cccAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCtgctgc                                                        >  1:1397575/1‑150 (MQ=255)
                                                                                                 gTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGAt                                                    >  2:1285247/1‑150 (MQ=255)
                                                                                                       cGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCt                                              <  1:113336/150‑1 (MQ=255)
                                                                                                       cGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGGTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCt                                              >  1:1102955/1‑150 (MQ=255)
                                                                                                              tctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTGCGGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAACCACGCTGCTGCGGATAACGCTATGTCGa                                       >  1:118882/1‑150 (MQ=255)
                                                                                                              tctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                       <  2:1351193/150‑1 (MQ=255)
                                                                                                              tctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGATAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                       <  2:492597/150‑1 (MQ=255)
                                                                                                                 gCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTc                                                                                                           >  2:648680/1‑79 (MQ=255)
                                                                                                                 gCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTc                                                                                                           <  1:648680/79‑1 (MQ=255)
                                                                                                                    gCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcgc                                 >  1:1223131/1‑150 (MQ=255)
                                                                                                                      tCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGa                               >  1:6310/1‑150 (MQ=255)
                                                                                                                        aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTa                                                                     >  1:809391/1‑110 (MQ=255)
                                                                                                                        aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTa                                                                     <  2:809391/110‑1 (MQ=255)
                                                                                                                        aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTAGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACa                             >  2:1215819/1‑150 (MQ=255)
                                                                                                                         aCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCg                                                                                              >  2:1308854/1‑84 (MQ=255)
                                                                                                                         aCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCg                                                                                              <  1:1308854/84‑1 (MQ=255)
                                                                                                                          cTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACCATATGCATGGCCCGCGACAGt                           >  2:1789503/1‑150 (MQ=255)
                                                                                                                           tGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTg                          <  2:921474/150‑1 (MQ=255)
                                                                                                                               ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                       <  1:1344686/133‑1 (MQ=255)
                                                                                                                               ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                       >  2:1344686/1‑133 (MQ=255)
                                                                                                                                cAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGa                     <  2:879841/150‑1 (MQ=255)
                                                                                                                                    gTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGa                 <  1:1137584/150‑1 (MQ=255)
                                                                                                                                         cTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGTTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACTCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCgg            >  2:1073638/1‑150 (MQ=255)
                                                                                                                                          tGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTCTGGTAGTTTCGTAAACACGCTGCGGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGa           >  1:974189/1‑150 (MQ=255)
                                                                                                                                           gCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCtttt                                                                           >  2:779692/1‑85 (MQ=255)
                                                                                                                                           gCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCtttt                                                                           <  1:779692/85‑1 (MQ=255)
                                                                                                                                           gCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGAc          >  1:1371007/1‑150 (MQ=255)
                                                                                                                                            cAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACa         <  1:1206528/150‑1 (MQ=255)
                                                                                                                                             aGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGa                                                     >  2:240141/1‑105 (MQ=255)
                                                                                                                                             aGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGa                                                     <  1:240141/105‑1 (MQ=255)
                                                                                                                                              gAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGc       >  2:665551/1‑150 (MQ=255)
                                                                                                                                                  cgctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGt   >  1:1715671/3‑150 (MQ=255)
                                                                                                                                                   tctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt  <  2:817592/150‑1 (MQ=255)
                                                                                                                                                   tctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt  <  1:757518/150‑1 (MQ=255)
                                                                                                                                                    ctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCtttt                                                                           <  1:1221248/76‑1 (MQ=255)
                                                                                                                                                    ctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCtttt                                                                           >  2:1221248/1‑76 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAAGGTGCCAGATTCGCTGTACAAACTGGCTGATGTCGAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTT  >  NC_002947/1400692‑1400988

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