Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 641,388 A→C 5.1% F155C (TTC→TGC) ‡ bdhA ← 2,3‑butanediol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947641,3880AC 5.1% 252.1 / 4.1 79F155C (TTC→TGC) ‡bdhA2,3‑butanediol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base C (2/2);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

CTGCAGCCTTGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCAGTAATAGCAGGTA  >  NC_002947/641242‑641522
                                                                                                                                                  |                                                                                                                                      
cTGCAGCCTTGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAg                                                                                                                                      <  5:312752/149‑1 (MQ=255)
 tGCAGCCTTGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGc                                                                                                                                     >  3:288900/1‑149 (MQ=255)
   cAGCCTTGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGccc                                                                                                                                   >  3:259444/1‑149 (MQ=255)
   cAGCCTTGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGcc                                                                                                                                    <  4:350796/148‑1 (MQ=255)
       cTTGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCg                                                                                                                               <  6:61702/149‑1 (MQ=255)
        ttGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGACTTCGGCAGGGAAGCCCTGCgg                                                                                                                              >  6:187914/1‑149 (MQ=255)
          ggCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCgg                                                                                                                              <  4:46390/147‑1 (MQ=255)
          ggCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTGCCGCAGGGCTTCCCTGccg                                                                                                                              >  2:195956/1‑145 (MQ=255)
           gcgcACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAg                                                                                                                           >  8:309666/1‑149 (MQ=255)
             gcACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGcc                                                                                                                         >  2:17237/1‑149 (MQ=255)
              cacaTGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCg                                                                                                                        <  5:152831/149‑1 (MQ=255)
               acaTGATGGTGTACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCCCCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCgg                                                                                                                       >  2:167999/1‑149 (MQ=255)
               acaTGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTGCCGCAGGGCTTCCCTGCCGAAGCCgg                                                                                                                       <  1:195956/149‑1 (MQ=255)
                caTGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAAAGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTGCCGCAGGGCTTCCCTGCCGAAGCCggg                                                                                                                      <  6:356474/149‑2 (MQ=255)
                 atgatgGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTa                                                                                                                     <  4:251955/149‑1 (MQ=255)
                 atgatgGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTa                                                                                                                     <  8:333603/149‑1 (MQ=255)
                       gTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCa                                                                                                                <  1:17237/148‑1 (MQ=255)
                       gTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAg                                                                                                               >  1:122604/1‑149 (MQ=255)
                          cacaGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGt                                                                                                             <  1:164838/148‑1 (MQ=255)
                             atgCCGATGGTGCCGGCCCCAACGACTACAACGGTTTGCCCTAGGAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTTCCCCGTCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGc                                                                                                         <  5:15886/147‑1 (MQ=255)
                               gCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCgg                                                                                                       >  6:58258/1‑149 (MQ=255)
                                  gATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCa                                                                                                    <  4:288900/149‑1 (MQ=255)
                                  gATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCa                                                                                                    >  2:115012/1‑149 (MQ=255)
                                   aTGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCt                                                                                                                                  <  1:253606/118‑1 (MQ=255)
                                   aTGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCt                                                                                                                                  >  2:253606/1‑118 (MQ=255)
                                    tGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACg                                                                                                  <  1:115012/149‑1 (MQ=255)
                                      gTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGtt                                                                                                <  3:327123/149‑1 (MQ=255)
                                         ccGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGAc                                                                                             <  6:288159/149‑1 (MQ=255)
                                         ccGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGACTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGAc                                                                                             <  5:187914/149‑1 (MQ=255)
                                            gCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGtt                                                                                                <  7:319084/143‑1 (MQ=255)
                                             ccccAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGc                                                                                         >  8:333477/1‑149 (MQ=255)
                                              cccAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCt                                                                                        >  4:170938/1‑149 (MQ=255)
                                                cAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCg                                                                                      >  6:269112/1‑149 (MQ=255)
                                                 aaCCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGTGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCgg                                                                                     >  7:249122/1‑149 (MQ=255)
                                                 aaCCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCgg                                                                                     >  8:131389/1‑149 (MQ=255)
                                                    cACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCg                                                                                   >  7:95385/1‑148 (MQ=255)
                                                     aCTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGaa                                                                                 <  7:196716/149‑1 (MQ=255)
                                                       tACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAgg                                                                               >  5:348914/1‑149 (MQ=255)
                                                         cAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTaga                                                                                                                   >  8:216386/1‑111 (MQ=255)
                                                         cAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTaga                                                                                                                   <  7:216386/111‑1 (MQ=255)
                                                         cAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCa                                                                             >  4:26803/1‑149 (MQ=255)
                                                          aaCGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCAc                                                                            >  5:216409/1‑149 (MQ=255)
                                                                  gCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCAc                                                                                                   >  7:49841/1‑118 (MQ=255)
                                                                  gCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCAc                                                                                                   <  8:49841/118‑1 (MQ=255)
                                                                    ccAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTtcat                                                                  >  1:314775/1‑149 (MQ=255)
                                                                    ccAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTtcat                                                                  <  5:161884/149‑1 (MQ=255)
                                                                      aaGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTtcatca                                                                <  5:161915/149‑1 (MQ=255)
                                                                         cAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGTGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGc                                                             <  8:249122/149‑1 (MQ=255)
                                                                          aGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCgttg                                                                       >  7:311282/1‑138 (MQ=255)
                                                                          aGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCgttg                                                                       <  8:311282/138‑1 (MQ=255)
                                                                                ccGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACATTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGaa                                                      >  1:2121/1‑149 (MQ=255)
                                                                                 cGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTtcatc                                                                 >  4:1271/1‑137 (MQ=255)
                                                                                 cGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTtcatc                                                                 <  3:1271/137‑1 (MQ=255)
                                                                                 cGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAAc                                                     <  1:40749/149‑1 (MQ=255)
                                                                                   gCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGc                                                   >  8:159658/1‑149 (MQ=255)
                                                                                     cTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGc                                                   >  3:121777/1‑147 (MQ=255)
                                                                                        tttCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAg                                                <  7:131389/147‑1 (MQ=255)
                                                                                        tttCACCGAGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCATCCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCAGGTGAACGCCAGGc                                              >  7:142680/1‑149 (MQ=255)
                                                                                            aCCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTc                                          <  6:348914/149‑1 (MQ=255)
                                                                                                 gTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGa                                     <  4:259444/149‑1 (MQ=255)
                                                                                                  tGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGAt                                    <  5:58258/149‑1 (MQ=255)
                                                                                                   gCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGAtg                                   >  2:300711/1‑149 (MQ=255)
                                                                                                     aTACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGAtgtt                                 <  8:95385/149‑1 (MQ=255)
                                                                                                     aTACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACATTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGAtgtt                                 <  2:2121/149‑1 (MQ=255)
                                                                                                          cACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTa                              >  6:217188/1‑147 (MQ=255)
                                                                                                            ccGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAgg                          <  2:314775/149‑1 (MQ=255)
                                                                                                             cGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGTGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTtcat                                                                  <  5:258732/108‑1 (MQ=255)
                                                                                                             cGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTtcat                                                                  >  6:258732/1‑108 (MQ=255)
                                                                                                                cAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCAt                      <  3:170938/149‑1 (MQ=255)
                                                                                                                  gCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATgg                    <  1:300711/149‑1 (MQ=255)
                                                                                                                     gCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGTCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACTCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGt                  >  8:200772/1‑148 (MQ=255)
                                                                                                                     gCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTg                 <  5:217188/149‑1 (MQ=255)
                                                                                                                         gATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCa               <  7:159658/147‑1 (MQ=255)
                                                                                                                         gATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCa               <  3:26803/147‑1 (MQ=255)
                                                                                                                            cAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCAGTaa           >  5:39572/1‑148 (MQ=255)
                                                                                                                                 ccccGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCAGTAATAGCa      <  7:333477/148‑1 (MQ=255)
                                                                                                                                 ccccGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCAGTAATAGCAg     >  6:24346/1‑149 (MQ=255)
                                                                                                                                 ccccGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCAGTAATAGCAg     <  5:269112/149‑1 (MQ=255)
                                                                                                                                    cGGCTGCCGCAGGGCTTCCCTGCCGAAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCAGTAATAGCAGGTa  >  4:121051/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                      
CTGCAGCCTTGGCGCACATGATGGTGCACAGGCCGATGGTGCCGGCCCCAACCACTACAACGGTTTGCCCAAGCAGGCTGCCGGCCTTTTTCACCGCGTGCATACCCACCGCCAGCGGCTCGATCAGTGCCCCGGCTTCGGCAGGGAAGCCCTGCGGCAGCCGGTAGAGCAGGTTGGCGGGCACGTTGACCAGCTCGGCGAAGGCACCGTTGTTCATCAGGCCGGTGAACGCCAGGCGTTCGCAGATGTTGTACAGGCCATGGGTGCAGTAATAGCAGGTA  >  NC_002947/641242‑641522

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