Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,681,667 A→T 6.4% E62D (GAA→GAT)  recJ → ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,681,6670AT 6.4% 223.7 / 4.4 78E62D (GAA→GAT) recJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base A (40/33);  new base T (3/2);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

CATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCT  >  NC_002947/1681519‑1681812
                                                                                                                                                    |                                                                                                                                                 
cATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGaa                                                                                                                                                   <  5:132624/149‑1 (MQ=255)
 aTTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTAGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAg                                                                                                                                                  <  5:137810/149‑1 (MQ=255)
      tGGGTAACCTGCCACCATTGCTGACCCGCCTTTACGGCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATc                                                                                                                                             <  7:352098/149‑1 (MQ=255)
        ggTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGa                                                                                                                                           >  2:220557/1‑149 (MQ=255)
            aCCTGCCAGCCTTGCTGACCCGCCTTTAGGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAcca                                                                                                                                        <  2:104646/148‑1 (MQ=255)
             ccTGCCAACCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAATGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagc                                                                                                                                      <  3:218731/149‑1 (MQ=255)
              cTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAcca                                                                                                                                        >  8:182477/1‑146 (MQ=255)
              cTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAcca                                                                                                                                        <  7:182477/146‑1 (MQ=255)
               tGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGCccagcgc                                                                                                                                    >  6:13967/1‑149 (MQ=255)
                 ccACGCTTGCTGACCCGCCTTTAGGCCGCGCGCGGCGTGCAGTGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCCGCTGGTCGATGGCTTCGAccagcag                                                                                                                                    <  4:112454/147‑3 (MQ=255)
                 ccACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgc                                                                                                                                    >  7:341800/1‑147 (MQ=255)
                 ccACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGAccagcgc                                                                                                                                    <  8:341800/147‑1 (MQ=255)
                       ttGCTGACGCGCCTTTACGGCGCGCGCGGCGTGCAGATCGAGGCCGAGCTTGACAATAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCATGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCATCGCAt                                                                                                                            <  6:289043/149‑1 (MQ=255)
                         gCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATc                                                                                                                           <  2:337735/148‑1 (MQ=255)
                          cTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCt                                                                                                                         >  1:342679/1‑149 (MQ=255)
                                 gCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgtc                                                                                                                  >  3:261478/1‑149 (MQ=255)
                                  ccTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtcgtcg                                                                                                                 >  8:2184/1‑149 (MQ=255)
                                   cTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg                                                                                                                <  1:38381/149‑1 (MQ=255)
                                        cgccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACt                                                                                                           <  4:130190/149‑1 (MQ=255)
                                        cgccgcGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGACGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCAAc                                                                                                            >  5:148118/1‑148 (MQ=255)
                                           cgcgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTc                                                                                                         >  3:17325/1‑148 (MQ=255)
                                             cgcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGAt                                                                                                      <  1:220557/149‑1 (MQ=255)
                                              gcgcgGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATg                                                                                                     <  2:373769/149‑1 (MQ=255)
                                                gcggcgTGCAGATCGAGGACGACCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCCAGGGCATCGAGGCCGGTGTGGACCTTCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGcc                                                                                                   <  2:114655/149‑1 (MQ=255)
                                                      tGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                                                                                             <  7:2184/149‑1 (MQ=255)
                                                       gCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg                                                                                                                >  1:376685/1‑129 (MQ=255)
                                                       gCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCgg                                                                                                                <  2:376685/129‑1 (MQ=255)
                                                       gCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGAGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTg                                                                                            >  3:248246/1‑149 (MQ=255)
                                                         agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTg                                                                                            <  2:147866/147‑1 (MQ=255)
                                                         agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTg                                                                                            >  1:147866/1‑147 (MQ=255)
                                                         agagCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGc                                                                                           >  6:70160/1‑148 (MQ=255)
                                                                 gCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAg                                                                                  >  5:67266/1‑149 (MQ=255)
                                                                     aGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCAcc                                                                              <  8:226411/149‑1 (MQ=255)
                                                                      gCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGAGTTCGATGCCGACGGTGCCCCCGCCCGCCCCg                                                                             >  6:252519/1‑149 (MQ=255)
                                                                        tGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTc                                                                           >  7:116762/1‑149 (MQ=255)
                                                                          gACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGCTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCgc                                                                         >  3:337496/1‑148 (MQ=255)
                                                                           aCAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                        >  1:314730/1‑149 (MQ=255)
                                                                           aCAAGAGCCGGGCAGGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCggg                                                                        <  5:252519/149‑1 (MQ=255)
                                                                             aaGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGt                                                                      >  5:234919/1‑149 (MQ=255)
                                                                                agCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCt                                                                   >  7:348146/1‑149 (MQ=255)
                                                                                  ccTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTgg                                                                 <  2:342679/149‑1 (MQ=255)
                                                                                        aCGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgc                                                           >  7:63507/1‑149 (MQ=255)
                                                                                        aCGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgc                                                           >  5:304829/1‑149 (MQ=255)
                                                                                         cGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcg                                                          >  6:380628/1‑149 (MQ=255)
                                                                                         cGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcg                                                          <  4:17325/149‑1 (MQ=255)
                                                                                            ctgctgCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGACACCGCCAGCACCGTCGGGGTGCTGGgcctgcgcct                                                       >  4:166366/1‑149 (MQ=255)
                                                                                                  cccTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTggg                                                 >  4:262103/1‑149 (MQ=255)
                                                                                                    cTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTg                                               >  3:39805/1‑149 (MQ=255)
                                                                                                       ccAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCgg                                            >  2:79504/1‑149 (MQ=255)
                                                                                                         agcagcTCAAGGGCATCGAGGCCGGGGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGcc                                          <  8:63507/149‑1 (MQ=255)
                                                                                                          gcagcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCa                                                                                   >  6:268104/1‑107 (MQ=255)
                                                                                                          gcagcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCa                                                                                   <  5:268104/107‑1 (MQ=255)
                                                                                                              ataaaGGGCTTCGAGGCGGGTGTGGACCTGCTGGGCGAAGCCATCGACCTGCGCCTGCGAATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGt                                     <  4:248246/146‑1 (MQ=255)
                                                                                                               tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTc                                    >  8:159689/1‑149 (MQ=255)
                                                                                                                 aaGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGa                                  >  7:275360/1‑149 (MQ=255)
                                                                                                                 aaGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACACGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGa                                  >  7:270363/1‑149 (MQ=255)
                                                                                                                   gggCATCGAGGCCGGTGTGGCCCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGCCTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGAt                                 >  6:110944/1‑148 (MQ=255)
                                                                                                                   gggCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGAt                                 >  7:333867/1‑148 (MQ=255)
                                                                                                                      cATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCTGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTAcc                             <  1:79504/149‑1 (MQ=255)
                                                                                                                       aTCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCt                            >  8:102010/1‑149 (MQ=255)
                                                                                                                        tCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTg                           <  5:380628/149‑1 (MQ=255)
                                                                                                                        tCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTg                           <  6:67266/149‑1 (MQ=255)
                                                                                                                        tCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTg                           <  4:39805/149‑1 (MQ=255)
                                                                                                                         cGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtc                                                                                                                     <  8:337168/58‑1 (MQ=255)
                                                                                                                         cGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCAtc                                                                                                                     >  7:337168/1‑58 (MQ=255)
                                                                                                                         cGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGCGTCCTGGGCCTGCGCCTGCGGGGTGCGGCACAAGACGATTACCAgg                          >  6:186175/1‑149 (MQ=255)
                                                                                                                               cGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGccc                     >  4:208307/1‑148 (MQ=255)
                                                                                                                                       tggcgctgGTCGATGGCTTCGACCAGCACCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTc            >  2:168583/5‑149 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                <  4:349808/113‑1 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGt                                                >  3:349808/1‑113 (MQ=255)
                                                                                                                                       aCCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTc            >  6:319210/1‑149 (MQ=255)
                                                                                                                                        ccTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTc            >  4:211802/1‑148 (MQ=255)
                                                                                                                                         gcgctgGTCGATGGCTTCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTgg                                                  >  1:353335/3‑109 (MQ=255)
                                                                                                                                         atgctgGTCGAAGCCATCGACCAGCGCCAGCGCATGCTCATCGTCGGTGTCTTCGATGGCGACGGTGCCACCGCCAGCACGGTCGGTGTGCTGGGCCTGAGCCTGCTGGGTGCTGCCGATGTCGATTACCTGGTTCCCAACCGCTTCGa          <  4:337496/148‑1 (MQ=255)
                                                                                                                                           gctgGTCGATGGCTTCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGt                         <  2:237318/132‑1 (MQ=255)
                                                                                                                                           gctgGTCGATGGCTTCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGt                         >  1:237318/1‑132 (MQ=255)
                                                                                                                                             tggagttaGCCATCCACCAGGGCAAGCGCATCCTCATCTTCGGCGACGTCGAGGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGTCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTAt      <  8:247034/142‑1 (MQ=255)
                                                                                                                                             tgGTCGAAGGCAGCGACCAGCGCCAGGGCATCCTCATCGTGGGCGACTTCGATGCCGAGGGTGCCACCGCCAGCACCGTCGGGGTGCGGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTAt      <  5:186175/149‑1 (MQ=255)
                                                                                                                                               gTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTCCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTAt      <  3:262103/147‑1 (MQ=255)
                                                                                                                                                  gAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGCCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCt  <  6:304829/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CATTTCTGGGTAACCTGCCACCCTTGCTGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGGCT  >  NC_002947/1681519‑1681812

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