Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,618,414 C→A 5.9% S162I (AGC→ATC)  PP_1418 ← tricarboxylate transport protein TctC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,618,4140CA 5.9% 217.1 / 3.5 68S162I (AGC→ATC) PP_1418tricarboxylate transport protein TctC
Reads supporting (aligned to +/- strand):  ref base C (30/34);  new base A (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.45e-01

GGAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCACCGC  >  NC_002947/1618269‑1618562
                                                                                                                                                  |                                                                                                                                                    
ggAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGc                                                                                                                                                   >  5:231717/1‑149 (MQ=255)
ggAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCAAGTAGCGCAGTTCGCGCGGGATGATGCCAGCGGACTTGGCGATCAACGCAGTCTGCTTCCAGTCATGGCTGc                                                                                                                                                   >  8:273325/1‑149 (MQ=255)
  aaaTGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGcc                                                                                                                                                  >  4:312103/1‑148 (MQ=255)
     tGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCAGCCAGGACTGGATGCAGAct                                                                                                                                              <  2:138453/149‑2 (MQ=255)
        cGGTTGAGCCGACCTGGATGTAGCCGCCCAGCAATGGCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCATGTCGCGCGGGTTGCATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACgg                                                                                                                                             <  1:286555/148‑1 (MQ=255)
                ccGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGa                                                                                                                                   <  1:208988/149‑1 (MQ=255)
                 cGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGTCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCgg                                                                                                                                    <  4:10531/147‑1 (MQ=255)
                   aCCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAAcc                                                                                                                                >  1:338282/1‑149 (MQ=255)
                     cTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCg                                                                                                                               >  2:170971/1‑148 (MQ=255)
                        gATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGAt                                                                                                                             <  1:62348/147‑1 (MQ=255)
                        gATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATc                                                                                                                            <  2:303549/148‑1 (MQ=255)
                        gATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCa                                                                                                                           >  7:54892/1‑149 (MQ=255)
                           gtgGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGa                                                                                                                        >  4:363720/1‑149 (MQ=255)
                             ggCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGAcc                                                                                                                      >  2:88950/1‑149 (MQ=255)
                                cgccCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTg                                                                                                                   <  1:42126/149‑1 (MQ=255)
                                   ccAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTc                                                                                                                <  4:70135/149‑1 (MQ=255)
                                    cAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGTCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACTGTACCGCCGGAACCGATCACGACCTTGCTCg                                                                                                               <  7:44467/149‑1 (MQ=255)
                                    cAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCg                                                                                                               >  4:252347/1‑149 (MQ=255)
                                    cAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCg                                                                                                               >  4:20321/1‑149 (MQ=255)
                                        aaTGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTAc                                                                                                                                          >  3:55902/1‑118 (MQ=255)
                                        aaTGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTAc                                                                                                                                          <  4:55902/118‑1 (MQ=255)
                                           gCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGtcttt                                                                                                        >  4:79142/1‑149 (MQ=255)
                                               tGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTttcttc                                                                                                    <  3:312103/149‑1 (MQ=255)
                                               tGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTttcttc                                                                                                    <  3:20321/149‑1 (MQ=255)
                                                 gCGATTTCGCCGCCGCCCTCCAGGGCAACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAg                                                                                                  >  1:5688/1‑149 (MQ=255)
                                                      ttCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTg                                                                                             >  1:364279/1‑149 (MQ=255)
                                                      ttCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTg                                                                                             >  1:201934/1‑149 (MQ=255)
                                                         gccgccgccCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGa                                                                                          <  5:41971/149‑1 (MQ=255)
                                                          ccgccgccCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGAc                                                                                         >  1:14443/1‑149 (MQ=255)
                                                             ccgccCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAg                                                                                      <  3:363720/149‑1 (MQ=255)
                                                                cccTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTGTCTTCAGCGCTGCGACCAGGt                                                                                    >  2:355789/1‑148 (MQ=255)
                                                                 ccTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGTCCTTGGAGCTAAACGCAGTCTGCATCCAGTCCTGGCTGCAGACGGTAACGCCGGAACCGATAAAGACCTTGCTACGGTCTTTCTTCAGAGCTGCTACCACGtc                                                                                   >  1:329144/1‑148 (MQ=255)
                                                                  cTCCAGGGCCACGTAGGGAAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGGCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGtcat                                                                                 <  7:273325/149‑1 (MQ=255)
                                                                  cTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGtcat                                                                                 <  2:201934/149‑1 (MQ=255)
                                                                   tCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAg                                                                                      >  5:68161/1‑143 (MQ=255)
                                                                   tCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAg                                                                                      <  6:68161/143‑1 (MQ=255)
                                                                     cAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGtcatc                                                                                <  3:252347/147‑1 (MQ=255)
                                                                     cAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCaa                                                                              <  1:170971/149‑1 (MQ=255)
                                                                       gggCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                            <  3:79142/149‑1 (MQ=255)
                                                                       gggCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                            <  1:88950/149‑1 (MQ=255)
                                                                       gggCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                            <  8:36037/149‑1 (MQ=255)
                                                                               tAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTa                                                                    <  1:355789/149‑1 (MQ=255)
                                                                                 gcgcAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGt                                                                     >  3:100567/1‑146 (MQ=255)
                                                                                 gcgcAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGt                                                                     <  4:100567/146‑1 (MQ=255)
                                                                                  cgcAGGTCGCGGGGGTTGATGCCTGCGGCCTTGGGGATCAAGGATTTGTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGTACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCATGGGTCTTGTAggg                                                                 <  8:119400/149‑1 (MQ=255)
                                                                                  cgcAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAggg                                                                 <  2:14443/149‑1 (MQ=255)
                                                                                                  tttaTGCCAGCGGCCTTGGCTATAAACGCAGTCTTAATCCAGTCCTGGCTGCTGACGTTACCGCTGGAACCGATCATGACCTTGCTCGGGTCTTTCTTCAGCGCTTCGACCAGGTCATCAATGGTCTTGTAGGGCGAGTCGCTCTTCAc                                                  <  2:329144/146‑1 (MQ=255)
                                                                                                   gATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACt                                                 <  8:54892/148‑1 (MQ=255)
                                                                                                         cAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATg                                           >  3:96805/1‑149 (MQ=255)
                                                                                                             ggCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCAc                                       <  6:334908/149‑1 (MQ=255)
                                                                                                                      gATCAACTCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTgg                              >  1:104089/1‑149 (MQ=255)
                                                                                                                      gATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTgg                              >  1:142935/1‑149 (MQ=255)
                                                                                                                           aCGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCa                                        <  5:321681/134‑1 (MQ=255)
                                                                                                                           aCGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCa                                        >  6:321681/1‑134 (MQ=255)
                                                                                                                            cGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGACGGAACCGATCACGACCTTGCTAGGGTCTTTCTTCAGCGCTGAGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACAGTAGCTGGTGCCGa                        >  1:243793/1‑149 (MQ=255)
                                                                                                                             gCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCAcc                                      >  2:52649/1‑134 (MQ=255)
                                                                                                                             gCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCAcc                                      <  1:52649/134‑1 (MQ=255)
                                                                                                                               aGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACgg                     <  2:104089/149‑1 (MQ=255)
                                                                                                                                gTCGGCAGCCAGGACTGGATGCAGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACgg                     >  1:63187/1‑148 (MQ=255)
                                                                                                                                 tcggCAGCCAGGACTGGATGCAGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                            >  7:32136/4‑92 (MQ=38)
                                                                                                                                 tcggCAGCCAGGACTGGATGCAGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAgg                                                                            <  8:32136/89‑1 (MQ=38)
                                                                                                                                 tCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGca                   <  2:364279/149‑1 (MQ=255)
                                                                                                                                 tCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGAGCCCGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCGTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGca                   <  6:231717/149‑1 (MQ=255)
                                                                                                                                   tGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGAACCGTAGCTGGTGCCGACGGcagc                 >  1:109130/1‑149 (MQ=255)
                                                                                                                                       tCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGcagccagc             >  7:93815/1‑149 (MQ=255)
                                                                                                                                          aGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCAt          <  2:5688/149‑1 (MQ=255)
                                                                                                                                            tCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCAtt         >  8:322580/1‑148 (MQ=255)
                                                                                                                                              cTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCa      >  2:50960/1‑149 (MQ=255)
                                                                                                                                                  cTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCACcgc  <  1:50960/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGAAATGTCGGTAGAGCCGACCTGGATGTGGCCGCCCAGCAATGCCGTGGCGATTTCGCCGCCGCCCTCCAGGGCCACGTAGCGCAGGTCGCGCGGGTTGATGCCAGCGGCCTTGGCGATCAACGCAGTCTGCATCCAGTCCTGGCTGCCGACGGTACCGCCGGAACCGATCACGACCTTGCTCGGGTCTTTCTTCAGCGCTGCGACCAGGTCATCAAGGGTCTTGTAGGGCGAGTCGCTCTTCACTGCGATGGCACCGTAGCTGGTGCCGACGGCAGCCAGCCATTTCACCGC  >  NC_002947/1618269‑1618562

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