Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 914,627 Δ1 bp 9.1% coding (1335/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,6270G. 9.1% 189.3 / 15.2 66coding (1335/1743 nt)fruAPTS fructose transporter subunit IIBC
Reads supporting (aligned to +/- strand):  ref base G (32/28);  new base . (4/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

GACAGCATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTT  >  NC_002947/914481‑914774
                                                                                                                                                   |                                                                                                                                                     
gACAGTATGTGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGAcc                                                                                                                                                    >  2:565687/1‑150 (MQ=255)
gACAGCATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGAcc                                                                                                                                                    >  1:496441/1‑150 (MQ=255)
   tgCATGGGCACCACCAAGGCCATCCTGCGGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATg                                                                                                                                                 <  2:1119644/149‑1 (MQ=255)
      aTGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGgcc                                                                                                                                              >  2:144668/1‑150 (MQ=255)
       tGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGgccg                                                                                                                                             >  1:583562/1‑150 (MQ=255)
        gggCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGgccgc                                                                                                                                            >  1:6657/1‑150 (MQ=255)
        gggCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGgccgc                                                                                                                                            >  1:952102/1‑150 (MQ=255)
        gggCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGgccgc                                                                                                                                            >  1:180103/1‑150 (MQ=255)
        gggCACAACCAATGCCATCCTGCTGGGGCTGTTGATGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGgccgc                                                                                                                                            >  1:214361/1‑150 (MQ=255)
          gCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGccgccg                                                                                                                                          >  2:1121588/1‑150 (MQ=255)
             gggccaGTGCCGTCCTGCTGGGGCTGTTGCTGGGGGGCGGGAGGTGCGTCGCCTGTGGCGGGCAGTTCAGCACGGCGGCCTCTGCGCTCTCGGTGGTGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGGCGGCa                                                                                                                                      <  2:807185/147‑1 (MQ=255)
               accaATGCCATCCCGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGCTGACGCTGCTGGCCTCGTCGAGCTGCGCGCGAATGGCGGCGACCATGGCCGCCGGCATg                                                                                                                                    >  1:561692/1‑150 (MQ=255)
                    tGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgcc                                                                                                                              >  1:45910/1‑150 (MQ=255)
                         tCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGa                                                                                                                         >  1:381578/1‑150 (MQ=255)
                         tCCTGCTGGGGCTGGTGCTGGGGGGCATGATGTGCGTCGACGTTGGGGGGCCGATCAACAAGGCCGCCTGTGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGa                                                                                                                         <  2:414593/150‑1 (MQ=255)
                           ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGTCTGCTATTCGTTCTAACAGCAGCTGCTGGCCTCGCCGAGCTACCCCCCCCTGGCGTCGACCATGGCCGCCGGCATGGTGCCGCCGAtt                                                                                                                       <  1:630131/150‑1 (MQ=255)
                           ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATCGGCGGCACCATGCCGGCGGCCATGGTCgccgccatt                                                                                                                        <  1:615021/150‑4 (MQ=37)
                           ctgctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGAtt                                                                                                                       <  2:381578/150‑1 (MQ=255)
                             gctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTg                                                                                                                                >  2:336757/1‑139 (MQ=255)
                             gctgGGGCTGTTGCTGGGCGGCATGATATTCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTg                                                                                                                                <  1:336757/139‑1 (MQ=255)
                              ctgGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGc                                                                                                                    <  1:277097/150‑1 (MQ=255)
                                 gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTg                                                                                                                 <  2:206079/150‑1 (MQ=255)
                                 gggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTg                                                                                                                 >  1:28530/1‑150 (MQ=255)
                                  ggCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGAATGGTGCCGCCGATTGGCTTgg                                                                                                                <  1:645674/150‑1 (MQ=255)
                                    tggTTGCTGGGCGGCATGATGTGCGTGGAGCGTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGt                                                                                                              <  2:521140/148‑1 (MQ=255)
                                            gggCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCAc                                                                                                      <  2:583562/150‑1 (MQ=255)
                                                     gatgTGCGTCGACCTTGGTGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgcc                                                                                                                              <  1:1074030/117‑1 (MQ=255)
                                                     gatgTGCGTCGACCTTGGTGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTgcc                                                                                                                              >  2:1074030/1‑117 (MQ=255)
                                                          gCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGgc                                                                                        <  2:1030125/150‑1 (MQ=255)
                                                            gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGGTGCTAGCCTCGTCGAGCTACGCTCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTAGGTATTGCCACCTTCCTGGCCCGgcac                                                                                      >  2:537576/1‑148 (MQ=255)
                                                            gTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTGCCTGGCCCGgcgc                                                                                      <  1:265816/150‑1 (MQ=255)
                                                                 ccTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGtt                                                                                 <  1:144668/150‑1 (MQ=255)
                                                                    tGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGc                                                                              <  2:952102/150‑1 (MQ=255)
                                                                        gggCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCATGGCGGGGCGCAAGTTCGCCCag                                                                          >  2:527273/1‑150 (MQ=255)
                                                                               taaacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgagc                                                                   >  2:775616/3‑150 (MQ=255)
                                                                                caacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGagcgagcg                                                                  >  1:958607/1‑150 (MQ=255)
                                                                                 aacaaGGCCGCCTATGCGTTCTCTGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                 >  2:556646/1‑150 (MQ=255)
                                                                                 aacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCAc                                                                                                      >  1:387487/1‑113 (MQ=255)
                                                                                 aacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCAc                                                                                                      <  2:387487/113‑1 (MQ=255)
                                                                                 aacaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAgcgc                                                                 >  1:396398/1‑150 (MQ=255)
                                                                                  acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGccgccgttg                                                                                                                       <  1:887855/95‑4 (MQ=255)
                                                                                  acaaGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGccgccgttg                                                                                                                       >  2:887855/1‑92 (MQ=255)
                                                                                      ggCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGc                                                            <  2:902688/150‑1 (MQ=255)
                                                                                        ccgccTATGCGTTCTCGACTGGGGTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCAGGGTGCAGCCGATTTGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                          <  1:607361/150‑1 (MQ=255)
                                                                                          gccTATGCGTTCTCGGTGGGGGTGGTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCCGCGCAAGTTCGCCCAGAGCGAGCGCGTGGCGGGc                                                        >  2:278783/1‑150 (MQ=255)
                                                                                           ccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGa                                                               <  1:657316/142‑1 (MQ=255)
                                                                                           ccTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGa                                                               >  2:657316/1‑142 (MQ=255)
                                                                                                  gTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAggcggc                                                <  2:180103/150‑1 (MQ=255)
                                                                                                     ctcGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAggcgcctct                                             >  1:1122846/1‑145 (MQ=255)
                                                                                                      tcGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCggcgctg                                            <  1:96994/150‑1 (MQ=255)
                                                                                                                   tgGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCtgt                               >  1:315354/1‑150 (MQ=255)
                                                                                                                     gCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGc                             <  2:396398/150‑1 (MQ=255)
                                                                                                                      ccTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCt                            >  2:737133/1‑150 (MQ=255)
                                                                                                                       cTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCtt                           >  1:1072770/1‑150 (MQ=255)
                                                                                                                          gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                          >  2:574971/1‑116 (MQ=255)
                                                                                                                          gtcgAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCgg                                                          <  1:574971/116‑1 (MQ=255)
                                                                                                                                 tACGCGCCAATGGCGGCGACCATGGCCGCCGGCCTGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAg                 >  2:1079986/1‑150 (MQ=255)
                                                                                                                                 tACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAg                 <  2:6657/150‑1 (MQ=255)
                                                                                                                                 tACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAg                 >  2:377287/1‑150 (MQ=255)
                                                                                                                                  aCGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCCAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAgg                <  2:496441/150‑1 (MQ=255)
                                                                                                                                       ccAATGGCGGCGACCATGGCCGCCGGCATGGTTCCGCCGATTGGCTTGTGCATTGCCACCTTCCTGGCCCGGCGCACGTTCGCCCAGTGCGTGCGCGATGCTGGCAATGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGACGGGGCg           <  2:214361/150‑1 (MQ=255)
                                                                                                                                        caatcggcggcACCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGa          >  2:984146/6‑150 (MQ=37)
                                                                                                                                         aatcggcggcACCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGAt         >  2:350441/5‑150 (MQ=37)
                                                                                                                                         aatcggcggcACCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGAt         >  1:560147/5‑150 (MQ=37)
                                                                                                                                          atcggcggcACCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgag                                                                    <  2:362492/87‑1 (MQ=12)
                                                                                                                                          atcggcggcACCATGCCGGCGGCCATGGTCGCCGCCATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAgagcgag                                                                    >  1:362492/4‑90 (MQ=12)
                                                                                                                                              gcggcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCTGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCg     >  1:180301/1‑150 (MQ=255)
                                                                                                                                                ggcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCaga                                                                         <  1:1130428/80‑1 (MQ=255)
                                                                                                                                                ggcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCaga                                                                         >  2:1130428/1‑80 (MQ=255)
                                                                                                                                                ggcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGtt   <  1:1121588/150‑1 (MQ=255)
                                                                                                                                                ggcgACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGtt   <  2:1072770/150‑1 (MQ=255)
                                                                                                                                                 gcgACCTTGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGttt  <  2:45910/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACAGCATGGGCACCACCAATGCCATCCTGCTGGGGCTGTTGCTGGGCGGCATGATGTGCGTCGACCTTGGCGGGCCGATCAACAAGGCCGCCTATGCGTTCTCGGTGGGGCTGCTGGCCTCGTCGAGCTACGCGCCAATGGCGGCGACCATGGCCGCCGGCATGGTGCCGCCGATTGGCTTGGGTATTGCCACCTTCCTGGCCCGGCGCAAGTTCGCCCAGAGCGAGCGCGAGGCGGGCAAGGCGGCGCTGGCGCTGGGGCTGTGCTTCATCTCCGAGGGGGCGATTCCGTTT  >  NC_002947/914481‑914774

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