Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,417,395 A→G 6.0% K215E (AAG→GAG)  PP_1239 → metallo‑beta‑lactamase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,417,3950AG 6.0% 196.7 / 2.8 67K215E (AAG→GAG) PP_1239metallo‑beta‑lactamase family protein
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base G (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCAACACCTCAAGCAAGCGGAGC  >  NC_002947/1417247‑1417535
                                                                                                                                                    |                                                                                                                                            
ccTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAgcagca                                                                                                                                              <  6:361720/149‑1 (MQ=255)
    ctgGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAg                                                                                                                                            <  7:211196/147‑1 (MQ=255)
       gCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGaacaa                                                                                                                                       >  3:54652/1‑149 (MQ=255)
         aCGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAAcc                                                                                                                                     >  1:163972/1‑149 (MQ=255)
               tCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCca                                                                                                                                >  6:130033/1‑148 (MQ=255)
                 aCTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCacac                                                                                                                             >  4:7798/1‑149 (MQ=255)
                 aCTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCacac                                                                                                                             >  2:137112/1‑149 (MQ=255)
                   tACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACtt                                                                                                                           <  4:35844/149‑1 (MQ=255)
                   tACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACtt                                                                                                                           >  2:36409/1‑149 (MQ=255)
                      ccGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCaa                                                                                                                                             >  4:311807/1‑128 (MQ=255)
                      ccGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCaa                                                                                                                                             <  3:311807/128‑1 (MQ=255)
                      ccGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGc                                                                                                                        >  8:152788/1‑149 (MQ=255)
                      ccGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGc                                                                                                                        >  3:47989/1‑149 (MQ=255)
                             ttCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCa                                                                                                                  <  4:54652/148‑1 (MQ=255)
                                tGAAGCAGCGGGTGGGTGGCATGAAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGc                                                                                                              >  3:117981/1‑149 (MQ=255)
                                    gcagcGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTg                                                                                                          >  7:298478/1‑149 (MQ=255)
                                        cGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGa                                                                                                      <  1:137112/149‑1 (MQ=255)
                                          ggtgggtgGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTg                                                                                                                          <  2:185315/127‑1 (MQ=255)
                                          ggtgggtgGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTg                                                                                                                          >  1:185315/1‑127 (MQ=255)
                                           gtgggtggCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCt                                                                                                             >  3:277296/1‑139 (MQ=255)
                                           gtgggtggCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCt                                                                                                             <  4:277296/139‑1 (MQ=255)
                                           gtgggtggCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGAGCTTTCTCCacac                                                                                                   >  8:253170/1‑149 (MQ=255)
                                                tggCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCg                                                                                              <  7:152788/149‑1 (MQ=255)
                                                  gCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCggg                                                                                            <  4:117981/149‑1 (MQ=255)
                                                    aTGCAAGGGCATTTGAGCAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGgtg                                                                                          >  3:45692/1‑149 (MQ=255)
                                                       cAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGa                                                                                       <  2:43362/149‑1 (MQ=255)
                                                        aaGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGAc                                                                                      <  1:36409/149‑1 (MQ=255)
                                                        aaGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGAc                                                                                      >  3:272614/1‑149 (MQ=255)
                                                              aTTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGAc                                                                                >  5:343372/1‑149 (MQ=255)
                                                                 tGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTg                                                                              <  6:189772/148‑1 (MQ=255)
                                                                    acaacCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTc                                                                          >  4:54585/1‑149 (MQ=255)
                                                                      aacCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACACCCTGCCACACTTGGCCCGCCAGTGCTTTGTCGCCCCCCGCGGGTGTGACCCGGACTGGCtgcg                                                                        >  4:250527/1‑146 (MQ=255)
                                                                                 tGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCGACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATc                                                             <  7:253170/149‑1 (MQ=255)
                                                                                   cgcgcCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAg                                                           <  5:130033/149‑1 (MQ=255)
                                                                                         tGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAacac                                                     >  6:241444/1‑149 (MQ=255)
                                                                                           gTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGa                                                                                 <  7:243808/119‑1 (MQ=255)
                                                                                           gTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGa                                                                                 >  8:243808/1‑119 (MQ=255)
                                                                                             gCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGc                                                 <  4:47989/149‑1 (MQ=255)
                                                                                                   gTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATc                                                             >  6:55552/1‑131 (MQ=255)
                                                                                                   gTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATc                                                             <  5:55552/131‑1 (MQ=255)
                                                                                                       ggCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCATCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGgc                                        <  7:210328/148‑1 (MQ=255)
                                                                                                          tGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCGAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGc                                                                  <  1:324993/119‑1 (MQ=255)
                                                                                                          tGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCGAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGc                                                                  <  1:319645/119‑1 (MQ=255)
                                                                                                          tGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCGAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGc                                                                  >  2:319645/1‑119 (MQ=255)
                                                                                                          tGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCGAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGc                                                                  >  2:324993/1‑119 (MQ=255)
                                                                                                           ggAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGa                                    >  4:1583/1‑148 (MQ=255)
                                                                                                               caacCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCg                               >  4:97491/1‑149 (MQ=255)
                                                                                                                 acCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGcc                             <  4:45692/149‑1 (MQ=255)
                                                                                                                  ccTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGcc                             >  4:287688/1‑148 (MQ=255)
                                                                                                                   cTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTa                           <  4:272614/149‑1 (MQ=255)
                                                                                                                     gCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGc                         >  6:7527/1‑149 (MQ=255)
                                                                                                                      cAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTa                           >  5:154006/1‑146 (MQ=255)
                                                                                                                      cAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTa                           <  6:154006/146‑1 (MQ=255)
                                                                                                                      cAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGcc                        >  4:167659/1‑149 (MQ=255)
                                                                                                                      cAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGc                         <  2:163972/148‑1 (MQ=255)
                                                                                                                       aaCACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTa                           >  4:5945/1‑145 (MQ=255)
                                                                                                                       aaCACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTa                           <  3:5945/145‑1 (MQ=255)
                                                                                                                         cacCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCaa                     <  8:298478/149‑1 (MQ=255)
                                                                                                                         cacCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCaa                     <  3:1583/149‑1 (MQ=255)
                                                                                                                         cacCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCaa                     <  3:54585/149‑1 (MQ=255)
                                                                                                                          acCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCa                                                            <  2:124472/109‑1 (MQ=255)
                                                                                                                           ccTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCCCGAAATCGCCTAGCccacca                   >  7:113057/1‑149 (MQ=255)
                                                                                                                             tGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAACTCGCCTAGCCCAACAcc                 >  8:330408/1‑149 (MQ=255)
                                                                                                                                  cTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCAACACCTcaa             <  3:287688/148‑1 (MQ=255)
                                                                                                                                      cccACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGgcgc                                      >  7:22864/1‑119 (MQ=255)
                                                                                                                                      cccACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGgcgc                                      <  8:22864/119‑1 (MQ=255)
                                                                                                                                           cTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCAACACCTCAAGCAAGCGGAg   <  3:97491/149‑1 (MQ=255)
                                                                                                                                            tCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCAACACCTCAAGCAAGCGGAGc  >  6:233333/1‑149 (MQ=255)
                                                                                                                                              agcagcAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGAc                                            <  3:283580/105‑1 (MQ=255)
                                                                                                                                              agcagcAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGAc                                            >  4:283580/1‑105 (MQ=255)
                                                                                                                                                    aaGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGc                                                 >  8:128877/1‑94 (MQ=255)
                                                                                                                                                    aaGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGc                                                 <  7:128877/94‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CCTGCTGGCACGTGGTCACTACCCGGCTTTTCTGAAGCAGCGGGTGGGTGGCATGCAAGGGCATTTGAACAACCACCAGGCTGCGCGCCTGGTGCACGAGTTGGGCTGGAACAACCTGCAACACCTGGTGCTGGCCCACCTCAGCAGCAAGAACAACCTGCCACACTTGGCCCGCCAGTGCTTTGTCGACACCCTCGGGTGTGACCCGGACTGGCTCCAGGTGGCGAATCAGGAACACGGGCTCGACTGGCGCGAAATCGCCTAGCCCAACACCTCAAGCAAGCGGAGC  >  NC_002947/1417247‑1417535

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