Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,053,631 T→A 9.4% intergenic (‑31/+35) 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,6310TA 9.4% 174.0 / 8.3 64intergenic (‑31/+35)PP_1828/PP_1829hypothetical protein/alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (4/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

TCGGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGGGG  >  NC_002947/2053483‑2053777
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                           >  5:200797/1‑145 (MQ=255)
                                             gCTCGGCCCCGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGcc                                                                                                           <  6:200797/145‑1 (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                                                                                           <  6:103186/149‑1 (MQ=255)
                                                         gcagggcGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGa                                                                                           <  7:101456/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                                                                                 <  2:256937/137‑1 (MQ=255)
                                                                               ttcttcGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTcac                                                                                 >  1:256937/1‑137 (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:104394/1‑149 (MQ=255)
                                                                                                                    atcaGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAg                                <  8:171558/149‑1 (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)
                                                                                                                                                 gTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTggg   >  6:213579/1‑149 (MQ=255)
                                                                                                                                                  tATCGGGAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATCAACGTCACAATTCGCTCCCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTgg    <  6:386225/147‑1 (MQ=255)
                                                                                                                                                  tATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTgggg  >  7:32675/1‑149 (MQ=255)
                                                                                                                                                  tATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTgggg  <  2:296663/149‑1 (MQ=255)
                                                                                                                                                  tATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTgg    <  7:200991/147‑1 (MQ=255)
                                                                                                                                                  tATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTgg    <  7:102708/147‑1 (MQ=255)
                                                                                                                                                    tCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctgc                    <  1:61369/129‑1 (MQ=255)
                                                                                                                                                    tCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCgctgc                    >  2:61369/1‑129 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCGGCAGTGCTGCCAGGCGTAGCAGCTGGAAATGTTCGGCGGTCAGCTCGGCCCAGGGCAGGGCGGGGGTTTCCTGGGGTTCTTCGGCCGGCTCGACGGCGGTGGGGGAGGTGGTCATCAGGGCTTGTCCTCAGCTATGCCTGCTGTATCGGCAGGCATGGCTGAGACCTGAGCCTGTTTGACTTACGCCGAAGGCAAATGCTGGATGAACGTCACAATTCGCTCGCCCAGCTCACGCGCCAGCGGCAGCTCGGGGTTGAGGTAGCTGTCGCGCTGCTGGCTCATGCCCTTGGGG  >  NC_002947/2053483‑2053777

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: