Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,702,726 T→C 8.5% intergenic (‑71/+79) kdpD ← / ← kdpC sensor protein/potassium‑transporting ATPase subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,702,7260TC 8.5% 492.2 / 23.0 165intergenic (‑71/+79)kdpD/kdpCsensor protein/potassium‑transporting ATPase subunit C
Reads supporting (aligned to +/- strand):  ref base T (81/68);  new base C (3/11);  total (84/81)
Fisher's exact test for biased strand distribution p-value = 2.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 166 total aligned reads displayed.

AGGAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGAT  >  NC_002947/4702579‑4702873
                                                                                                                                                   |                                                                                                                                                   
  gAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACATTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTAccc                                                                                                                                                 <  2:4685940/150‑1 (MQ=255)
    acacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgc                                                                                                                                               <  1:4694721/150‑1 (MQ=255)
    acacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgc                                                                                                                                               <  1:5115492/150‑1 (MQ=255)
    acacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgc                                                                                                                                               <  2:1872772/150‑1 (MQ=255)
      acCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCAGGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGa                                                                                                                                             >  2:2348245/1‑150 (MQ=255)
       ccTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                            <  2:22108/150‑1 (MQ=255)
       ccTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                            <  2:2750278/150‑1 (MQ=255)
       ccTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTAAACCcgctcc                                                                                                                                            <  2:1704838/150‑4 (MQ=255)
        cTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAAt                                                                                                                                           >  1:1473118/1‑150 (MQ=255)
         ttCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATa                                                                                                                                          <  2:2501434/150‑1 (MQ=255)
         acgcaccggccggggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATa                                                                                                                                          >  1:2273438/6‑150 (MQ=255)
           gcaccggccggggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACg                                                                                                                                        >  2:404500/4‑150 (MQ=255)
           gcaccggccggggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACg                                                                                                                                        >  1:829424/4‑150 (MQ=255)
            aGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGt                                                                                                                                       <  2:1392778/150‑1 (MQ=255)
              ccggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCa                                                                                                                                     <  1:1779781/150‑1 (MQ=255)
              ccggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTGTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTCACGTGTAGGAGCGGGTTCACCCGCGAATGCggtg                                                                                                                                     <  2:3397891/150‑4 (MQ=255)
                ggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGt                                                                                                                                   <  1:2008165/150‑1 (MQ=255)
                 gccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTg                                                                                                                                  >  2:1508593/1‑150 (MQ=255)
                 gccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATACGTCAGTg                                                                                                                                  <  2:4540800/150‑1 (MQ=255)
                   cccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGc                                                                                                                                >  2:4497940/1‑150 (MQ=255)
                    ccggccggTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGca                                                                                                                               <  2:3143848/150‑1 (MQ=255)
                     cggccggTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcac                                                                                                                              >  2:5175318/1‑150 (MQ=255)
                      ggccggTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcacc                                                                                                                             >  1:843605/1‑150 (MQ=255)
                                   gAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTcg                                                                                                                >  1:3288784/1‑150 (MQ=255)
                                    aGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTcgc                                                                                                               <  2:5066579/150‑1 (MQ=255)
                                     ggTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTcgcg                                                                                                              <  1:3639038/150‑1 (MQ=255)
                                      gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                             <  2:649690/150‑1 (MQ=25)
                                      gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                             >  2:448147/1‑150 (MQ=255)
                                      gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                             <  1:4685037/150‑1 (MQ=255)
                                      gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                             <  1:2994429/150‑1 (MQ=25)
                                      gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCgg                                                                                                             <  1:1103490/150‑1 (MQ=25)
                                      gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                             <  1:2745718/150‑1 (MQ=255)
                                         ggCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTg                                                                                                          >  1:4273475/1‑150 (MQ=255)
                                           ccAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGaa                                                                                                        >  2:2076042/1‑150 (MQ=255)
                                            cAACAGCGCATCGGCGCGGGCGGAGTCACTCCTGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAAc                                                                                                       <  2:1123805/150‑1 (MQ=21)
                                            cAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAAc                                                                                                       <  2:4003622/150‑1 (MQ=21)
                                                  cgcATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTc                                                                                                 <  2:2530221/150‑1 (MQ=255)
                                                             gggcggAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTc                                                                                      >  1:4048673/1‑150 (MQ=255)
                                                               gcggAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCt                                                                                    >  1:3431308/1‑150 (MQ=255)
                                                                   aGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCa                                                                                <  2:2499952/150‑1 (MQ=255)
                                                                       aCTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAgcgc                                                                            >  2:670915/1‑150 (MQ=255)
                                                                           aTGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCg                                                                        >  2:3359538/1‑150 (MQ=255)
                                                                              tttCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGAtt                                                                     >  1:1972313/1‑150 (MQ=255)
                                                                               ttCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTa                                                                    >  1:946200/1‑150 (MQ=255)
                                                                                tCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAc                                                                   <  1:3481145/150‑1 (MQ=255)
                                                                                tCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAc                                                                   >  1:4303040/1‑150 (MQ=255)
                                                                                 cATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAcc                                                                  <  1:1262619/150‑1 (MQ=255)
                                                                                     cTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGAc                                                              >  2:4068993/1‑150 (MQ=255)
                                                                                     cTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAAccc                                                                                                     <  2:1407103/111‑1 (MQ=14)
                                                                                     cTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAAccc                                                                                                     >  1:1407103/1‑111 (MQ=14)
                                                                                      tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAgc                                                                              <  2:3548305/133‑1 (MQ=255)
                                                                                      tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAgc                                                                              >  1:3548305/1‑133 (MQ=255)
                                                                                           tttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGc                                                        >  1:234005/1‑150 (MQ=255)
                                                                                           tttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGc                                                        >  2:4271191/1‑150 (MQ=255)
                                                                                           tttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCGTCGATTACCTGACAATGGc                                                        >  2:3436986/1‑150 (MQ=255)
                                                                                            ttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCg                                                       >  2:1604688/1‑150 (MQ=255)
                                                                                              gTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAg                                                     <  1:4311310/150‑1 (MQ=255)
                                                                                               tCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGc                                                    >  2:3807653/1‑150 (MQ=255)
                                                                                                  aaCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCttt                                                 >  2:3454731/1‑150 (MQ=255)
                                                                                                  aaCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCttt                                                 >  1:1042412/1‑150 (MQ=255)
                                                                                                  aaCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCttt                                                 >  2:2022384/1‑150 (MQ=255)
                                                                                                  aaCAGTCAGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCACCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCttt                                                 >  2:2735349/1‑150 (MQ=255)
                                                                                                   aCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCtttt                                                <  2:253531/150‑1 (MQ=255)
                                                                                                    cAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTc                                               >  2:1865781/1‑150 (MQ=255)
                                                                                                    cAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTc                                               >  2:1681798/1‑150 (MQ=255)
                                                                                                    cAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTc                                               <  2:4891245/150‑1 (MQ=255)
                                                                                                     aGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTcc                                              <  2:4446675/150‑1 (MQ=255)
                                                                                                     aGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTcc                                              <  1:2981762/150‑1 (MQ=255)
                                                                                                       tCCGCCGCTTCAGCGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCaa                                            <  1:507279/150‑1 (MQ=255)
                                                                                                        ccgccgCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAg                                           >  2:1891739/1‑150 (MQ=255)
                                                                                                          gccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGc                                         >  1:4919570/1‑150 (MQ=255)
                                                                                                          gccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGc                                         >  1:4919584/1‑150 (MQ=255)
                                                                                                          gccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGc                                         >  2:2658881/1‑150 (MQ=255)
                                                                                                          gccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCCAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGc                                         >  1:3163471/1‑150 (MQ=255)
                                                                                                                 cAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTc                                  <  1:201948/150‑1 (MQ=255)
                                                                                                                   gggAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCaa                                <  1:1397161/150‑1 (MQ=255)
                                                                                                                    ggAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCtt                                                                          <  1:3013958/107‑1 (MQ=255)
                                                                                                                    ggAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAAc                               <  2:3291193/150‑1 (MQ=255)
                                                                                                                    ggAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAAc                               <  2:4273475/150‑1 (MQ=255)
                                                                                                                      aaaGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGc                             <  1:5190839/150‑1 (MQ=255)
                                                                                                                        aGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCa                           >  1:1169449/1‑150 (MQ=255)
                                                                                                                          cccccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTa                                                                    >  1:247759/1‑107 (MQ=255)
                                                                                                                          cccccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTa                                                                    <  2:247759/107‑1 (MQ=255)
                                                                                                                          cccccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGc                         >  2:4660868/1‑150 (MQ=255)
                                                                                                                          cccccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGc                         >  2:4660856/1‑150 (MQ=255)
                                                                                                                            cccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCa                                                                                                                    <  2:3150368/57‑1 (MQ=255)
                                                                                                                               ttAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGtt                    <  1:1604189/150‑1 (MQ=255)
                                                                                                                               ttAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGtt                    <  1:1604688/150‑1 (MQ=255)
                                                                                                                                 aaCGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAAc                               <  2:765384/137‑1 (MQ=255)
                                                                                                                                 aaCGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAAc                               >  1:765384/1‑137 (MQ=255)
                                                                                                                                 aaCGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGa                  <  1:3935094/150‑1 (MQ=255)
                                                                                                                                     tgtAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGaccac              <  1:1098010/150‑1 (MQ=255)
                                                                                                                                      gtAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGaccacc             >  2:1562379/1‑150 (MQ=255)
                                                                                                                                         ggAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggc          >  1:551578/1‑150 (MQ=255)
                                                                                                                                           aGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggcgg        <  1:1646747/150‑1 (MQ=255)
                                                                                                                                            gCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCggg       >  2:3095100/1‑150 (MQ=255)
                                                                                                                                              gggTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGcc     >  1:1837418/1‑150 (MQ=21)
                                                                                                                                                gTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGa   >  2:1339223/1‑150 (MQ=255)
                                                                                                                                                gTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGa   <  1:4390535/150‑1 (MQ=255)
                                                                                                                                                 tttACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGAt  <  1:4625779/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGGAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGAT  >  NC_002947/4702579‑4702873

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