Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,989,441 T→C 28.2% 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,4410TC28.2% 295.9 / 148.2 211intergenic (+32/+57)PP_5242/ycaC‑IIGAF domain‑containing sensor histidine kinase/putative hydrolase
Reads supporting (aligned to +/- strand):  ref base T (58/84);  new base C (23/34);  total (91/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-05

Only 100 of 220 total aligned reads displayed.

GATGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTT  >  NC_002947/5989292‑5989590
                                                                                                                                                     |                                                                                                                                                     
gATGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                                       <  2:5425469/150‑1 (MQ=255)
 aTGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTc                                                                                                                                                      <  2:5325980/150‑1 (MQ=255)
     gCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                                  <  2:6250544/150‑1 (MQ=255)
       ttGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                                       <  2:6965920/143‑1 (MQ=255)
             aTCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCt                                                                                                                                          <  2:6881511/150‑1 (MQ=255)
             aTCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCt                                                                                                                                          <  2:4155487/150‑1 (MQ=255)
              tCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCACTTCGCGGGTAAAACCCCGCCAACCGTCGGcgg                                                                                                                                         >  1:6930798/1‑148 (MQ=255)
                tGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTa                                                                                                                                       <  2:146229/150‑1 (MQ=255)
                 gTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAg                                                                                                                                      >  1:6945566/1‑150 (MQ=39)
                 gTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAg                                                                                                                                      <  2:6945566/150‑1 (MQ=39)
                  tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAgg                                                                                                                                     <  2:4982863/150‑1 (MQ=255)
                  tAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGGAAACCCGCACCAACCGTCGGCGGTAgg                                                                                                                                     >  1:5678916/1‑150 (MQ=255)
                   aGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGCCGACCGTAgga                                                                                                                                    <  1:5080820/150‑2 (MQ=255)
                     ttCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGc                                                                                                                                  <  1:935778/150‑1 (MQ=38)
                      tCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTACTGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCg                                                                                                                                 <  2:3119665/150‑1 (MQ=37)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                <  2:6154199/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                <  1:1271858/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                <  2:5739155/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCgg                                                                                                                                <  1:80129/150‑1 (MQ=255)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                <  2:3877734/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                <  1:888370/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                <  1:331294/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                <  2:1565541/150‑1 (MQ=38)
                       cGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCgg                                                                                                                                <  2:2939618/150‑1 (MQ=38)
                        ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                               <  2:2434826/150‑1 (MQ=255)
                        ggTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCggg                                                                                                                               <  1:3008365/150‑1 (MQ=255)
                         gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                              <  1:4151947/150‑1 (MQ=255)
                         gTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGt                                                                                                                              <  1:6410843/150‑1 (MQ=255)
                          tGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGtt                                                                                                                             <  2:1656644/150‑1 (MQ=255)
                             cTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAc                                                                                                                          <  2:3581707/150‑1 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAcc                                                                                                                         >  2:7414671/1‑150 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCGACGGTGGGCGGTAGGGGCGGGTTTAcc                                                                                                                         >  1:5891952/1‑150 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCCCGCCTACGGTCGGCGGTAGGGGCGGGTTTaac                                                                                                                         >  1:3636079/1‑148 (MQ=255)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTACTGGCCTCTTCGCGGGTAAACCCCCCCCCACAGCCGACCGGAGGAGCGGGTTTAcc                                                                                                                         >  1:2177865/1‑150 (MQ=25)
                              tCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGACTGCCATTGAAACACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGATTAcc                                                                                                                         <  2:4243444/150‑1 (MQ=255)
                               cAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTAccc                                                                                                                        <  2:6025299/150‑1 (MQ=255)
                                   cccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCGGTAGGGGCGGGTTTACCCGCGa                                                                                                                    >  2:795423/1‑150 (MQ=255)
                                   cccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTAGCGGGTAAACCCCCTCCTACCGTAGGCAGTATGGGCGGGTTTAACCGCGa                                                                                                                    >  1:6160309/1‑150 (MQ=255)
                                    ccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCcgc                                                                                                                      >  2:4130557/1‑147 (MQ=255)
                                    ccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACGGTCGGCGGTAGGGGCGGGGTTACCcgcgca                                                                                                                   >  1:603708/1‑148 (MQ=255)
                                    ccACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTACTGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGc                                                                                                                                  >  2:1367899/1‑135 (MQ=25)
                                      aCGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCACCAACGGTCGGCTGTAGGGGCGGGTTTTCCCGCGAaga                                                                                                                 >  1:1804248/1‑150 (MQ=255)
                                         gcggGCTGCTGCCGGCGGTGGTGTCGCTGTTGGGTGGATGCCGGATTGGGGTCCGTCTGCCTTTTAATCACGCCTAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTGGGCTGTAGGGGGGGTTTTACCCGCGAAGAGGc                                                                                                              <  2:3865807/150‑1 (MQ=37)
                                           gggCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCa                                                                                                            >  1:7326323/1‑150 (MQ=255)
                                              cAGCTGCAGGCGCTGGTGTCGCTGTTGGGTGGTTGCCGTATTCGTGTCAGTGTGCCTTTTAATCACTCATCGCGTGTGCCGGCCTCTGCGGGGGAAAACCCGCTCCTCGGTTGGGGTGTTGGGGGGGTTTTACCCGCGAAGAGGCCAGTa                                                                                                         <  1:5459233/150‑1 (MQ=255)
                                                 cggccggCGCTGGTGTCGCTGTTGGGTGGCTGCGGGTTTCGTTTCAGTCTGCCTTTTAATCCCGCATAGCCTGTGCCGGCCACTGCGGGGGTAAACCCGCTCCTCCGGTCGGCTGTCGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAg                                                                                                      <  2:2076422/145‑1 (MQ=255)
                                                 cTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCGACGGTCGGCTGTAGGGGCGGGGTTACCCGCGACGAGGCCAGTACAg                                                                                                      >  2:6463619/1‑150 (MQ=255)
                                                  tGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCACCGACGGTCGGCTGGAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAgg                                                                                                     >  1:2897538/1‑150 (MQ=255)
                                                         gcgggTGGCGCCGTTGGGTGGAGGCCGGATGCGTGTCCGTCTGCCTTTGAACCCCGCATCGCGTGTGCCGGCCGCTGGGCGGGTAAACCCGCCCCTCGGGGCGGGGGTCGGGGCGGTTTTCCCCGCGAAGAGGCCAGTACAGGTTTACAt                                                                                              <  2:1982107/147‑1 (MQ=255)
                                                         gcTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACAt                                                                                              >  2:5449670/1‑150 (MQ=255)
                                                         gcTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACCGTCCGCGGGAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACaa                                                                                              >  2:1362935/1‑149 (MQ=38)
                                                           tGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCa                                                                                            >  2:1078155/1‑150 (MQ=255)
                                                                           tGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCggg                                                                                                                               <  2:5872363/99‑1 (MQ=21)
                                                                                   gttttCGTGTCCGTCGGCATTTGACTCCCGCAGCGCCTGGGCCGGCCTCTGCGGGGTTAACCCCGCCCCACCGGTCGGCTGTGGGGGGGGTTTCACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCt                                                                    <  1:1020083/147‑1 (MQ=255)
                                                                                   ggATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCAACGGTCGGCGGGTGGGGCGGGTTTACCCGCGGCGAGGCCACTACAGGTTTACATCAGCCCGAGGCCAGCGCTGCGGAGcg                                                                    >  2:2712766/1‑149 (MQ=38)
                                                                                           gtCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGGAAACCCCCGCCGACGGTCGGCGGGAGGGGCGGGTTTACCCCCGAAGAGGCCACGACAGGTTTACAACAGCCCCCGGCCCGCCCGGCGTAGCGGttaaaca                                                            >  1:4832254/1‑145 (MQ=255)
                                                                                           ggctgccTGCATTTGAATCAGGCTTAGCCTGTGCGGGCCTCTCCGGGGGTAAACCGGCTCCTACGGTCGGCTGTTGGGGGGGTTTTACCCGCGAAGAGGCCAGTTCAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtcatca                                                            <  2:1937866/144‑1 (MQ=255)
                                                                                            tcagtcagCCTTTGAATCACGCATAGCCTGTGCCGGCCTCGTCGCGGGTAAACCCGCTCCTACGCTCGGCGGGAGGGGCGGGGTTACCCCCCACCAGGCCACTACCGGGTTACATCAGCGCGCGGCCACCGCTGCGGAGCGGTTCATCAg                                                           >  2:2826648/1‑150 (MQ=25)
                                                                                                  tGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCTCCTAAGGGCGGCGGTAGGGGCGGGTTTACCCGCGTAGAGGCCAGTACCGGTTTTCATAAGCGCGCGGCCAGCGCTGCGGAGCTGGTCATCAGGTGGCg                                                     >  2:5604638/1‑150 (MQ=25)
                                                                                                     cTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTa                                                  >  2:6781059/1‑150 (MQ=255)
                                                                                                      tttGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGGGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAg                                                 >  2:2364748/1‑150 (MQ=255)
                                                                                                          aaacACGCACCGCCTGGGCCGGCCTCTTCGGGTTTAACCCGGCTCATCGGGCGGGCTGTCGGGGCGGTTTAACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                             <  1:1407420/147‑1 (MQ=255)
                                                                                                          aaTCACGCATCGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                             <  2:2897538/150‑1 (MQ=255)
                                                                                                          aaTCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTgg                                             <  2:1804248/150‑1 (MQ=255)
                                                                                                                        tgtgCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAc                               <  2:16895/150‑1 (MQ=255)
                                                                                                                         gtgCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTTGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAca                              <  1:1729426/150‑1 (MQ=255)
                                                                                                                          tgCCGGCCTCTTCGCGGGTAAAACCCCGCCAACGGTCCGGGGTAGGGGCGGGTTTAACCCCCAAGAGGCCAGTACAGGTTTACAACACCGCGAGGCCAGCGCAGCGGAGCGGTTAAACAGGTTGCGGGAGTTTGGGAAAAGCTGCGCcac                             >  2:5134833/1‑149 (MQ=255)
                                                                                                                            ccGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagca                           <  2:2960960/150‑1 (MQ=255)
                                                                                                                            ccGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagca                           <  2:7304734/150‑1 (MQ=255)
                                                                                                                             cGGCCTCTTCGCGGGTAAACCCGCTCCTCCGGGCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGAcagcag                          <  2:1432707/150‑1 (MQ=255)
                                                                                                                             cGGCCTCTTCGCGGGTAAACCCGCTACTAAGGTAGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTATTTGGGGATACGCTGCGAcagcag                          >  2:538666/1‑150 (MQ=255)
                                                                                                                               gCCTCTTCGCGGGTAAACCCGCTCCAACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATAAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTTGGGATACGCTGCGACAGCAGGt                        >  1:5526398/1‑150 (MQ=255)
                                                                                                                                 ctctTCGCGGCTAACCCCGCCCCTCCGGTCGGCTTTTGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTa                      <  1:414238/150‑1 (MQ=255)
                                                                                                                                      cgcgGGTAACCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAg                 <  1:3560173/150‑1 (MQ=255)
                                                                                                                                       gcgGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGc                >  1:5011272/1‑150 (MQ=255)
                                                                                                                                       gcgGGTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGc                >  2:7264610/1‑150 (MQ=37)
                                                                                                                                         gtttAACCCCGCGCCTTGGGGCGGGTGTTGGGGGGGGTTCACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGccc              <  1:1362935/147‑1 (MQ=255)
                                                                                                                                         gggTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGccc              >  1:525338/1‑150 (MQ=37)
                                                                                                                                          ggTAAACCCGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCt             >  1:4125119/1‑150 (MQ=37)
                                                                                                                                            tAAACCCGCTCCTCCGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTc           <  2:5011272/150‑1 (MQ=255)
                                                                                                                                              accccGCTCCTCCGGGGGGCTGTTGGGGGGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGa         <  1:3995623/148‑1 (MQ=255)
                                                                                                                                               aCCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt        <  2:7209341/150‑1 (MQ=255)
                                                                                                                                                cccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTtc                                                                <  1:3122161/93‑1 (MQ=255)
                                                                                                                                                cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATg       >  2:3125282/1‑150 (MQ=37)
                                                                                                                                                cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                            <  1:2036451/113‑1 (MQ=25)
                                                                                                                                                cccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTggg                                            >  2:2036451/1‑113 (MQ=25)
                                                                                                                                                 ccGCTCCTCCGGTCGGCTGTCGGGGCGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      <  1:804302/150‑1 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      >  2:1599544/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      >  2:333002/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      >  1:5397415/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      >  1:1125305/1‑150 (MQ=255)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt        >  1:70016/1‑148 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGAt        <  2:70016/148‑1 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      >  1:6555981/1‑150 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      >  1:4932602/1‑150 (MQ=38)
                                                                                                                                                 ccGCCCCTACAGCCGACCGTAGGAGCGGGTTTACCCGAGAAGAGGCCAGTACAGGTTTACATCAGCGAGAGTCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGt      >  2:3315130/1‑150 (MQ=37)
                                                                                                                                                  cGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTc     <  1:4900571/150‑1 (MQ=255)
                                                                                                                                                  cGCTCATACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCTCGAGGCCAGCGCTGCGTAGCTCTACAACAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGCTTACACAGCCCTACGAAGTc     >  2:4773372/1‑150 (MQ=255)
                                                                                                                                                    ctcctcCGGTGGGCTGTGGGGGGGGTTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGt   <  1:1953950/150‑1 (MQ=39)
                                                                                                                                                     tCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGtt  >  1:7254460/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATGGGCTTGGCGATCTGTAGTTCGGTGGCTCAGGCCCACGGCGGGCAGCTGCAGGCGCTGGTGTCGCAGTTGGGTGGATGCCGGATTCGTGTCAGTCTGCCTTTGAATCACGCATAGCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGTCGGCTGTAGGGGCGGGTTTACCCGCGAAGAGGCCAGTACAGGTTTACATCAGCGCGAGGCCAGCGCTGCGTAGCTGTTCATCAGGTTGCGGTAGTTGGGGATACGCTGCGACAGCAGGTTACCCAGCCCTTCGATGTCGTT  >  NC_002947/5989292‑5989590

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: