Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,073,444 G→A 13.6% P271L (CCA→CTA) ‡ PP_5327 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,073,4440GA13.6% 164.7 / 14.2 66P271L (CCA→CTA) ‡PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (5/4);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

CGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAG  >  NC_002947/6073298‑6073583
                                                                                                                                                  |                                                                                                                                           
cGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGc                                                                                                                                           >  2:146565/1‑149 (MQ=255)
     aTGATCAGGGGCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAAc                                                                                                                                      <  6:279465/149‑1 (MQ=255)
        aTCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACgg                                                                                                                                    <  1:146565/148‑1 (MQ=255)
         tCAGGGTCATCATCACGGTACCGAAGGTGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCa                                                                                                                                  >  1:183574/1‑149 (MQ=255)
         tCAGGGTCATCATCACGGTACCGAAGGTGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCa                                                                                                                                  >  1:189448/1‑149 (MQ=255)
           aGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATg                                                                                                                                <  5:35866/149‑1 (MQ=255)
              gTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACg                                                                                                                                     >  4:288466/1‑141 (MQ=255)
              gTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACg                                                                                                                                     <  3:288466/141‑1 (MQ=255)
                catcatcaCGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGTCCAACTTGGTGAACGCCGACATGCCGTTCAGCCAAACGGCATgtcggc                                                                                                                           <  7:282269/149‑6 (MQ=255)
                        cGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCTCGCGGGTCGTCTCTGAGGAACTCCCAGACCTTGCCGAAGTAGGGGCCCAACTTGGTGAACGACGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTc                                                                                                                   <  8:85272/149‑1 (MQ=255)
                             ccGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGt                                                                                                               <  3:234962/148‑1 (MQ=255)
                              cGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATg                                                                                                                                >  7:168706/1‑130 (MQ=255)
                              cGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATg                                                                                                                                <  8:168706/130‑1 (MQ=255)
                                     ggCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAgg                                                                                                                >  2:168354/1‑139 (MQ=255)
                                     ggCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAgg                                                                                                                <  1:168354/139‑1 (MQ=255)
                                          ggAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTcn                                                                                                 <  2:189448/149‑2 (MQ=255)
                                          ggAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCa                                                                                                 <  2:183574/149‑1 (MQ=255)
                                                tACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCg                                                                                           >  1:311465/1‑149 (MQ=255)
                                                   cGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTg                                                                                        <  6:314360/149‑1 (MQ=255)
                                                    gCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGc                                                                                       >  4:272312/1‑149 (MQ=255)
                                                       cTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCAt                                                                                    <  2:311465/149‑1 (MQ=255)
                                                            tgtTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                               <  6:241331/149‑1 (MQ=255)
                                                            tgtTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                               <  2:51190/149‑1 (MQ=255)
                                                                ggCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACgg                                                                                              <  8:159350/130‑1 (MQ=255)
                                                                ggCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACgg                                                                                              >  7:159350/1‑130 (MQ=255)
                                                                ggCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGcca                                                                           >  6:8134/1‑149 (MQ=255)
                                                                   ttCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGt                                                                                                               >  1:78746/1‑110 (MQ=255)
                                                                   ttCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGt                                                                                                               <  2:78746/110‑1 (MQ=255)
                                                                     cACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGAt                                                                                                           >  4:208906/1‑112 (MQ=255)
                                                                     cACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGAt                                                                                                           <  3:208906/112‑1 (MQ=255)
                                                                      aCGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGc                                                                     >  2:95043/1‑149 (MQ=255)
                                                                                tcgCTGAGGAACTCCCAGACATTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGAATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCAGTTGCCATACCCGGCCACCAGGCTGTAGCggg                                                            >  4:293329/1‑148 (MQ=255)
                                                                                   cTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACt                                                                                                                            <  7:259136/81‑1 (MQ=255)
                                                                                   cTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACt                                                                                                                            >  8:259136/1‑81 (MQ=255)
                                                                                           cTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                >  8:274481/1‑149 (MQ=255)
                                                                                           cTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                >  2:34107/1‑149 (MQ=255)
                                                                                           cTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                <  7:274481/149‑1 (MQ=255)
                                                                                            tCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTg                                               >  6:246645/1‑149 (MQ=255)
                                                                                              ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGc                                              <  5:246645/148‑1 (MQ=255)
                                                                                                 gACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTg                                                  >  3:142690/1‑141 (MQ=255)
                                                                                                 gACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTg                                                  <  4:142690/141‑1 (MQ=255)
                                                                                                 gACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGc                                          <  1:95043/149‑1 (MQ=255)
                                                                                                   ccTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCg                                        <  5:8134/149‑1 (MQ=255)
                                                                                                    cTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGATAAAGCCTGGTGCAGGCCgg                                       >  7:228089/1‑149 (MQ=255)
                                                                                                         cGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTg                                                  <  1:161911/133‑1 (MQ=255)
                                                                                                         cGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTg                                                  >  2:161911/1‑133 (MQ=255)
                                                                                                           aaGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTAACTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCg                                <  8:117058/149‑1 (MQ=255)
                                                                                                                        aCTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGcc                   <  1:34107/149‑1 (MQ=255)
                                                                                                                        aCTTGGTGAACGCCGACATGCCGTTTGGCTCAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGcc                   >  1:31235/1‑149 (MQ=255)
                                                                                                                             gTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTg                >  7:162601/1‑147 (MQ=255)
                                                                                                                                    ccGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCa                                             >  5:70722/1‑111 (MQ=255)
                                                                                                                                    ccGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCa                                             <  6:70722/111‑1 (MQ=255)
                                                                                                                                     gtgcATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCgg                                                             <  4:9866/91‑1 (MQ=38)
                                                                                                                                     gtgcATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCgg                                                             >  3:9866/4‑94 (MQ=38)
                                                                                                                                     gtgcATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGt      >  8:31343/4‑149 (MQ=255)
                                                                                                                                      gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCg           >  4:11687/1‑143 (MQ=255)
                                                                                                                                      gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCg           <  3:11687/143‑1 (MQ=255)
                                                                                                                                      gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCg           <  3:6634/143‑1 (MQ=255)
                                                                                                                                      gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCg           >  4:6634/1‑143 (MQ=255)
                                                                                                                                      gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGtt     <  3:272312/149‑1 (MQ=255)
                                                                                                                                       aCATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGtt     >  4:275099/1‑148 (MQ=255)
                                                                                                                                        cATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGt                                                                  >  6:20444/1‑86 (MQ=38)
                                                                                                                                        cATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGt                                                                  <  5:15272/86‑1 (MQ=38)
                                                                                                                                        cATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGt                                                                  <  5:20444/86‑1 (MQ=38)
                                                                                                                                        cATGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGt                                                                  >  6:15272/1‑86 (MQ=38)
                                                                                                                                         aTGCCGTTCAGCCAAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAg  >  3:262028/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
CGGCCATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAG  >  NC_002947/6073298‑6073583

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: