Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,984,041 C→T 7.5% intergenic (+35/+155) ilvE → / ← PP_3512 branched‑chain‑amino acid aminotransferase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,984,0410CT 7.5% 194.6 / 4.2 66intergenic (+35/+155)ilvE/PP_3512branched‑chain‑amino acid aminotransferase/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (2/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

CCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGT  >  NC_002947/3983895‑3984188
                                                                                                                                                  |                                                                                                                                                   
ccACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCaa                                                                                                                                                   >  6:155882/1‑149 (MQ=255)
  aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                                 >  1:166077/1‑149 (MQ=255)
    gACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGc                                                                                                                                               <  5:350471/149‑1 (MQ=255)
         ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGcc                                                                                                                                          >  8:241189/1‑149 (MQ=255)
         ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGcc                                                                                                                                          >  4:326787/1‑149 (MQ=255)
          gAAAAGGTCGGCCCGGTTAGCGAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGAGATGTAGAGGCACCGGCAGGGTGGATCGTCAAGGTCGCCTGAAACGCATCAGGTGTATGATAACGGGGCATGCCAGCAATGGCATGc                                                                                                                                           <  1:55134/147‑1 (MQ=255)
          gAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGcc                                                                                                                                          <  3:349308/148‑1 (MQ=255)
          gAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGccc                                                                                                                                         >  3:226147/1‑149 (MQ=255)
                   ggCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGtttttt                                                                                                                                 >  5:246172/1‑148 (MQ=255)
                    gCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCa                                                                                                                               <  5:414659/149‑1 (MQ=255)
                     cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAt                                                                                                                              <  1:263909/149‑1 (MQ=255)
                          ttACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGtt                                                                                                                                     <  6:190330/137‑1 (MQ=255)
                          ttACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGtt                                                                                                                                     >  5:190330/1‑137 (MQ=255)
                             cccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                                 <  4:384649/122‑1 (MQ=255)
                             cccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                                 >  3:384649/1‑122 (MQ=255)
                             cccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTac                                                                                                                      <  3:178323/149‑1 (MQ=255)
                              ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                     <  7:208798/149‑1 (MQ=255)
                              ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                     >  7:222373/1‑149 (MQ=255)
                                        tACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTc                                                                                                           >  6:417168/1‑149 (MQ=255)
                                         aCAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                          >  6:291839/1‑149 (MQ=255)
                                          cAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGt                                                                                                         <  5:155882/149‑1 (MQ=255)
                                                  tGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGCGCgg                                                                                                 >  8:95821/1‑149 (MQ=255)
                                                    aCCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCggg                                                                                                <  7:241189/148‑1 (MQ=255)
                                                    aCCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCgg                                                                                                 <  3:326787/147‑1 (MQ=255)
                                                      cGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGtt                                                                                              >  4:151573/1‑148 (MQ=255)
                                                      cGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGtt                                                                                              >  1:239470/1‑148 (MQ=255)
                                                      cGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGtt                                                                                              >  2:134662/1‑148 (MQ=255)
                                                        gTATCCACAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCAGCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATCGCCCGATTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTAc                                                                                           >  4:347473/1‑149 (MQ=255)
                                                             cAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgc                                                                                       >  5:11582/1‑148 (MQ=255)
                                                              agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgcg                                                                                      >  8:196555/1‑148 (MQ=255)
                                                              agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTAAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgcg                                                                                      >  1:298711/1‑148 (MQ=255)
                                                                agCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAAt                                                                                   <  4:226147/149‑1 (MQ=255)
                                                                  cggcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAg                                                                                 >  6:310432/1‑149 (MQ=255)
                                                                    gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                               <  1:134662/149‑1 (MQ=255)
                                                                    gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                               <  6:246172/149‑1 (MQ=255)
                                                                    gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                               <  6:11582/149‑1 (MQ=255)
                                                                    gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                               <  3:151573/149‑1 (MQ=255)
                                                                             agGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTAATGTGTTTACACGGCCCCCTCGTAGGAGCTGGTTTACCCGCGAATAGGCCACACCCgc                                                                      >  8:130323/1‑149 (MQ=255)
                                                                                     gCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaa                                                               >  2:345478/1‑148 (MQ=255)
                                                                                       aGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCATCCaataa                                                            <  5:291839/149‑1 (MQ=255)
                                                                                        ggCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaataat                                                           <  2:298711/149‑1 (MQ=255)
                                                                                         gCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATg                                                          <  7:196555/149‑1 (MQ=255)
                                                                                         gCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATg                                                          <  6:208424/149‑1 (MQ=255)
                                                                                                 gTCAAGGTCGCCTGAAACGCATGAGCTGTATGAAAACGGGGCATGCGAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGTCCACACCCTCCAGCCAATAATGCCTGGga                                                   <  5:200095/148‑1 (MQ=255)
                                                                                                 gTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGgag                                                  <  5:417168/149‑1 (MQ=255)
                                                                                                       gTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCAcc                                            <  8:27695/149‑1 (MQ=255)
                                                                                                           ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTg                                        >  6:76582/1‑149 (MQ=255)
                                                                                                                   gCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTg                                  >  3:255158/1‑147 (MQ=255)
                                                                                                                     aTCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCCGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGc                              >  6:123764/1‑149 (MQ=255)
                                                                                                                             tATGAAAACGGGGCATGCCAGCAATGGCGTGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAg                      >  1:24109/1‑149 (MQ=255)
                                                                                                                             tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                   >  8:275172/1‑120 (MQ=255)
                                                                                                                             tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                   <  7:275172/120‑1 (MQ=255)
                                                                                                                             tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGtt                                         >  1:61272/1‑130 (MQ=255)
                                                                                                                             tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGtt                                         <  2:61272/130‑1 (MQ=255)
                                                                                                                              aTGAGAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTa                                                                                            >  6:168477/1‑78 (MQ=255)
                                                                                                                              aTGACAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTa                                                                                            <  5:168477/78‑1 (MQ=255)
                                                                                                                               tGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGc                     <  2:166077/148‑1 (MQ=255)
                                                                                                                                 aaaaCGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                               >  4:376186/1‑88 (MQ=37)
                                                                                                                                 aaaaCGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                               <  3:376186/88‑1 (MQ=37)
                                                                                                                                 aaaaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                   <  4:255158/148‑1 (MQ=255)
                                                                                                                                 aaaaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                   <  5:76582/148‑1 (MQ=255)
                                                                                                                                   aaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGc                 >  6:74528/1‑148 (MQ=255)
                                                                                                                                        ggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc           <  2:7998/149‑1 (MQ=255)
                                                                                                                                             gCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGcc      >  5:166869/1‑149 (MQ=255)
                                                                                                                                                 gCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGt  <  1:345478/149‑1 (MQ=255)
                                                                                                                                                 gCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGt  >  5:318857/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGT  >  NC_002947/3983895‑3984188

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: ‑