Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,073,447 T→C 13.8% Q270R (CAG→CGG) ‡ 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,4470TC13.8% 153.3 / 16.2 65Q270R (CAG→CGG) ‡PP_5327phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base C (5/4);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

ATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTC  >  NC_002947/6073303‑6073595
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                                 <  3:288466/141‑1 (MQ=255)
         gTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACg                                                                                                                                                 >  4:288466/1‑141 (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                                                                                           <  2:51190/149‑1 (MQ=255)
                                                       tgtTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCgg                                                                                           <  6:241331/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                                                            <  7:274481/149‑1 (MQ=255)
                                                                                      cTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                            >  2:34107/1‑149 (MQ=255)
                                                                                      cTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGt                                                            >  8:274481/1‑149 (MQ=255)
                                                                                       tCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTg                                                           >  6:246645/1‑149 (MQ=255)
                                                                                         ccAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGc                                                          <  5:246645/148‑1 (MQ=255)
                                                                                            gACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTg                                                              <  4:142690/141‑1 (MQ=255)
                                                                                            gACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTg                                                              >  3:142690/1‑141 (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                       <  3:11687/143‑1 (MQ=255)
                                                                                                                                 gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCg                       <  3:6634/143‑1 (MQ=255)
                                                                                                                                 gACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCg                       >  4:11687/1‑143 (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)
                                                                                                                                                tGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTc  >  1:260956/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATGATCAGGGTCATCATCACGGTACCGAAGATGGCCGGGAAGATACCGCCCTCGGTGTTGGCTTCACGCGGGTCGTCGCTGAGGAACTCCCAGACCTTGCCGAAGTAGGTGCCCAACTTGGTGAACGCCGACATGCCGTTTGGCTGAACGGCATGCACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTC  >  NC_002947/6073303‑6073595

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: