Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 912,283 C→A 5.0% L930M (CTG→ATG)  fruB → PTS fructose transporter subunit EI/HPr/IIA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947912,2830CA 5.0% 314.8 / 3.0 98L930M (CTG→ATG) fruBPTS fructose transporter subunit EI/HPr/IIA
Reads supporting (aligned to +/- strand):  ref base C (44/49);  new base A (2/3);  total (46/52)
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.48e-01

GGGCGTGTGTGGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGCCCAGGGCAG  >  NC_002947/912138‑912426
                                                                                                                                                 |                                                                                                                                               
gggCGTGTGTGGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGgc                                                                                                                                             <  1:1253196/150‑1 (MQ=255)
gggCGTGTGTGGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGgc                                                                                                                                             <  1:1254371/150‑1 (MQ=255)
  gCGTGTGTGGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGgcgc                                                                                                                                           <  1:676476/150‑1 (MQ=255)
    gtgtgtGGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGgcgcgc                                                                                                                                         >  1:652602/1‑150 (MQ=255)
     tgtgtgGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCACAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCAGGgcgcggg                                                                                                                                        <  2:2293947/150‑3 (MQ=255)
        gtgGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAgg                                                                                                                                     <  1:814073/150‑1 (MQ=255)
          ggCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCAGGgcgcgggctgtc                                                                                                                                   <  2:2223880/150‑8 (MQ=255)
            cGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCTGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCt                                                                                                                                 <  2:1600501/150‑1 (MQ=255)
             gAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTg                                                                                                                                >  1:570930/1‑150 (MQ=255)
                  ggCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACa                                                                                                                           <  1:2139058/150‑1 (MQ=255)
                   gCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTAAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAg                                                                                                                          <  2:1851930/150‑1 (MQ=255)
                     aGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGAACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGa                                                                                                                        >  1:1476753/1‑150 (MQ=255)
                         gACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGcc                                                                                                                                   <  2:1771742/135‑1 (MQ=255)
                         gACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGcc                                                                                                                                   >  1:1771742/1‑135 (MQ=255)
                           cccGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGCCGACCTCAGCCTCTCCGCGCCCTTCATGACCGGAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAgcgc                                                                                                                  <  1:1860751/150‑1 (MQ=255)
                           cccGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAgcgc                                                                                                                  <  1:1876141/150‑1 (MQ=255)
                             cGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACa                                                                                                                           <  1:1255117/139‑1 (MQ=255)
                             cGCAGGCGGTGGCCGTGCTGCTGGGGCTAGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCGTCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACa                                                                                                                           >  2:1255117/1‑139 (MQ=255)
                              gCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGcgccgc                                                                                                               >  2:2389183/1‑150 (MQ=255)
                              gCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGcgccgc                                                                                                               >  2:2383433/1‑150 (MQ=255)
                               cAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCgca                                                                                                              <  1:1472767/150‑1 (MQ=255)
                                 ggCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGAGCCgcagc                                                                                                            >  2:298490/1‑150 (MQ=255)
                                  gCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCgcagca                                                                                                           <  2:555294/150‑1 (MQ=255)
                                   cGGGGGCCGTGCTGCTGGGGCTGGACGTGTACTAACTCAGCGTCTCCGCGCGCAACATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGcagcag                                                                                                          <  1:834799/150‑1 (MQ=255)
                                     gTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGtt                                                                                                        <  2:1008258/150‑1 (MQ=255)
                                      tGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTc                                                                                                       <  2:2059064/150‑1 (MQ=255)
                                          cGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTcgcgc                                                                                                   <  1:1423796/150‑1 (MQ=255)
                                              ctgctgGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCtg                                                                                               <  1:2062803/150‑1 (MQ=255)
                                                gctgGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCtggt                                                                                             <  2:2120390/150‑1 (MQ=255)
                                                 ctgGGGCTGGACGTGGACGCACTCAGCGTCTCCCCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCtggtg                                                                                            <  1:839639/150‑1 (MQ=255)
                                                 ctgGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGccc                                                                                                                                  >  2:387665/1‑112 (MQ=255)
                                                 ctgGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGccc                                                                                                                                  >  2:387600/1‑112 (MQ=255)
                                                 ctgGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGccc                                                                                                                                  <  1:387600/112‑1 (MQ=255)
                                                 ctgGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGccc                                                                                                                                  <  1:387665/112‑1 (MQ=255)
                                                     ggCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGc                                                                                        <  1:1662771/150‑1 (MQ=255)
                                                     ggCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGc                                                                                        <  1:1423201/150‑1 (MQ=255)
                                                      gCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCg                                                                                       >  1:143151/1‑150 (MQ=255)
                                                      gCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCAGGgcgcgg                                                                                                                                         <  1:970023/100‑2 (MQ=255)
                                                      gCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCAGGgcgcgg                                                                                                                                         >  2:970023/1‑99 (MQ=255)
                                                      gCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCAATCGCTCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCg                                                                                       >  1:1207826/1‑150 (MQ=255)
                                                        tggacgtggacgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCACCAGACAGCCCGCGCACTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTAGCGCGCTGGTGGAGCGtt                                                                                     >  1:730752/1‑150 (MQ=255)
                                                          gacgtggacgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTAc                                                                                   >  2:234193/1‑150 (MQ=255)
                                                           acgtggacgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACt                                                                                  >  2:596265/1‑150 (MQ=255)
                                                             gtggacgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGa                                                                                >  1:705239/1‑150 (MQ=255)
                                                               ggacgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAAt                                                                              <  1:357900/150‑1 (MQ=255)
                                                                gacgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATg                                                                             <  2:1555141/150‑1 (MQ=255)
                                                                gacgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATg                                                                             <  2:1746890/150‑1 (MQ=255)
                                                                  cgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGc                                                                           >  2:811255/1‑150 (MQ=255)
                                                                  cgAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGc                                                                           >  1:1011880/1‑150 (MQ=255)
                                                                   gAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGcc                                                                          >  1:80285/1‑150 (MQ=255)
                                                                   gAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGcc                                                                          >  2:2318913/1‑150 (MQ=255)
                                                                     aCTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGa                                                                                                                        >  1:793863/1‑102 (MQ=255)
                                                                     aCTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGa                                                                                                                        <  2:793863/102‑1 (MQ=255)
                                                                      cTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAg                                                                       >  1:1961675/1‑150 (MQ=255)
                                                                       tCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCgcagca                                                                                                           >  1:423634/1‑113 (MQ=255)
                                                                       tCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCgcagca                                                                                                           <  2:423634/113‑1 (MQ=255)
                                                                       tCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGa                                                                      >  1:2223225/1‑150 (MQ=255)
                                                                          gCGTCTCCGCGCGCAGCATTGCCGAAGTCCAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTTGTGCGAGAGGCCCTGCAACAGGACAGCGCCGACGCCGTTCCCGCGCTGGTGGAGCGTTTCTTGATGGCCAAGATcc                                                                   <  2:521629/150‑1 (MQ=255)
                                                                          ccGTCTCCGCGCGCAGCATCGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATcc                                                                   <  2:828341/149‑1 (MQ=255)
                                                                            gTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGc                                                                                                                                    >  1:1069844/1‑83 (MQ=255)
                                                                            gTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGc                                                                                                                                    <  2:1069844/83‑1 (MQ=255)
                                                                                       cagcaTTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCt                                                      <  2:80285/150‑1 (MQ=255)
                                                                                          caTTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGac                                                   >  2:979073/1‑149 (MQ=255)
                                                                                          caTTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGaa                                                   <  2:652602/150‑1 (MQ=255)
                                                                                           aTTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAAc                                                  >  1:1191436/1‑150 (MQ=255)
                                                                                              gCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCt                                               <  1:203066/150‑1 (MQ=255)
                                                                                               ccGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTg                                              >  1:861914/1‑150 (MQ=255)
                                                                                               ccGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGAGCTGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTg                                              >  1:149052/1‑150 (MQ=255)
                                                                                                 gAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTgcg                                            >  1:512587/1‑150 (MQ=255)
                                                                                                   aGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCt                                          <  2:1961675/150‑1 (MQ=255)
                                                                                                       aaGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTc                                                                                                       >  1:361499/1‑85 (MQ=255)
                                                                                                       aaGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTc                                                                                                       <  2:361499/85‑1 (MQ=255)
                                                                                                       aaGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGat                                      <  1:811255/150‑1 (MQ=255)
                                                                                                       aaGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGat                                      >  1:33307/1‑150 (MQ=255)
                                                                                                          gCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATc                                   <  2:570930/150‑1 (MQ=255)
                                                                                                           cACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCAc                                                               <  1:621420/121‑1 (MQ=255)
                                                                                                           cACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCAc                                                               >  2:621420/1‑121 (MQ=255)
                                                                                                            aCTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATgg                                                                            <  2:175750/107‑1 (MQ=255)
                                                                                                            aCTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATgg                                                                            >  1:175750/1‑107 (MQ=255)
                                                                                                                 ttCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCg                            <  1:963181/150‑1 (MQ=255)
                                                                                                                        ggCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGa                     <  1:2383433/150‑1 (MQ=255)
                                                                                                                        ggCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGa                     <  2:861914/150‑1 (MQ=255)
                                                                                                                              tcatcaGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACActct               <  1:596265/150‑1 (MQ=255)
                                                                                                                                atcaGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTgc             >  1:1806508/1‑150 (MQ=255)
                                                                                                                                atcaGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTga             >  1:232199/1‑149 (MQ=255)
                                                                                                                                atcaGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGAGCTGGTGGAGAGTTACTGCATGACAAAGATTCTACTTCACAGCGTGAGACCAGCGATGGAAATCTCCATCGTCCTGGACACTCTgc             >  2:745150/1‑150 (MQ=255)
                                                                                                                                 tcaGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCa                                  <  2:443734/128‑1 (MQ=255)
                                                                                                                                 tcaGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCa                                  >  1:443734/1‑128 (MQ=255)
                                                                                                                                   aGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGcc          >  1:2353747/1‑150 (MQ=255)
                                                                                                                                    gACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGccc         >  2:862476/1‑150 (MQ=255)
                                                                                                                                    gACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCCCTGGATATCACCATCGGCCTGGACACTCTGCGccc         >  2:2229444/1‑150 (MQ=255)
                                                                                                                                          ccGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGTATAGCACCATCTGCCTGGACACTCTGCGCCCAGGGCa   >  2:1592591/1‑150 (MQ=255)
                                                                                                                                          ccGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGCCCAGGGCa   >  2:670494/1‑150 (MQ=255)
                                                                                                                                          ccGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGCCCAGGGCa   <  2:508187/150‑1 (MQ=255)
                                                                                                                                           cgcgCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGCCCAGGGCAg  >  1:1547962/1‑150 (MQ=255)
                                                                                                                                           cgcgCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGCCCAGGGCAg  <  1:2417840/150‑1 (MQ=255)
                                                                                                                                            gcgcCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATcc                                                                   >  1:374614/1‑84 (MQ=255)
                                                                                                                                            gcgcCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATcc                                                                   <  2:374614/84‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GGGCGTGTGTGGCGAGCTGGCAGCCGACCCGCAGGCGGTGGCCGTGCTGCTGGGGCTGGACGTGGACGAACTCAGCGTCTCCGCGCGCAGCATTGCCGAAGTCAAGGCACTGGTTCGCCAGGCCGATCATCAGACAGCCCGCGCCCTGGCGCGCGAGGCCCTGCAACAGGACAGCGCCGCAGCAGTTCGCGCGCTGGTGGAGCGTTACTGAATGGCCAAGATCCTCACTCTCACCCTGAACCCTGCGCTGGATATCACCATCGGCCTGGACACTCTGCGCCCAGGGCAG  >  NC_002947/912138‑912426

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: