Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,989,441 T→C 20.6% intergenic (+32/+57) PP_5242 → / ← ycaC‑II GAF domain‑containing sensor histidine kinase/putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,989,4410TC20.6% 362.6 / 90.9 205intergenic (+32/+57)PP_5242/ycaC‑IIGAF domain‑containing sensor histidine kinase/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (56/91);  new base C (22/16);  total (95/110)
Fisher's exact test for biased strand distribution p-value = 4.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-04

Only 100 of 206 total aligned reads displayed.

GATGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCG  >  NC_002947/5989292‑5989588
                                                                                                                                                     |                                                                                                                                                     
 aTGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                                      <  1:482916/150‑1 (MQ=255)
 aTGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                                      <  1:5050279/150‑1 (MQ=255)
  tGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTcc                                                                                                                                                     <  1:1695290/150‑1 (MQ=255)
     gCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                                  <  2:825536/150‑1 (MQ=255)
     gCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCGCCTAc                                                                                                                                                  >  1:413206/1‑150 (MQ=255)
       ttGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGGCGGCCTCTTCGCGGGTAAACCCGCTCCTACgg                                                                                                                                                >  2:5087404/1‑150 (MQ=255)
         ggCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTc                                                                                                                                              <  1:5732257/150‑1 (MQ=255)
          gCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCg                                                                                                                                             <  2:6811180/150‑1 (MQ=255)
           cGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCgg                                                                                                                                            <  2:2691516/150‑1 (MQ=255)
           cGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCgg                                                                                                                                            <  1:5466470/150‑1 (MQ=255)
            gATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCACCAACGGTCGGc                                                                                                                                           >  2:2885054/1‑150 (MQ=255)
              tCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGGCTGTGCCGGCCGCTGCGCGGGGACCCCCGCTGCGACggggggcgg                                                                                                                                         >  1:3484091/1‑143 (MQ=255)
                tGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTAGCGGGTGAACCCCCGCCAACCggcggcgggt                                                                                                                                       >  1:151775/1‑146 (MQ=255)
                     ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGc                                                                                                                                  <  1:2921541/150‑1 (MQ=255)
                     ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGc                                                                                                                                  <  1:4943809/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  1:5980912/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  1:5223759/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  2:4088778/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  2:122586/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  2:3806415/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  1:331302/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  1:3590178/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                               <  1:2692610/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGCCGACCGTAGGAGCgg                                                                                                                               <  2:856756/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                               <  2:5932395/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                               <  2:2207025/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                               <  1:2597292/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                               <  2:4343571/150‑1 (MQ=38)
                        ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                              <  2:4053718/150‑1 (MQ=255)
                        ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                              <  1:2806140/150‑1 (MQ=255)
                         gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                             <  2:1368652/150‑1 (MQ=255)
                         gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGt                                                                                                                             <  2:4942039/150‑1 (MQ=38)
                           ggCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGttt                                                                                                                           <  1:1398269/150‑1 (MQ=37)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTaac                                                                                                                        >  2:5504369/1‑148 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCGGCTGGAGGGGCGGGTTTAcc                                                                                                                        >  1:2223064/1‑150 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTACTACGGTCGGCGGTAGGGGCGGGgtttgc                                                                                                                        >  2:1704412/1‑147 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCGCCTACGGTCGGCTGTAGGGGGGGGTTTaac                                                                                                                        >  1:6640609/1‑148 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGAGTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAcc                                                                                                                        >  2:5625661/1‑150 (MQ=255)
                               cAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAccc                                                                                                                       <  1:2474308/150‑1 (MQ=255)
                               cAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAccc                                                                                                                       <  1:999301/150‑1 (MQ=255)
                                 ggCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCcgc                                                                                                                     <  1:5376542/150‑1 (MQ=255)
                                 ggCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCACCCACCGGCGGGGGTTGGGGGGGGGTTAccccc                                                                                                                     >  1:3156986/1‑148 (MQ=255)
                                  gCCCACGGCGGGCCGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCACGTGGCGGGGGAACCCCCACCAACGGGCGGGGGTTGGGGCGGGGTTTCCcgcc                                                                                                                    >  1:5161076/1‑149 (MQ=255)
                                  gCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCcgcg                                                                                                                    <  1:6243587/150‑1 (MQ=255)
                                  gCCCACGGCGGGCAGCTGCAGGCGCGGGTGTCGCAGTTGGGGGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATCGCCTGTGCCGGCCTCGTCGCGGGTAAACCCGCACCTACGGTCGGCTGTAGGGGCGGGTTTACCcgcg                                                                                                                    <  2:4750734/150‑1 (MQ=255)
                                      ccggcggGCAGCGGCCGGCGCGGGTGTCGCAGTTGGGTGGATGCGGGATTCGTGTCAGTGTGCCTTGTACTCACGCAGAGCCTGTGCCGGCCGCTTCGCGGGTAAACCGGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAaga                                                                                                                <  2:6404964/149‑1 (MQ=37)
                                      aCGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAaga                                                                                                                <  2:5592551/150‑1 (MQ=255)
                                      aCGGCGGGCAGCGGCCGGCGCTGGTGTCGCAGTTGGGTGTATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTACGTTCGGCTGTAGGGGCGGGTTTACCCGCGAaga                                                                                                                <  2:5315406/150‑1 (MQ=39)
                                             gCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGt                                                                                                         >  2:4879512/1‑150 (MQ=255)
                                             gCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                               >  2:6855157/1‑128 (MQ=37)
                                              cAGCTGCCGGCGCTGGTGTCGCTGTTGGGTGGTTGCCGGATTCGTGTCAGTCTGCCTTTGAAACCCGCATAGCCTGTGCCGGCCTCTTCGCGGTTAAACCCGCTCCTCCGGTCGGCTGTCGGGGGGGGTTTACCCGCGAAGAGGCCAGTa                                                                                                        <  2:3585675/150‑1 (MQ=38)
                                                  tGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCACCAACGGTGGGCGGTAGGGGCGGGTTTAACCCCGAAGTGGCCAGAACAgg                                                                                                    >  1:1519322/1‑150 (MQ=255)
                                                     aGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCACTACAGGttt                                                                                                 >  2:2062660/1‑150 (MQ=255)
                                                         gcTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCACCTACGGTCGGCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACACGTTTACaa                                                                                             >  2:1312983/1‑149 (MQ=38)
                                                          cTGGTGTCGCTGTTGGGTGGTTGCCGGTTTCGTGTCAGTCTGCCTTTGAACCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTCGGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATc                                                                                            <  2:315412/150‑1 (MQ=38)
                                                          cTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCggg                                                                                                                              >  2:5004281/1‑116 (MQ=25)
                                                           tGGTGTCGCAGTTGGGGGGATGCCGGATTCGTGTCAGTCTGCCTTGGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCa                                                                                           <  2:2667846/150‑1 (MQ=255)
                                                           tGGTGTCGCAGTGGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCa                                                                                           <  2:1448825/150‑1 (MQ=255)
                                                                 cGCTGTTGGTTGGATGCCGTATTCGTGTCAGTCTGCCTTTTAAACACGCATAGCCTGGGCCGGCCTCTTCGCGGTTAACCCCGCTCCTCGGGTCGGCTGTTGGGGGGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGa                                                                                     <  2:1139472/150‑1 (MQ=255)
                                                                    aGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGc                                                                                  <  2:6171651/150‑1 (MQ=255)
                                                                        gggTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAgc                                                                              >  2:3346216/1‑150 (MQ=255)
                                                                         ggTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGGAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAgcg                                                                             >  2:6796759/1‑150 (MQ=255)
                                                                           tGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAACCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCt                                                                           >  2:5335904/1‑150 (MQ=255)
                                                                            ggATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTg                                                                          >  2:2219937/1‑150 (MQ=255)
                                                                             gATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCATCGCGGGTAAAACCCCACCTACCGTCCGCGGTAGGGGCGGGGTTAACCGCGAAGAGGCCACTACAGGTTTACAAAAGCCCGAGGACACCGCTGc                                                                         >  2:2336992/1‑150 (MQ=255)
                                                                                gCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTCGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTa                                                                      <  1:1004090/150‑1 (MQ=255)
                                                                                gCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCACCTACGGGGGGGGGGGGGGGCGGGTTTTCCCGCGAAGGGGCCCGGACAGGTTTACAACACCGCGCGGCCCGCGCTGCgga                                                                      >  1:1101215/1‑148 (MQ=25)
                                                                                          tgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCCCGAAGAGGCCACTACAGGTTTACATAAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatc                                                            >  2:3207296/1‑150 (MQ=255)
                                                                                          tgtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCAAAGGTCGGCTGTAGGGGCGGGTTTACCCGCGCAGAGGCCAGTACAGGGTTACATCAGCGCCAGGCCAGCGCTGCGTAGCTGGtcctc                                                            >  2:1835001/1‑150 (MQ=37)
                                                                                             cAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAgg                                                        <  2:2002080/150‑1 (MQ=255)
                                                                                              aGTCTGCGTTTGAATCACGCTTAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCGAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGt                                                       <  2:6441775/150‑1 (MQ=255)
                                                                                                 cTGCCTTTGAATCCCGCATCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGc                                                    <  2:4246158/150‑1 (MQ=255)
                                                                                                   gCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACGGTCGGCGGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCAGCGTAGCTGTTCATCAGGTTGCgg                                                  >  2:2164103/1‑150 (MQ=255)
                                                                                                         gAATCCCGCATCGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGGCGGCTGTTGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTg                                            <  2:4590994/150‑1 (MQ=38)
                                                                                                          aaTCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                           <  1:5087404/150‑1 (MQ=255)
                                                                                                           aTCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                          <  2:4828661/150‑1 (MQ=255)
                                                                                                           aTCACGCATAGCCTGTGCCGGCCTCTTAGCGGGTAAACCCCCTCCAACGGTCGGCGGTAGGGGCGGGGTTTACCCCGCAGAGGCCAGTACAGGTTTTAAACAGCCCGCGGCCAGCGCGGCGTAGCTGTTCATCAGGTTGCGGGGGTTggg                                          >  1:4992246/1‑150 (MQ=255)
                                                                                                            tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAACCCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgggg                                         <  1:5607518/150‑1 (MQ=255)
                                                                                                            tCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAAACGCTCCAACGGTCGGGGGTAGGGGGGGGTTTAACCGCGAAGAGGCCACTAAAGTTTTAAATCATCCCGAGGCCACAGCAGCGTAGAAGTTCATCATGGTTGGGGAGTTgggg                                         >  2:3686537/1‑150 (MQ=255)
                                                                                                                         gtgCCGGCCTCTTCGCGGGTAAACCCGCTCCTGCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAca                            <  2:1548226/150‑1 (MQ=255)
                                                                                                                                 ctctTCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTa                    <  2:413206/150‑1 (MQ=255)
                                                                                                                                     tCGCGGGTCACCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATGCGCTGCGACAGCAGGTTACCCa                <  1:3908507/150‑1 (MQ=255)
                                                                                                                                        cgGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGcc             >  2:5720787/1‑150 (MQ=255)
                                                                                                                                        cgGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGcc             >  2:2498286/1‑150 (MQ=255)
                                                                                                                                            tAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTc         >  2:2630761/1‑150 (MQ=37)
                                                                                                                                              aaCCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGa       >  1:79612/1‑150 (MQ=38)
                                                                                                                                               aCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt      >  2:6409841/1‑150 (MQ=255)
                                                                                                                                                cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg     >  1:3480704/1‑150 (MQ=255)
                                                                                                                                                cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg     >  1:4461283/1‑150 (MQ=38)
                                                                                                                                                 ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt    >  1:6417350/1‑150 (MQ=255)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt    >  2:5354415/1‑150 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt    >  1:4258082/1‑150 (MQ=37)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt    >  2:6081705/1‑150 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt    >  2:253850/1‑150 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt    >  1:1927644/1‑150 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt    >  2:696543/1‑150 (MQ=38)
                                                                                                                                                  cGCTCCTCGGGTGGGCTGTTGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTc   <  1:1237812/150‑1 (MQ=39)
                                                                                                                                                   gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCg  >  1:1662025/1‑150 (MQ=255)
                                                                                                                                                   gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGGTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTGGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCg  >  2:1108597/1‑150 (MQ=255)
                                                                                                                                                   gCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGGTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGGTCATCAGGTGGCGGGAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCg  >  2:2080097/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCG  >  NC_002947/5989292‑5989588

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: