Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,235,038 A→G 6.6% N50S (AAC→AGC)  uvrB → uvrABC excinuclease protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,235,0380AG 6.6% 175.6 / 2.1 60N50S (AAC→AGC) uvrBuvrABC excinuclease protein B
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base G (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTCCGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTACGACTACTACCAGCCGGAAG  >  NC_002947/2234892‑2235184
                                                                                                                                                  |                                                                                                                                                  
gTCCGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACg                                                                                                                                                  >  6:189530/1‑149 (MQ=255)
 tCCGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGt                                                                                                                                                 >  8:50327/1‑149 (MQ=255)
 tCCGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGt                                                                                                                                                 >  7:211011/1‑149 (MQ=255)
   cGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAGCGTc                                                                                                                                                >  7:7635/1‑148 (MQ=255)
   cGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAGCGTCa                                                                                                                                               >  7:3166/1‑149 (MQ=255)
   cGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCa                                                                                                                                               <  3:272187/149‑1 (MQ=255)
     aGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATc                                                                                                                                             <  2:7302/149‑1 (MQ=255)
     aGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCAGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCGCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATc                                                                                                                                             >  6:284397/1‑149 (MQ=255)
                     acGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGcc                                                                                                                             >  3:253577/1‑149 (MQ=255)
                      cGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGccc                                                                                                                            <  5:284397/149‑1 (MQ=255)
                         ttCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGcc                                                                                                                             <  6:99741/145‑1 (MQ=255)
                         ttCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGcc                                                                                                                             >  5:99741/1‑145 (MQ=255)
                            cAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCt                                                                                                                      >  1:187167/1‑149 (MQ=255)
                                cGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTg                                                                                                                  >  7:263699/1‑149 (MQ=255)
                                  gCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGACGGCTTCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCAATGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCt                                                                                                                <  7:171697/149‑1 (MQ=255)
                                   cTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTg                                                                                                               <  3:102519/149‑1 (MQ=255)
                                      gCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGgcg                                                                                                            >  1:291323/1‑149 (MQ=255)
                                           cAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGcc                                                                                                                             >  5:257375/1‑127 (MQ=255)
                                           cAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGcc                                                                                                                             <  6:257375/127‑1 (MQ=255)
                                           cAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCcaa                                                                                                       <  4:67399/149‑1 (MQ=255)
                                            aGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCcaa                                                                                                       >  4:89583/1‑148 (MQ=255)
                                                cGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGa                                                                                                  <  2:187167/149‑1 (MQ=255)
                                                           gCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGgccgcc                                                                                       >  7:252357/1‑149 (MQ=255)
                                                                    tCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTAc                                                                              >  3:58593/1‑149 (MQ=255)
                                                                      gAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACg                                                                             <  5:189530/148‑1 (MQ=255)
                                                                          gCATCGAGGCTGGTCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTGCGGGCAAGACCTTCATCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGTGAg                                                                        <  4:88285/149‑1 (MQ=255)
                                                                           cATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGt                                                                       <  3:89583/149‑1 (MQ=255)
                                                                            aTCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGt                                                                                   >  1:126565/1‑136 (MQ=255)
                                                                               gAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCaa                                                                   <  4:58593/149‑1 (MQ=255)
                                                                                        cTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCAttctt                                                          <  8:252357/149‑1 (MQ=255)
                                                                                             gCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTcccc                                                      <  4:253577/148‑1 (MQ=255)
                                                                                             gCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGAACTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCTGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTcccc                                                      <  8:263699/148‑1 (MQ=255)
                                                                                                  aGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTcc                                                        >  5:257030/1‑141 (MQ=255)
                                                                                                  aGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTcc                                                        <  6:257030/141‑1 (MQ=255)
                                                                                                         gCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAGCGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTg                                          <  8:7635/148‑1 (MQ=255)
                                                                                                         gCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAGCGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTg                                          <  8:3166/148‑1 (MQ=255)
                                                                                                          cTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTg                                          >  7:34130/1‑147 (MQ=255)
                                                                                                             ggggTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTa                                     <  2:291323/149‑1 (MQ=255)
                                                                                                                   aCCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCAGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGt                               >  5:216829/1‑149 (MQ=255)
                                                                                                                       gTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTcc                           >  2:186459/1‑149 (MQ=255)
                                                                                                                       gTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTcc                           >  1:168556/1‑149 (MQ=255)
                                                                                                                        ttCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCt                          <  8:211011/149‑1 (MQ=255)
                                                                                                                          cGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTa                         <  7:50327/148‑1 (MQ=255)
                                                                                                                            ggCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTa                      <  5:12655/149‑1 (MQ=255)
                                                                                                                             gCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCCACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTAc                     >  6:178291/1‑149 (MQ=255)
                                                                                                                                  aCCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGa                                        >  7:186816/1‑125 (MQ=255)
                                                                                                                                  aCCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGa                                        <  8:186816/125‑1 (MQ=255)
                                                                                                                                  aCCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTAc                     <  7:208471/144‑1 (MQ=255)
                                                                                                                                  aCCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTAc                     >  8:208471/1‑144 (MQ=255)
                                                                                                                                    cTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTACGactac               >  1:200257/1‑148 (MQ=255)
                                                                                                                                     ttCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTa                      <  6:128173/140‑1 (MQ=255)
                                                                                                                                     ttCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTa                      >  5:128173/1‑140 (MQ=255)
                                                                                                                                     ttCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTAc                     >  3:158483/1‑141 (MQ=255)
                                                                                                                                     ttCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTAc                     <  4:158483/141‑1 (MQ=255)
                                                                                                                                     ttCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTACGActacta             >  3:39748/1‑149 (MQ=255)
                                                                                                                                          cATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCAttcttc                                                         >  1:186356/1‑100 (MQ=255)
                                                                                                                                          cATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCAttcttc                                                         <  2:186356/100‑1 (MQ=255)
                                                                                                                                            ttGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTACGACTACTACCAGCCg      >  6:69330/1‑149 (MQ=255)
                                                                                                                                               ccAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTACGACTACTACCAGCCGGa    >  1:113798/1‑148 (MQ=255)
                                                                                                                                                cAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTACGACTACTACCAGCCGGAAg  >  6:112145/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTCCGAGTTCCAGCTCGTCACACGTTTCCAGCCGGCTGGCGACCAGCCCGAGGCTATCCGCCAGTTGGTCGAAGGCATCGAGGCGGGGCTGTCGCACCAGACGTTGCTCGGGGTGACCGGTTCGGGCAAGACCTTCAGCATTGCCAACGTCATCCAGCAGGTGCAGCGCCCGACGCTGGTGCTGGCGCCCAACAAGACGCTGGCCGCCCAGTTGTACGGGGAGTTCAAGGCATTCTTCCCCAACAACGCGGTGGAGTACTTCGTTTCCTATTACGACTACTACCAGCCGGAAG  >  NC_002947/2234892‑2235184

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: