Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,253,883 G→A 7.1% R251C (CGC→TGC)  PP_4630 ← MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,253,8830GA 7.1% 206.9 / 3.0 70R251C (CGC→TGC) PP_4630MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base A (3/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GACGGGCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAG  >  NC_002947/5253736‑5254030
                                                                                                                                                   |                                                                                                                                                     
gACGGGCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGgcgc                                                                                                                                                      >  5:272729/1‑149 (MQ=255)
 aCGGGCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCt                                                                                                                                                     >  2:264763/1‑149 (MQ=255)
     gCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGaa                                                                                                                                                 <  1:341388/149‑1 (MQ=255)
       cAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGt                                                                                                                                               >  6:106474/1‑149 (MQ=255)
            aCTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTg                                                                                                                                          >  2:131771/1‑149 (MQ=255)
                ccGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGc                                                                                                                                        >  7:339251/1‑147 (MQ=255)
                ccGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGa                                                                                                                                      <  2:27614/149‑1 (MQ=255)
                 cGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAg                                                                                                                                     <  1:360170/149‑1 (MQ=255)
                    aaGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCCGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTc                                                                                                                                 >  7:125312/1‑149 (MQ=255)
                      gcgcTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTctg                                                                                                                               <  8:20880/149‑1 (MQ=255)
                      gcgcTGGCGCAGCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCGGGCGAGTTctg                                                                                                                               <  6:21198/149‑1 (MQ=255)
                        gcTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTctgct                                                                                                                             <  6:354009/149‑1 (MQ=255)
                         cTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTctgct                                                                                                                             <  1:332650/148‑1 (MQ=255)
                         cTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTctgct                                                                                                                             >  2:332650/1‑148 (MQ=255)
                            gcgcATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCtgtg                                                                                                                          <  5:24096/148‑1 (MQ=255)
                              gcATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCATTAACAACGTGCATTCGAGCATCGACGCCAGGCGCTTGAAGTTCCTGGCGAGTTCTGCTGTGGGc                                                                                                                       <  2:117290/149‑1 (MQ=255)
                                   acacTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGAGCGACt                                                                                                                   >  4:127520/1‑148 (MQ=255)
                                    cacTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGa                                                                                                                 <  5:106474/149‑1 (MQ=255)
                                     acTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGAt                                                                                                                >  8:65289/1‑149 (MQ=255)
                                     acTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGAGCGACTGAt                                                                                                                <  3:127520/149‑1 (MQ=255)
                                      cTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATg                                                                                                               <  3:47506/149‑1 (MQ=255)
                                                 gcCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCa                                                                                                    <  1:131771/149‑1 (MQ=255)
                                                   ccAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAGCGGGCATTCCAGCATCGACGCCAGGCGCTCGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCAc                                                                                                   >  3:49513/1‑148 (MQ=255)
                                                         aggcTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCg                                                                                            >  2:279366/1‑149 (MQ=255)
                                                               cAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCgg                                                                                       <  8:125312/148‑1 (MQ=255)
                                                               cAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCgg                                                                                       <  8:339251/148‑1 (MQ=255)
                                                                 gggCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTGACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCtt                                                                                    >  5:306124/1‑149 (MQ=255)
                                                                  ggCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTc                                                                                   <  1:318114/149‑1 (MQ=255)
                                                                  ggCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTc                                                                                   >  1:220995/1‑149 (MQ=255)
                                                                     tGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCACTTCAAgc                                                                                                                                               <  3:83055/87‑2 (MQ=255)
                                                                     tGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCACTTCAAgc                                                                                                                                               >  4:83055/1‑86 (MQ=255)
                                                                        cagccagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGt                                                                             >  2:3513/1‑149 (MQ=255)
                                                                         agccagGCTGTCCTGTACCGCTTCCGCTGCTTCAAAAGCCAGTAACAACGGGCATTCCAGCATCGACGCCCGGCGCTTGAAGTGGCTGGCGAGTTCTGCTGTGGGCGACTGATGGGTAAACAGCACCAGGGCGTCCGGCTTCAGCCGtt                                                                            >  7:188192/1‑149 (MQ=255)
                                                                            cagGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCg                                                                          <  1:279366/149‑1 (MQ=255)
                                                                                cTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTGACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGAc                                                                     <  6:306124/149‑1 (MQ=255)
                                                                                    ccTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                                                 <  4:261314/149‑1 (MQ=255)
                                                                                     cTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAGCGGGCATTCCAGCATCGACGCCAGGCGCTCGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgc                                                                 <  4:49513/148‑1 (MQ=255)
                                                                                       gTACAGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAgcgc                                                               >  4:90650/1‑148 (MQ=255)
                                                                                        tACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCa                                                                        <  6:308472/138‑1 (MQ=255)
                                                                                        tACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACACCACCAGGGCGTCCGGCTTCAGCCGTTCGCa                                                                        >  5:308472/1‑138 (MQ=255)
                                                                                           cGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCt                                                          >  8:304112/1‑149 (MQ=255)
                                                                                                  gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTggg                                                                                                                        <  8:35383/80‑1 (MQ=255)
                                                                                                  gctgctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTggg                                                                                                                        >  7:35383/1‑80 (MQ=255)
                                                                                                      ctTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCg                                                <  7:65289/147‑1 (MQ=255)
                                                                                                        tCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTg                                            >  2:208759/1‑149 (MQ=255)
                                                                                                         cACCAGCCAGTAACAACGGGCATTCCAGCATCGCCGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCgg                                               >  6:167871/1‑145 (MQ=255)
                                                                                                         cACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCgg                                               <  5:167871/145‑1 (MQ=255)
                                                                                                         cACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCt                                             <  1:79330/147‑1 (MQ=255)
                                                                                                           ccagccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGccc                                         >  1:150000/1‑149 (MQ=255)
                                                                                                            cagccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCg                                        >  4:228690/1‑149 (MQ=255)
                                                                                                             agccagTAACAACGGGCATGCCAGCATCGACGCCAGGCGCTTGTAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGt                                       <  6:272729/149‑1 (MQ=255)
                                                                                                              gccagTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGtt                                      <  2:150000/149‑1 (MQ=255)
                                                                                                                   tAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCa                                 <  2:220995/149‑1 (MQ=255)
                                                                                                                        acGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCggg                            <  1:208759/149‑1 (MQ=255)
                                                                                                                          gggCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACTGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGAt                          >  2:356392/1‑149 (MQ=255)
                                                                                                                          gggCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCaa                                                             <  1:144292/115‑1 (MQ=255)
                                                                                                                          gggCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCaa                                                             >  2:144292/1‑115 (MQ=255)
                                                                                                                           ggCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATc                         <  3:228690/149‑1 (MQ=255)
                                                                                                                               ttCGAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGTCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCAcgc                     <  5:335651/149‑1 (MQ=255)
                                                                                                                               ttCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCAcgc                     <  7:304112/149‑1 (MQ=255)
                                                                                                                                tCCAGCATCGACGCCAGGCGCTTGAAGTGCGTGGGGAGTTGTGCTGTGGGCGACAGATGGCAAAACACCACCAGGGCGTCCGGCTTCAGCCGTTAGCACACCAGCGCCAGCAGCAGCAGCGGCAGCCCGTTAGCCAGCGGGCTcccccg                    >  2:342347/1‑144 (MQ=255)
                                                                                                                                tCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCAcgcg                    >  2:336983/1‑149 (MQ=255)
                                                                                                                                  cAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCgtg                  >  4:240341/1‑149 (MQ=255)
                                                                                                                                  cAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCgt                   <  1:3513/148‑1 (MQ=255)
                                                                                                                                         gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtc           <  1:336983/149‑1 (MQ=255)
                                                                                                                                         gACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTtcgtc           <  1:342347/149‑1 (MQ=255)
                                                                                                                                          tcGCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCgg                             >  6:56079/2‑130 (MQ=255)
                                                                                                                                          tcGCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCgg                             <  5:56079/129‑1 (MQ=255)
                                                                                                                                            gCCAGGCACTTCAAGCGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCCGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCt        >  2:187650/1‑149 (MQ=255)
                                                                                                                                             ccAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTg       >  8:880/1‑149 (MQ=255)
                                                                                                                                              cAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGCCCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGc      >  8:10109/1‑149 (MQ=255)
                                                                                                                                                ggggcTTGCAGTGGCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAATCAGCACGATGGCGTCGGGCTTGATCCGTTGGCAGAGCATCGGAAGCTGCTTCAGCGGCTTCCGGTTGGGCAGCGTGATCACGCTTGTTGGGTGGCTGccg    <  8:172910/146‑2 (MQ=255)
                                                                                                                                                 gcgcTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGcc     >  1:152519/1‑147 (MQ=255)
                                                                                                                                                  cgcTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCACGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCa   <  1:264763/149‑1 (MQ=255)
                                                                                                                                                  cgcTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAg  <  2:152519/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACGGGCCAGAAACTGCCGCAAGCGCTGGCGCATCACACTGCCCTGCGCGCCCAGGCAGGCTACAGGGCTGCCAGCCAGGCTGTCCTGTACCGCTTCCGCTGCTTCACCAGCCAGTAACAACGGGCATTCCAGCATCGACGCCAGGCGCTTGAAGTGCCTGGCGAGTTCTGCTGTGGGCGACTGATGGGAAAACAGCACCAGGGCGTCCGGCTTCAGCCGTTCGCAGACCAGCGCAAGCTCCTGCAGCGGCTGCCCGTTCGCCAGCGGGATCACGCGTGTTTCGTCGCTGCCCAG  >  NC_002947/5253736‑5254030

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