Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,776,107 A→C 5.7% K152Q (AAG→CAG)  sotB → sugar efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,776,1070AC 5.7% 223.7 / 2.4 77K152Q (AAG→CAG) sotBsugar efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (34/38);  new base C (3/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAT  >  NC_002947/2775963‑2776252
                                                                                                                                                   |                                                                                                                                                   
aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCCAgcg                                                                                                                                                 >  3:250569/1‑147 (MQ=255)
aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGgc                                                                                                                                                 <  7:275374/149‑1 (MQ=255)
aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGgc                                                                                                                                                 >  6:299209/1‑149 (MQ=255)
      aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTgg                                                                                                                                           <  6:127128/148‑1 (MQ=255)
        ccACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCC‑‑‑GCCCGCCAAGGCGCTggggc                                                                                                                                        >  2:292535/1‑147 (MQ=255)
          aCTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCctg                                                                                                                                      <  1:192160/149‑1 (MQ=255)
           cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCctgc                                                                                                                                     <  6:249164/149‑1 (MQ=255)
              gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCCAGCGCCTTGGCctgctg                                                                                                                                   <  4:250569/148‑1 (MQ=255)
              gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCctgctg                                                                                                                                   <  1:382778/148‑1 (MQ=255)
                  tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCCGGCCAAGGCGCTGGGCCTGCAGGcccc                                                                                                                              >  4:248930/1‑147 (MQ=255)
                      ggCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                          >  8:255741/1‑148 (MQ=255)
                       gCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGt                                                                                                                         <  3:55618/148‑1 (MQ=255)
                             gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTAccc                                                                                                                  >  8:234179/1‑149 (MQ=255)
                             gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTAccc                                                                                                                  >  3:282877/1‑149 (MQ=255)
                             gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTAcc                                                                                                                   >  4:42319/1‑148 (MQ=255)
                                 aGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑CCAGGCCAAGGCGGGGGGCCTGCGGGCCCCCCCTACTACCCAGGc                                                                                                              >  1:144689/1‑149 (MQ=255)
                                  gCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGcc                                                                                                             <  8:335056/149‑1 (MQ=255)
                                   cTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCa                                                                                                            <  2:222189/149‑1 (MQ=255)
                                    ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCctgc                                                                                                                                     <  2:306864/124‑1 (MQ=255)
                                    ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCctgc                                                                                                                                     >  1:306864/1‑124 (MQ=255)
                                      cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATgg                                                                                                         >  2:276093/1‑149 (MQ=255)
                                      cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATgg                                                                                                         >  3:65662/1‑149 (MQ=255)
                                         cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCaa                                                                                                                                                    <  8:3598/105‑1 (MQ=255)
                                         cATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCaa                                                                                                                                                    >  7:3598/1‑105 (MQ=255)
                                          aTGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑ACAGGCCAACTCGCTGGGCCAGCAGGCCACCCGTACTACCCTAGACAAGGTc                                                                                                       >  8:351978/1‑147 (MQ=255)
                                           tgctgcTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTg                                                                                                    <  8:231238/149‑1 (MQ=255)
                                            gctgctGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGGCTGCTGGCCACCGGGACTAACCTGGCCAGGGTCCTgg                                                                                                   >  6:205194/1‑149 (MQ=255)
                                             ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGTCCATGGTCCTgg                                                                                                   <  8:326896/148‑1 (MQ=255)
                                             ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                   <  2:350869/148‑1 (MQ=255)
                                             ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                   <  3:302063/148‑1 (MQ=255)
                                              tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTggg                                                                                                  >  6:355109/1‑148 (MQ=255)
                                                   gTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtactac                                                                                                                    >  6:322157/1‑125 (MQ=255)
                                                   gTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtactac                                                                                                                    <  5:322157/125‑1 (MQ=255)
                                                      aGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCt                                                                                         >  4:135578/1‑149 (MQ=255)
                                                       gCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTg                                                                                        <  7:47857/149‑1 (MQ=255)
                                                          gCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCTTCCCGCTGGGc                                                                                     >  2:34617/1‑149 (MQ=255)
                                                                 cATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAAAGGCAGGCCAAAGCGCGAGGGCCTGCTGGCCACCGGTACTAACCTGGCCATGGTCCTGGGCATCCCGCTGGGCCgt                                                                                  >  1:122126/1‑149 (MQ=255)
                                                                  aTTGCCCTGGCGCATGCGGTGTTCATGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑CCATGCCAAGGCGCAGGGGCTGCTGGCCACCCGTACTAACCAGGCCAAGGAACATGGCAAACCGCAGGGACgtgtt                                                                               >  6:211495/1‑146 (MQ=255)
                                                                     gCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCTTCCCGCTGGGCCGTGtggtgg                                                                           <  1:34617/148‑1 (MQ=255)
                                                                      cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACGACCCAGGCCATGGTCCGGGGCCTCCCGCTGGGCCGGGTGGTGGGc                                                                         >  8:349760/1‑149 (MQ=255)
                                                                      cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGCCCTGCTGGCCCCCGCGACAACCCTCGCCCTGGTCCTGGGCATACCGCAGGGCCCTCtggtggcc                                                                         >  3:288052/1‑147 (MQ=255)
                                                                        cTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                       <  7:351978/149‑1 (MQ=255)
                                                                        cTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                       <  3:135578/149‑1 (MQ=255)
                                                                           gcgcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGta                                                                                                                        <  3:70002/97‑1 (MQ=255)
                                                                           gcgcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGta                                                                                                                        >  4:70002/1‑97 (MQ=255)
                                                                           gcgcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAg                                                                      <  3:42319/147‑1 (MQ=255)
                                                                               aTGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTGCTGGGCATACCGCTGGGCCGTGTGGTGGGCGCGGCACTg                                                                >  3:182473/1‑149 (MQ=255)
                                                                                   ggTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCt                                                                                         <  5:235522/120‑1 (MQ=255)
                                                                                            gTCCATCTCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTa                                                      <  1:133916/146‑1 (MQ=255)
                                                                                            gTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCGGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCa                                                   >  4:357946/1‑149 (MQ=255)
                                                                                               cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCCCGt                                                >  5:277594/1‑149 (MQ=255)
                                                                                               cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGt                                                <  4:65662/149‑1 (MQ=255)
                                                                                                aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTggg                                                              <  3:93371/134‑1 (MQ=255)
                                                                                                aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGtt                                               >  5:361612/1‑149 (MQ=255)
                                                                                                 tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTc                                              <  4:282877/149‑1 (MQ=255)
                                                                                                   aCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCt                                            <  5:355109/149‑1 (MQ=255)
                                                                                                     cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                          <  2:174271/149‑1 (MQ=255)
                                                                                                     cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                          <  7:349760/149‑1 (MQ=255)
                                                                                                       ccTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                         <  5:120104/132‑1 (MQ=255)
                                                                                                           gCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCAtt                                     <  4:244952/148‑1 (MQ=255)
                                                                                                               gCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                <  2:122126/149‑1 (MQ=255)
                                                                                                                       tgtgGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcgc                        <  3:248930/149‑1 (MQ=255)
                                                                                                                           gcgcCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAg                     >  1:86055/1‑149 (MQ=255)
                                                                                                                           gcgcCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAg                     <  2:86055/148‑1 (MQ=255)
                                                                                                                             gccgccGGGCCA‑‑‑GCAGGCCAACCCGCGGGGCCGGCTGGCCCCCCGTACGACCCTGGCCCAGGGCCGGGGGCTCCCCCGGGGGCGTGTGGTGGGCGCGGCACGGGGCGGGCGTTTCACGTTGCGGCGCATTGCCGGGGTAGGGGTAGccc                  >  6:6162/1‑148 (MQ=255)
                                                                                                                               cgccgGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacc                <  7:255741/149‑1 (MQ=255)
                                                                                                                                    ggCCA‑‑‑GCAGGCCCAGCGCCTTGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGGGGCACTGGGCTGGCGTATCACGTTCCt                                            <  2:175897/116‑1 (MQ=38)
                                                                                                                                    ggCCA‑‑‑GCAGGCCCAGCGCCTTGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGGGGCACTGGGCTGGCGTATCACGTTCCt                                            >  1:175897/1‑116 (MQ=38)
                                                                                                                                       ca‑‑‑gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGGATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTg        >  6:259302/1‑149 (MQ=255)
                                                                                                                                        a‑‑‑gcagGCCCAGCGCCTTGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc       >  2:26696/1‑149 (MQ=255)
                                                                                                                                        a‑‑‑gcagGCCAAGGCGCTGGGCCGGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATGACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc       <  1:276093/149‑1 (MQ=255)
                                                                                                                                            gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCCTCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc      >  7:143893/1‑149 (MQ=255)
                                                                                                                                            gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc      <  5:205194/149‑1 (MQ=255)
                                                                                                                                                gCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                         >  3:151096/1‑94 (MQ=255)
                                                                                                                                                gCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                         <  4:151096/94‑1 (MQ=255)
                                                                                                                                                gCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTCGCGCTAGCCACCATGCTTGCCTGAt  >  8:155987/1‑148 (MQ=255)
                                                                                                                                                gCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt  >  7:356142/1‑149 (MQ=255)
                                                                                                                                                gCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt  <  7:234179/149‑1 (MQ=255)
                                                                                                                                                 ccAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt  >  4:340539/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAT  >  NC_002947/2775963‑2776252

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