Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,363,111 G→T 5.6% L235I (CTC→ATC)  glmM ← phosphoglucosamine mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,363,1110GT 5.6% 223.0 / 2.3 71L235I (CTC→ATC) glmMphosphoglucosamine mutase
Reads supporting (aligned to +/- strand):  ref base G (35/32);  new base T (2/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

ATCAGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTT  >  NC_002947/5362963‑5363255
                                                                                                                                                    |                                                                                                                                                
atcaGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGAGCGCCGTCACCGTCAAAGGCAATGCCgaga                                                                                                                                                 >  2:1364891/1‑150 (MQ=255)
    gCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCgg                                                                                                                                             <  1:627345/150‑1 (MQ=255)
    gCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCgg                                                                                                                                             >  1:1194702/1‑150 (MQ=255)
     cGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGTGAATTCCGAGATCGGc                                                                                                                                            <  2:808515/150‑1 (MQ=255)
     cGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGc                                                                                                                                            <  1:898917/150‑1 (MQ=255)
       tACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCAt                                                                                                                                          <  2:1416897/150‑1 (MQ=255)
         ccGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGATCTCGGCAttg                                                                                                                                        <  2:1243604/150‑3 (MQ=255)
         ccGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGt                                                                                                                                        >  2:651399/1‑150 (MQ=255)
         ccGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGt                                                                                                                                        >  2:1513307/1‑150 (MQ=255)
           gACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGATCTCGGCAttgcc                                                                                                                                      <  1:373933/150‑5 (MQ=255)
           gACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTcc                                                                                                                                      <  1:1513307/150‑1 (MQ=255)
             ccaccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGa                                                                                                                                    >  1:559889/1‑150 (MQ=255)
             ccaccCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGa                                                                                                                                    <  1:1364891/150‑1 (MQ=255)
                ccccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGAGGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCa                                                                                                                                 >  1:1240158/1‑150 (MQ=255)
                 cccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCAt                                                                                                                                          <  1:419091/140‑1 (MQ=255)
                 cccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCAt                                                                                                                                          >  2:419091/1‑140 (MQ=255)
                 cccGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAg                                                                                                                                >  1:1508801/1‑150 (MQ=255)
                       ttGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGc                                                                                                                          >  1:626664/1‑150 (MQ=255)
                          cAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGc                                                                                                                       >  2:644137/1‑150 (MQ=255)
                              ttGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGc                                                                                                                   >  2:454994/1‑150 (MQ=255)
                                   aTGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGa                                                                                                              >  2:659855/1‑150 (MQ=255)
                                     gTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGAtt                                                                                                            >  2:1444860/1‑150 (MQ=255)
                                        ccTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGGCCGGTTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGa                                                                                                         <  2:1206676/150‑1 (MQ=255)
                                         cTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAt                                                                                                        <  1:77758/150‑1 (MQ=255)
                                          tGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtg                                                                                                       >  2:1082571/1‑150 (MQ=255)
                                           gCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                      <  1:436070/150‑1 (MQ=255)
                                           gCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                      >  1:945988/1‑150 (MQ=255)
                                           gCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgt                                                                                                      <  2:559889/150‑1 (MQ=255)
                                             aGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAtgtgt                                                                                                    <  1:659855/150‑1 (MQ=255)
                                                          tgatgaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTTCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAAc                                                                                       >  1:669726/1‑150 (MQ=255)
                                                          tgatgaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAAc                                                                                       >  2:124956/1‑150 (MQ=255)
                                                          tgatgaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAAc                                                                                       >  2:1241575/1‑150 (MQ=255)
                                                             tgaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCtt                                                                                    >  1:1211245/1‑150 (MQ=255)
                                                              gaACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTc                                                                                   >  2:701339/1‑150 (MQ=255)
                                                                aCAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGt                                                                                 <  1:1444860/150‑1 (MQ=255)
                                                                aCAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGAGGTGTGTCGAGCCGCAACCTTCGt                                                                                 <  2:1508801/150‑1 (MQ=255)
                                                                aCAGCAGCTCGTCACCGTCGACGATGGCGACCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGt                                                                                 <  2:1194702/150‑1 (MQ=255)
                                                                 cagcagCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGtt                                                                                <  1:644137/150‑1 (MQ=255)
                                                                 cagcagCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGtt                                                                                >  2:32832/1‑150 (MQ=255)
                                                                          gTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAg                                                                       >  2:1496250/1‑150 (MQ=255)
                                                                          gTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAg                                                                       <  1:124956/150‑1 (MQ=255)
                                                                          gTCACCGTCGACGATCGCGCTCGTGTGGTCGACCCTCAGAACGCGATCGCCGTCACCGTCAAAGGCAATGCCTAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAg                                                                       <  2:85046/150‑1 (MQ=255)
                                                                           tCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAgg                                                                      <  2:945988/150‑1 (MQ=255)
                                                                            cACCGTCGACGATGGCTCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACATCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGc                                                                     <  1:633434/150‑1 (MQ=255)
                                                                                gTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTc                                                                 >  2:623040/1‑150 (MQ=255)
                                                                                 tCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATNGCCGTCACCGTCAAAGGCAATGCCGGGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTAGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTAAGGCCGTcc                                                                >  1:1068789/1‑150 (MQ=255)
                                                                                 tCGACGATGGCGCACGCCTGGTCGAGCATCGGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTcc                                                                >  2:145554/1‑150 (MQ=255)
                                                                                     cgaCGGCGGCCGGGTGGTCGACCATCAGCACGCGATCGCCCTCAACGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGGTGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGtt                                                            <  2:669726/150‑1 (MQ=255)
                                                                                      gaTGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGcc                                                                    <  2:1130962/141‑1 (MQ=255)
                                                                                      gaTGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGcc                                                                    >  1:1130962/1‑141 (MQ=255)
                                                                                        tGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTgcg                                                         >  2:156136/1‑150 (MQ=255)
                                                                                            gcCCGTGGGGTCGACCAGCAGCCCGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATg                                                     <  2:71589/150‑1 (MQ=255)
                                                                                                   tgGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCTGTGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTTTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACg                                              <  2:1135060/150‑1 (MQ=255)
                                                                                                   tgGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCgg                                                              <  2:1127499/134‑1 (MQ=255)
                                                                                                   tgGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCgg                                                              >  1:1127499/1‑134 (MQ=255)
                                                                                                     gTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGt                                            <  1:623040/150‑1 (MQ=255)
                                                                                                      tCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTc                                           >  1:10056/1‑150 (MQ=255)
                                                                                                      tCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTc                                           >  1:10299/1‑150 (MQ=255)
                                                                                                          ccATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGt                                       >  2:151495/1‑150 (MQ=255)
                                                                                                                 cACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCAcc                                >  1:1428281/1‑150 (MQ=255)
                                                                                                                   cgcgATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCa                              <  2:626664/150‑1 (MQ=255)
                                                                                                                         cGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCAc                        <  1:1084204/150‑1 (MQ=255)
                                                                                                                                  cgtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGc               >  2:1578964/1‑150 (MQ=255)
                                                                                                                                   gtcaAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCt              <  1:1496250/150‑1 (MQ=255)
                                                                                                                                      aaaGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTggg           >  2:486054/1‑150 (MQ=255)
                                                                                                                                      aaaGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGAACGTCCGGTTGCGAAAGACACACAGTCAAGTCGGCAACCAGCTCACGGAATACGCTggg           >  1:876009/1‑150 (MQ=255)
                                                                                                                                      aaaGGAAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTggg           <  1:651399/150‑1 (MQ=255)
                                                                                                                                       aaGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGgc          <  1:486054/150‑1 (MQ=255)
                                                                                                                                         ggacaTGCCGATCTCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGgcgc        >  1:1009651/5‑150 (MQ=255)
                                                                                                                                             aTGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAAcc    <  1:1241575/150‑1 (MQ=255)
                                                                                                                                              tGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCt   >  2:1092894/1‑150 (MQ=255)
                                                                                                                                               gCCGATCTCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGGTAAAGATTTCCGGTTGCGCATGCAGCACGGTCACGTCGGCCCCGCGCTTACGGAAGACGCCGGGCGCAACCtt  >  1:1593985/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
ATCAGCGTACCGACCACCCCGCCTTGCAGCTTGCCATGTTCCTGCAGGTCGCGAGCAATGATGAACAGCAGCTCGTCACCGTCGACGATGGCGCCCGTGTGGTCGACCATCAGCACGCGATCGCCGTCACCGTCAAAGGCAATGCCGAGATCGGCATGTCCGACCAGTACAGCGGCCTGCAGCGATTCGATGTGTGTCGAGCCGCAACCTTCGTTGATGTTCAGGCCGTCCGGTTGCGCATGCAGCACGGTCACGTCGGCACCCAGCTCACGGAATACGCTGGGCGCAACCTT  >  NC_002947/5362963‑5363255

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: