Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 914,647:1 +C 6.2% coding (1355/1743 nt) fruA → PTS fructose transporter subunit IIBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947914,6471.C 6.2% 343.5 / 15.3 81coding (1355/1743 nt)fruAPTS fructose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base . (37/39);  new base C (3/2);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

AATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGACCCGTTACGGG  >  NC_002947/914499‑914796
                                                                                                                                                       |                                                                                                                                                     
aaTGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGCCGGCGGCCATGGtc                                                                                                                                                       <  1:1319622/150‑2 (MQ=255)
  tGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgcc                                                                                                                                                    >  1:1628976/1‑150 (MQ=255)
  tGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgcc                                                                                                                                                    >  2:80444/1‑150 (MQ=255)
  tGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgcc                                                                                                                                                    <  1:521205/150‑1 (MQ=255)
  tGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgcc                                                                                                                                                    >  1:1393019/1‑150 (MQ=255)
    ccATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGCCGGCGGCCATGGtcgccg                                                                                                                                                   <  2:192859/150‑6 (MQ=255)
     cATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgccgcc                                                                                                                                                 >  1:410361/1‑150 (MQ=255)
       tCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGa                                                                                                                                               <  1:470723/150‑1 (MQ=255)
       tCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGa                                                                                                                                               <  1:1353318/150‑1 (MQ=255)
       tCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGa                                                                                                                                               <  2:1628976/150‑1 (MQ=255)
               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgccgc                                                                                                                                                  >  2:806510/1‑139 (MQ=255)
               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgccgc                                                                                                                                                  <  1:806510/139‑1 (MQ=255)
               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTg                                                                                                                                       <  1:1426219/150‑1 (MQ=255)
               gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGGTCTCGGTGGGCGTGCTGGCCTCATCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCgattggattg                                                                                                                                       >  2:710068/1‑150 (MQ=255)
                    gTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTAt                                                                                                                                  >  1:1213802/1‑150 (MQ=255)
                       gCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGc                                                                                                                               >  1:874575/1‑150 (MQ=255)
                           ggcggcATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTGCGCGCCAATGGGGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCAAcc                                                                                                                           >  1:606013/1‑150 (MQ=255)
                           ggcggcATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTgg                                                                                                                                      <  2:842313/139‑1 (MQ=255)
                           ggcggcATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTgg                                                                                                                                      >  1:842313/1‑139 (MQ=255)
                              ggCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTC                                                                                                                        >  1:322616‑M1/1‑150 (MQ=255)
                              ggCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTC                                                                                                                        <  1:852563‑M1/150‑1 (MQ=255)
                              ggCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTC                                                                                                                        <  2:1393019‑M1/150‑1 (MQ=255)
                                 atgatgTGCGTCGACCTTGGCGGGCCGATCAACGAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTg                                                                                                                     <  1:68720/150‑1 (MQ=255)
                                     tgtgCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGAt                                                                                                                                              >  1:1490521/1‑121 (MQ=255)
                                     tgtgCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGAt                                                                                                                                              <  2:1490521/121‑1 (MQ=255)
                                          gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGTGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTgg                                                                                                                                      <  1:508352/124‑1 (MQ=255)
                                          gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGTGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTgg                                                                                                                                      >  2:508352/1‑124 (MQ=255)
                                           tCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCTCGCCGGTGCTCCCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGGATTGCCACCTTCCTGGCCCGGCGCa                                                                                                           >  1:1790690/1‑150 (MQ=255)
                                            cGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCaa                                                                                                          <  2:572316/150‑1 (MQ=255)
                                            cGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCaa                                                                                                          <  2:410361/150‑1 (MQ=255)
                                             gACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAg                                                                                                         <  2:1221364/150‑1 (MQ=255)
                                                  tGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGc                                                                                                    <  1:1740615/150‑1 (MQ=255)
                                                        gCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCagag                                                                                              >  2:1366384/1‑150 (MQ=255)
                                                          cGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcg                                                                                            >  1:755573/1‑150 (MQ=255)
                                                           gATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcga                                                                                           >  1:102886/1‑150 (MQ=255)
                                                            aTCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgag                                                                                          >  2:411968/1‑150 (MQ=255)
                                                            aTCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgag                                                                                          >  1:770695/1‑150 (MQ=255)
                                                            aTCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgag                                                                                          <  1:80444/150‑1 (MQ=255)
                                                               aacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                                       <  1:1060930/150‑1 (MQ=255)
                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATCGGCGGCACCATGCCGGCGGCCATGGTCGCCGCCATTGGcgcgta                                                                                                                                     >  1:1072742/1‑98 (MQ=9)
                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATCGGCGGCACCATGCCGGCGGCCATGGTCGCCGCCATTGGcgcgta                                                                                                                                     <  2:1072742/103‑6 (MQ=9)
                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTtgcccgcctcgcgctcgctctgggcgaacttg                                                                                      >  1:1776316‑M1/1‑118 (MQ=255)
                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                                      >  1:1763337/1‑150 (MQ=255)
                                                                acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGcgcg                                                                                      >  2:1110893/1‑150 (MQ=255)
                                                                  aaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAg                                                                                    >  2:1310228/1‑150 (MQ=255)
                                                                   aGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAgg                                                                                   >  2:1644200/1‑150 (MQ=255)
                                                                     gccgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGTATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGCGCGCGAGGCg                                                                                 >  1:1539506/1‑150 (MQ=255)
                                                                      ccgccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTtgcccgcctcgcgctcgctctgggcgaacttgcgccgg                                                                                <  2:1776316‑M1/150‑39 (MQ=255)
                                                                        gccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGc                                                                              <  2:429823/150‑1 (MQ=255)
                                                                                 ttCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGgcggcg                                                                     <  2:1221146/150‑1 (MQ=255)
                                                                                   ctcGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgct                                                                   <  1:752625/150‑1 (MQ=255)
                                                                                    tcGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctg                                                                  <  1:540786/150‑1 (MQ=255)
                                                                                       gTGGGGCTGCTGGCCTCGTCGAGCTACGCGCGAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCGTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAAGCGGGCAAGGCggcggtggcg                                                               >  2:1844906/1‑150 (MQ=255)
                                                                                           ggCTGCTGGCCTCGTCGAGCTCCGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCTGCGCAAGTTCGCCCAGAGCCAGCGCGAGGCGGGCAAGGCGGcgctggcgctgg                                                           <  2:675827/150‑1 (MQ=255)
                                                                                           ggCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGTCCGCCTGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGcgctggcgctgg                                                           <  2:583697/150‑1 (MQ=255)
                                                                                              tgctgGCCACGTCGAGCTACGCGCCAATGGCGGCGCCAATGGCCGCCGGGATGGTGCCGCCGATTGGCTTGTGCATTGCCACCTTCCTGGTCGGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGc                                                        <  2:553781/150‑1 (MQ=255)
                                                                                               gctcGCCTCGCCGACCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGCGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCt                                                       <  1:1844906/150‑1 (MQ=255)
                                                                                                 tgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                                     <  2:1763337/150‑1 (MQ=255)
                                                                                                 tgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                                                     >  2:1826393/1‑150 (MQ=255)
                                                                                                   gCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGc                                                   <  2:1790690/150‑1 (MQ=255)
                                                                                                    ccTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCaga                                                                                               >  1:208401/1‑105 (MQ=255)
                                                                                                    ccTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCaga                                                                                               <  2:208401/105‑1 (MQ=255)
                                                                                                    ccTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCt                                                  >  1:1075361/1‑150 (MQ=255)
                                                                                                     cTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCtt                                                 >  2:656870/1‑150 (MQ=255)
                                                                                                      tcgtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTc                                                <  1:1310228/150‑1 (MQ=255)
                                                                                                        gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTtgcccgcctcgcgctcgctctgggcgaacttgcgccgggccaggaaggtg                                                                    <  1:261627‑M1/128‑51 (MQ=255)
                                                                                                        gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTtgcccgcctcgcgctcgctctgggcgaacttgcgccgggccaggaaggtg                                                                    >  2:261627‑M1/1‑78 (MQ=255)
                                                                                                               tACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAg                                       >  1:1103370/1‑150 (MQ=255)
                                                                                                                 caagcCAATCGGCGGCACCATGCCGGCGGCCATGGTCGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAgg                                      >  2:367180/4‑150 (MQ=37)
                                                                                                                 caagcCAATCGGCGGCACCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAggg                                     >  2:1126609/4‑150 (MQ=25)
                                                                                                                  gcgcCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAgggg                                    >  2:1748978/1‑150 (MQ=255)
                                                                                                                  gcgcCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAgggg                                    >  1:270032/1‑150 (MQ=255)
                                                                                                                  aagcCAATCGGCGGCACCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAgggg                                    <  1:1126609/148‑1 (MQ=25)
                                                                                                                   cgcCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCGCCGAggggg                                   <  1:1119723/150‑1 (MQ=255)
                                                                                                                    gcCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGc                                  >  2:1915645/1‑150 (MQ=255)
                                                                                                                      cAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGa                                >  1:1376207/1‑150 (MQ=255)
                                                                                                                      cAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGa                                >  1:1376216/1‑150 (MQ=255)
                                                                                                                              ggcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGACAGATCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGtt                         >  1:121057/1‑150 (MQ=255)
                                                                                                                               gcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTTCCACCTTCCTGGCCCGGCGCAAGTTCGCCCGGAGGGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttg                        <  1:376079/150‑2 (MQ=255)
                                                                                                                               gcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt                        <  2:755573/150‑1 (MQ=255)
                                                                                                                               gcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt                        <  1:1748978/150‑1 (MQ=255)
                                                                                                                                 gACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgc                      >  2:1816618/1‑150 (MQ=255)
                                                                                                                                    cATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGgcg                                                                                                             <  1:694965/60‑1 (MQ=255)
                                                                                                                                    cATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGgcg                                                                                                             >  2:694965/1‑60 (MQ=255)
                                                                                                                                     aTGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTgccgcc                  >  1:1281355/1‑150 (MQ=255)
                                                                                                                                         ccgccgGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCAGCTCCGAGGGGGCGATTCCGTTTGCCGCCAAgg              <  2:874575/150‑1 (MQ=255)
                                                                                                                                             cgGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGAccc          <  2:1213802/150‑1 (MQ=255)
                                                                                                                                             cgGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGAccc          >  1:1871912/1‑150 (MQ=255)
                                                                                                                                             cgGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGAccc          <  1:412079/150‑1 (MQ=255)
                                                                                                                                                   tGGTGCGGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGACCCGTTAc     <  1:903843/150‑1 (MQ=255)
                                                                                                                                                      tGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGACCCGTTACggg  <  1:1366384/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
AATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTTGCCGCCAAGGACCCGTTACGGG  >  NC_002947/914499‑914796

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