Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,776,107 A→C 5.3% 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.3% 251.3 / 2.0 82K152Q (AAG→CAG) sotBsugar efflux transporter
Reads supporting (aligned to +/- strand):  ref base A (34/43);  new base C (2/3);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATC  >  NC_002947/2775962‑2776253
                                                                                                                                                    |                                                                                                                                                   
cATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAg                                                                                                                                                   <  3:341152/148‑1 (MQ=255)
 aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAg                                                                                                                                                   <  5:30180/147‑1 (MQ=255)
 aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAg                                                                                                                                                   >  6:30180/1‑147 (MQ=255)
 aTTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGgc                                                                                                                                                 >  1:230243/1‑149 (MQ=255)
     tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCt                                                                                                                                              <  5:382040/148‑1 (MQ=255)
     tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTg                                                                                                                                             <  3:31270/149‑1 (MQ=255)
     tCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTg                                                                                                                                             <  6:239592/149‑1 (MQ=255)
        gCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                          <  8:220434/149‑1 (MQ=255)
         ccACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGcc                                                                                                                                         >  7:318007/1‑149 (MQ=255)
            cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgc                                                                                                                                      <  8:143944/149‑1 (MQ=255)
            cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgc                                                                                                                                      <  2:239572/149‑1 (MQ=255)
             ttGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCATGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAACGCCCAGGGCcggct                                                                                                                                     >  2:328346/1‑149 (MQ=255)
               gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                          >  5:62128/1‑142 (MQ=255)
               gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGc                                                                                                                                          <  6:62128/142‑1 (MQ=255)
               gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCctgctg                                                                                                                                    <  3:179968/148‑1 (MQ=255)
                   tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCAc                                                                                                                               <  3:345142/149‑1 (MQ=255)
                     cTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCAGGCCCAGGCGCGGGGCCGGCGGGCccccc                                                                                                                             >  1:347582/1‑148 (MQ=255)
                     cTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGGGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCg                                                                                                                             <  6:229136/149‑1 (MQ=255)
                      tGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCACATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCgg                                                                                                                            >  2:67046/1‑149 (MQ=255)
                        gCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGt                                                                                                                           <  8:137861/148‑1 (MQ=255)
                          tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgc                                                                                                                                      >  4:262123/1‑135 (MQ=255)
                          tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGCGCCTTGGCctgc                                                                                                                                      <  3:262123/135‑1 (MQ=255)
                          tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtac                                                                                                                         >  3:78670/1‑148 (MQ=255)
                          tttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCGGCTGGCCAcccggac                                                                                                                         >  4:19179/1‑143 (MQ=255)
                           ttCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAACAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtacta                                                                                                                       >  7:134282/1‑149 (MQ=255)
                             cGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGtactac                                                                                                                      <  7:242322/148‑1 (MQ=255)
                              gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTAccc                                                                                                                    >  1:64292/1‑149 (MQ=255)
                                aaaGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAAGCGCTGGGCCAGCTGGCCACCGGTACTTCCCTg                                                                                                                  >  7:319200/1‑149 (MQ=255)
                                    cTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCa                                                                                                              <  4:52338/149‑1 (MQ=255)
                                     ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCATGGCCATCTGGCCCCCGGTACTACCCAGGCCAt                                                                                                             >  5:324609/1‑149 (MQ=255)
                                      tCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGCCCAACGCGCTGGGCCGGCTGGCCACCGGTACAAACCTGGCCATg                                                                                                            >  2:384766/1‑149 (MQ=255)
                                       cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATg                                                                                                            <  4:239107/148‑1 (MQ=255)
                                         ccATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCCCCGGTACTACCCTGGCCATGGTc                                                                                                         >  7:227553/1‑149 (MQ=255)
                                              ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                     <  2:203532/148‑1 (MQ=255)
                                              ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                     <  4:355604/148‑1 (MQ=255)
                                              ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTgg                                                                                                     <  2:230243/148‑1 (MQ=255)
                                               tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCGCAGGCCAAGGCGCGGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTggg                                                                                                    >  6:7250/1‑148 (MQ=255)
                                                gctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCCCGCCAAGGCGCAGGGGCGGCTGGCCCCCCGTACAACCCTGGCCATGGTACTGGG‑Ca                                                                                                 >  5:212675/1‑149 (MQ=255)
                                                 ctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGGCATGGTCCTGGG‑CAt                                                                                                >  3:157595/1‑149 (MQ=255)
                                                     tCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCCCCCCGCCAAGGCGCGGGGGCGGCTGGCCCCCGGGACTACCCTGGCCCGGGGCCGGGG‑CAGCCCg                                                                                            >  3:10808/1‑149 (MQ=255)
                                                     tCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCt                                                                                                       >  6:150672/1‑139 (MQ=255)
                                                     tCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCt                                                                                                       <  5:150672/139‑1 (MQ=255)
                                                       aGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAACAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCt                                                                                          <  8:134282/149‑1 (MQ=255)
                                                         ccGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCCTGTGGCACCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTgg                                                                                        <  7:202804/149‑1 (MQ=255)
                                                           gCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGc                                                                                      >  2:283165/1‑149 (MQ=255)
                                                                ggCATTGCCCTGGCGCATGCGGTGTGCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTggagacga                                                                                  <  3:19179/148‑7 (MQ=255)
                                                                ggCATTGCCCTGGCGCATGCGGGGGTCTGGCCAATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCgtgt                                                                                 <  8:319200/149‑1 (MQ=255)
                                                                   aTTGCCGTGGCGCATGCGGTGTGCTGGTCCATCA‑‑‑CCGCCTCGGTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTgtggt                                                                              <  5:84618/149‑1 (MQ=255)
                                                                      gCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGttgtgg                                                                            <  5:305762/148‑1 (MQ=255)
                                                                       cccTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTggg                                                                           <  1:67046/148‑1 (MQ=255)
                                                                         cTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGa                                                                        <  1:283165/149‑1 (MQ=255)
                                                                                   cGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTggg                                                               <  1:384766/148‑1 (MQ=255)
                                                                                    ggTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCACCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCAAGGTCCTGGG‑CATACCGCGGGGCCGTGGGGGGGGCGCGGCGCTGGGcg                                                             >  8:299680/1‑148 (MQ=255)
                                                                                      tgtTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTgg                                                           >  3:263056/1‑149 (MQ=255)
                                                                                       gtTCTGGTACATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCTCCGCCGGGCCAGCAGGCCAAGGCGCTGGGGCTGCTGGCCACCGGTACTACGCTGTCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGGGGTGGGCGAGGCACTGGGCTGGc                                                          <  6:242821/149‑1 (MQ=255)
                                                                                        ttCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGACCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTgtccataa                                                          >  5:238781/1‑141 (MQ=255)
                                                                                         tCTGGTCCATCATAGCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTgg                                                           <  1:328346/149‑1 (MQ=255)
                                                                                         tCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGt                                                        >  5:313030/1‑149 (MQ=255)
                                                                                         tCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAACGCGCTGGGCCTGCTGGCCACCGGTACTACACTGGCCATGGTCCTGGG‑CCTCCCGCAGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGt                                                        >  6:281439/1‑149 (MQ=255)
                                                                                          cTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTgg                                                                                                                 >  2:50361/1‑92 (MQ=255)
                                                                                          cTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTgg                                                                                                                 <  1:50361/92‑1 (MQ=255)
                                                                                             gTCCATCT‑‑‑CGGCCTCGCTGGGGGTGCGTGTGGCGGCGCCGGTCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATc                                                     <  4:52145/148‑1 (MQ=255)
                                                                                              tCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCAc                                                   <  8:318007/149‑1 (MQ=255)
                                                                                              gccctCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCAc                                                   <  4:20792/145‑1 (MQ=255)
                                                                                                 aTCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGtt                                                >  4:289481/1‑149 (MQ=255)
                                                                                                  tCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCCCCGGTACTACCCTGGCCATGGTCCTGGG‑CCTCCCGCTGGGCCGTGTGGTGGGCGAGGCACGGGGCTGGCGTATCACGTTc                                               >  5:89071/1‑149 (MQ=255)
                                                                                                           ccTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAg                                          >  2:277845/1‑148 (MQ=255)
                                                                                                               gCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTg                                     <  5:7250/149‑1 (MQ=255)
                                                                                                                   gCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCCAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                 <  4:10808/149‑1 (MQ=255)
                                                                                                                   gCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                 <  6:313030/149‑1 (MQ=255)
                                                                                                                    cGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCggg                                <  1:277845/149‑1 (MQ=255)
                                                                                                                           tgtgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcgc                         <  7:299680/149‑1 (MQ=255)
                                                                                                                             tgGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTa                       <  3:74851/149‑1 (MQ=255)
                                                                                                                                 gccgccGGGCCACCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTAATAACCTGGCCCTGGTCCAGGG‑CATCCCGCAGGGCCGGGTGGGGGGGGAGGGACTGGGGGGGCGGATCAAGTTCCTGAGCATTGCCGGGGTAGCGCTAGcca                   >  6:221629/1‑149 (MQ=255)
                                                                                                                                        ggCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CTTCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGTGGTAGCGCTAGCCACCATGCt            <  6:324609/149‑1 (MQ=255)
                                                                                                                                         gCCAGCAGGCCAAGGCGCTGGGCCTGCTGTCCACCGGTACTACCCTGGCCATGGTCCTGGGACATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGGTAGCCACCATGCt            >  1:111443/1‑149 (MQ=255)
                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc        <  3:289481/149‑1 (MQ=255)
                                                                                                                                            agcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc        >  5:310546/1‑149 (MQ=255)
                                                                                                                                             gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc       <  2:347582/149‑1 (MQ=255)
                                                                                                                                             gcagGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc       >  8:219461/1‑149 (MQ=255)
                                                                                                                                                ggCCAAGGCGCGGGGCCTGCGGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGa    <  6:212675/149‑1 (MQ=255)
                                                                                                                                                   cAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTCGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATc  >  8:35209/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATTGTCAGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCA‑‑‑CCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAGCAGGCCAAGGCGCTGGGCCTGCTGGCCACCGGTACTACCCTGGCCATGGTCCTGGG‑CATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATC  >  NC_002947/2775962‑2776253

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