Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,053,627 T→C 10.3% intergenic (‑27/+39) PP_1828 ← / ← PP_1829 hypothetical protein/alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,053,6270TC10.3% 152.4 / 8.3 58intergenic (‑27/+39)PP_1828/PP_1829hypothetical protein/alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base C (4/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TCGGTCGGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCC  >  NC_002947/2053479‑2053771
                                                                                                                                                    |                                                                                                                                                 
gcgggcggCTGGGCTGCCAGGCGTAGCAGCTGGAAAGGTTCGGCGGGCAGCTCGGCCCCGGGCAGGGCGGGGGTTGCCTGGGGTTCTTCGGCCGGCTCGCCGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCctgct                                                                                                                                                   <  4:161703/148‑1 (MQ=255)
    gcggCAGTGCTGCCTGGCGTTGCAGCTGGAAATGTGCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTAt                                                                                                                                               <  3:46161/148‑1 (MQ=255)
       gCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCgg                                                                                                                                            >  1:289023/1‑149 (MQ=255)
        cAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCgg                                                                                                                                            >  1:31783/1‑148 (MQ=255)
          gTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGGGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCggcag                                                                                                                                         <  5:34598/149‑1 (MQ=255)
                ccAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATgg                                                                                                                                   >  1:302406/1‑149 (MQ=255)
                   ggCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCt                                                                                                                                 >  8:101456/1‑148 (MQ=255)
                    gCGTAGCAGCTGGAAATGTTCGGCGGTCTGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCCATGCCTGCCGATACAGCAGGCATAGCTga                                                                                                                               <  2:256098/149‑1 (MQ=38)
                     cGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTgag                                                                                                                              <  2:31783/149‑1 (MQ=255)
                      gTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTgaga                                                                                                                             >  5:12149/1‑149 (MQ=255)
                           agcTGTAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGa                                                                                                                        <  5:301/149‑1 (MQ=255)
                                      ttgggcggTCCGCGCGGCCCCGGGCAGGGCGGGGGTTTCCTGGGGTGCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGAc                                                                                                            <  8:76447/146‑1 (MQ=255)
                                      ttCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGAc                                                                                                            >  6:17345/1‑149 (MQ=255)
                                      ttCGGCGGTCAGCCCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGAc                                                                                                            >  7:355770/1‑149 (MQ=255)
                                         ggcggGCCGCGCGGCCCCGGGCTGGGCGGGGGTTTCCTGGGGTTCTGCGGCCGGCGCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTa                                                                                                         <  3:220385/149‑1 (MQ=255)
                                         ggcggGCAGCGCGGCCCAGGGCTGGGGGGGGGTTTCGTGGGGTTCTTCGGCCGGCGCGACGGCGGTGGGGGTGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTa                                                                                                         <  3:215392/149‑1 (MQ=255)
                                               cAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCg                                                                                                    >  7:311374/1‑148 (MQ=255)
                                                aGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGaa                                                                                                  <  2:340447/149‑1 (MQ=255)
                                                 gCTCGGCCCCGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGcc                                                                                                     <  6:200797/145‑1 (MQ=255)
                                                 gCTCGGCCCCGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGcc                                                                                                     >  5:200797/1‑145 (MQ=255)
                                                  cTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAgg                                                                                                >  2:383641/1‑149 (MQ=255)
                                                   tCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAgg                                                                                                >  5:103186/1‑148 (MQ=255)
                                                   tCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTAGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAgg                                                                                                >  5:386225/1‑148 (MQ=255)
                                                            ggcagggcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTg                                                                                       >  8:102708/1‑148 (MQ=255)
                                                             gcagggcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGa                                                                                     <  7:101456/149‑1 (MQ=255)
                                                             gcagggcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGa                                                                                     <  6:103186/149‑1 (MQ=255)
                                                              ctggggggggGTTTCCTGGGGTTCTTCGGCCGGCTCGAGGGCGGGGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGAt                                                                                    <  2:289023/147‑1 (MQ=255)
                                                                  gggggggTTTCCTGGTGTTCTTCGGCCGGATCTACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCATCTATGCATTCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAAc                                                                                <  2:132751/147‑1 (MQ=255)
                                                                   cGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACg                                                                               >  5:144230/1‑149 (MQ=255)
                                                                    gggggTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGt                                                                              >  3:200866/1‑149 (MQ=255)
                                                                    gggggTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGt                                                                              <  6:12149/149‑1 (MQ=255)
                                                                       ggTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTcac                                                                           <  5:17345/149‑1 (MQ=255)
                                                                             cTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATtcg                                                                     >  7:171558/1‑149 (MQ=255)
                                                                                gggTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATtcgctc                                                                  >  8:200991/1‑149 (MQ=255)
                                                                                   ttcttcGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTcac                                                                           >  1:256937/1‑137 (MQ=255)
                                                                                   ttcttcGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTcac                                                                           <  2:256937/137‑1 (MQ=255)
                                                                                     cgacggccggcGCGACGGCGGGGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCTCCCa                                                             <  2:302406/146‑1 (MQ=255)
                                                                                      ttcGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAg                                                            >  6:226415/1‑149 (MQ=255)
                                                                                           ccggcTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCAc                                                       <  8:86537/149‑1 (MQ=255)
                                                                                                     ggcggTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCgg                                             <  4:165645/149‑1 (MQ=255)
                                                                                                                   tggtCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGa                               <  6:144230/149‑1 (MQ=255)
                                                                                                                       catcaGGGGTTGTCCTCAGATATTCCTGCTGTATCTGCAGGCATTGCTGAGACCTGAGACTTTTTGACTTAGGCCGAAGTCAAATGCTGGATTAACGTCACAATTCGCTCGCCCTGCTCACGCGCCAGGGGAAGCTGGGGGTTGAGGt                            <  6:27342/148‑1 (MQ=255)
                                                                                                                        atcaGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAg                          <  8:171558/149‑1 (MQ=255)
                                                                                                                        atcaGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAg                          >  8:104394/1‑149 (MQ=255)
                                                                                                                            gggCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGt                      >  1:81289/1‑149 (MQ=255)
                                                                                                                             ggCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCGCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGt                      >  1:296663/1‑149 (MQ=255)
                                                                                                                                    ccTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctgc              >  6:297355/1‑149 (MQ=255)
                                                                                                                                     cTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctgct             >  7:109488/1‑149 (MQ=255)
                                                                                                                                     cTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctgc              <  8:311374/148‑1 (MQ=255)
                                                                                                                                      tCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCGGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCCGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGctgctg            >  1:297640/1‑149 (MQ=255)
                                                                                                                                      tCAGCCATGCCTGCCGATACAGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGctgctg            >  7:45366/1‑149 (MQ=255)
                                                                                                                                       cAGCCATGCCTGCCGATACAGCAGGCATAGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTcgc                   >  7:119595/1‑141 (MQ=38)
                                                                                                                                       cAGCCATGCCTGCCGATACAGCAGGCATAGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTcgc                   <  8:119595/141‑1 (MQ=38)
                                                                                                                                       cAGCCATGCCTGCCGATACAGCAGGCATAGCTGAGACCTGAGCCTGTTTGACTTACGCCGAACGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTgg           >  8:274832/1‑149 (MQ=38)
                                                                                                                                        aGCCATGCCTGCCGATACAGCAGGCATAGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGc          >  7:252730/1‑149 (MQ=255)
                                                                                                                                            ctGCCTGCTGTATCGCCAGGCATGGCTGAGACCTGAGCCGGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCAt      <  8:347027/148‑1 (MQ=255)
                                                                                                                                            aTGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGTACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAgg                             <  1:259304/126‑1 (MQ=255)
                                                                                                                                            aTGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAgg                             >  2:259304/1‑126 (MQ=255)
                                                                                                                                               ccTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGcc   >  6:22518/1‑149 (MQ=255)
                                                                                                                                                ctgctgTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAACTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGccc  >  1:312948/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCGGTCGGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCC  >  NC_002947/2053479‑2053771

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