Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,705,753 T→C 5.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,705,7530TC 5.5% 1179.0 / 33.6 363intergenic (‑71/+79)kdpD/kdpCsensor protein/potassium‑transporting ATPase subunit C
Reads supporting (aligned to +/- strand):  ref base T (152/190);  new base C (4/16);  total (157/206)
Fisher's exact test for biased strand distribution p-value = 3.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

Only 100 of 369 total aligned reads displayed.

CCAGGAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGATCA  >  NC_002947/4705604‑4705902
                                                                                                                                                     |                                                                                                                                                     
 cAGGAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTa                                                                                                                                                      >  1:1182637/1‑150 (MQ=255)
 cAGGAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTa                                                                                                                                                      >  2:2525159/1‑150 (MQ=255)
    gAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTAccc                                                                                                                                                   <  2:7746457/150‑1 (MQ=255)
     aaCACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcg                                                                                                                                                  <  2:1521247/150‑1 (MQ=255)
      acacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgc                                                                                                                                                 >  2:883790/1‑150 (MQ=255)
      acacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgc                                                                                                                                                 <  2:4880153/150‑1 (MQ=255)
       cacCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCcgcg                                                                                                                                                <  1:5554256/150‑1 (MQ=255)
         ccTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                                              <  1:8648725/150‑1 (MQ=255)
          cTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAAt                                                                                                                                             >  1:8031577/1‑150 (MQ=255)
          cTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAAt                                                                                                                                             >  1:5031720/1‑150 (MQ=255)
          cTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAAt                                                                                                                                             >  2:4096764/1‑150 (MQ=255)
           ttCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATa                                                                                                                                            >  1:2420005/1‑150 (MQ=255)
            tCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATAc                                                                                                                                           <  1:8216346/150‑1 (MQ=255)
            tCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATAc                                                                                                                                           <  2:4321995/150‑1 (MQ=255)
              aGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGt                                                                                                                                         <  1:1907062/150‑1 (MQ=255)
              aGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGt                                                                                                                                         <  2:5513161/150‑1 (MQ=255)
               gCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTc                                                                                                                                        <  2:5911331/150‑1 (MQ=255)
               gCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCGTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTc                                                                                                                                        >  1:4918699/1‑150 (MQ=255)
                ccggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCa                                                                                                                                       <  2:4803474/150‑1 (MQ=255)
                 cggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAg                                                                                                                                      >  1:6112512/1‑150 (MQ=255)
                  ggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGt                                                                                                                                     <  1:7412779/150‑1 (MQ=255)
                  ggccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGt                                                                                                                                     <  2:1407715/150‑1 (MQ=255)
                    ccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTgg                                                                                                                                   <  2:4560281/150‑1 (MQ=255)
                    ccccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTgg                                                                                                                                   <  2:8381658/150‑1 (MQ=255)
                     cccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGc                                                                                                                                  >  1:4948618/1‑150 (MQ=255)
                     cccggccGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGc                                                                                                                                  <  1:1395104/150‑1 (MQ=255)
                            ggTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGcaccacca                                                                                                                           <  2:1720362/150‑1 (MQ=255)
                                 gTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCa                                                                                                                      <  2:2498691/150‑1 (MQ=255)
                                      aGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTcgc                                                                                                                 <  2:2906768/150‑1 (MQ=255)
                                        gtcggCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                               <  1:759622/147‑1 (MQ=255)
                                        gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                               <  1:7269284/150‑1 (MQ=255)
                                        gTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTGAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                               <  2:5765163/150‑1 (MQ=255)
                                        gTTGGCCAACAGCGCATCGGCGCGCGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCgg                                                                                                               <  2:3795111/150‑1 (MQ=255)
                                              cAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAAc                                                                                                         <  2:3234663/150‑1 (MQ=255)
                                              cAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAAc                                                                                                         <  2:3017279/150‑1 (MQ=255)
                                               aaCAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAAcc                                                                                                        <  1:758343/150‑1 (MQ=255)
                                               aaCAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTGACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTCACCCGCGAATGCGGTGGTGGTGCCACTGACGTATTCGCGGGTAAAcc                                                                                                        <  1:81725/150‑1 (MQ=21)
                                                 cAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCg                                                                                                      <  2:828010/150‑1 (MQ=255)
                                                  aGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGc                                                                                                     <  1:4742085/150‑1 (MQ=255)
                                                    cgcATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTc                                                                                                   >  1:5290321/1‑150 (MQ=255)
                                                    cgcATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTc                                                                                                   >  2:5229749/1‑150 (MQ=255)
                                                       aTCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTa                                                                                                <  1:1973332/150‑1 (MQ=255)
                                                         cGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACa                                                                                              <  1:2111247/150‑1 (MQ=255)
                                                          ggCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAg                                                                                             <  1:1468767/150‑1 (MQ=255)
                                                           gcgcGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAgg                                                                                            <  1:1531305/150‑1 (MQ=255)
                                                              cgggcggAGTCACTCATGTTTCATCCTTGACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCATAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCtt                                                                                         <  2:8074591/150‑1 (MQ=255)
                                                               gggcggAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTc                                                                                        <  2:1229326/150‑1 (MQ=255)
                                                                   ggAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGc                                                                                    >  1:2652091/1‑150 (MQ=255)
                                                                     aGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCa                                                                                  <  1:3156670/150‑1 (MQ=255)
                                                                         aCTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAgcgc                                                                              >  1:927351/1‑150 (MQ=255)
                                                                         aCTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAgcgc                                                                              >  2:756645/1‑150 (MQ=255)
                                                                          cTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCt                                                                             <  2:7983397/150‑1 (MQ=255)
                                                                           tCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCtt                                                                            <  1:5159669/150‑1 (MQ=255)
                                                                           tCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCtg                                                                            >  2:1620625/1‑149 (MQ=255)
                                                                                tttCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGAtt                                                                       >  1:4795352/1‑150 (MQ=255)
                                                                                 ttCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTa                                                                      <  1:7186479/150‑1 (MQ=255)
                                                                                  tCATCCTTTACTTTGTCCAACAGTCCGCCTCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGTGTTTACCCGCGAATACGTCAGTGGCACCCCCACCGAATTCGCGTGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAc                                                                     <  1:3095922/150‑1 (MQ=255)
                                                                                   cATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAcc                                                                    >  2:568680/1‑150 (MQ=255)
                                                                                   cATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTAcc                                                                    <  1:3892897/150‑1 (MQ=255)
                                                                                      ccTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGa                                                                 <  2:6789307/150‑1 (MQ=255)
                                                                                       cTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGAc                                                                >  1:7081879/1‑150 (MQ=255)
                                                                                       cTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACCTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGAc                                                                >  2:2676878/1‑150 (MQ=255)
                                                                                        tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACa                                                               >  2:1940848/1‑150 (MQ=255)
                                                                                        tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACa                                                               >  1:7702645/1‑150 (MQ=255)
                                                                                        tttACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACa                                                               >  2:6394024/1‑150 (MQ=255)
                                                                                         ttACTTTGTACAACAGTCCGCCGCATCAGGGAAAGCCCCCTTAAAGTGTATGAGCGGGTTTACACGCCAAAACGTCAGTGGCACCACCACCCCAATCGCGGGTGAAACCGCTCCTAAACGCTTCCTGCCACCGCTTCGATAACCAGACaa                                                              >  2:2075379/1‑150 (MQ=255)
                                                                                           aCTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATg                                                            <  1:5234489/150‑1 (MQ=255)
                                                                                             tttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCCGCGCTTCGATTACCTGACAATGGc                                                          >  1:3449762/1‑150 (MQ=255)
                                                                                             tttGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGc                                                          <  2:2284634/150‑1 (MQ=255)
                                                                                               tGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGa                                                        >  2:377001/1‑150 (MQ=255)
                                                                                                gTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAg                                                       >  1:7854112/1‑150 (MQ=255)
                                                                                                   cAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCtt                                                    >  1:17665/1‑150 (MQ=255)
                                                                                                     aCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCtttt                                                  <  2:598231/150‑1 (MQ=255)
                                                                                                     aCAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCtttt                                                  <  1:5827948/150‑1 (MQ=255)
                                                                                                       aGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTcc                                                <  2:5403148/150‑1 (MQ=255)
                                                                                                         tCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCaa                                              <  2:1079789/150‑1 (MQ=255)
                                                                                                         tCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCaa                                              <  2:4161496/150‑1 (MQ=255)
                                                                                                          ccgccgCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAg                                             <  1:2963366/150‑1 (MQ=255)
                                                                                                            gccgcTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGc                                           >  2:8682603/1‑150 (MQ=255)
                                                                                                                cTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTgg                                       >  1:5374708/1‑150 (MQ=255)
                                                                                                                  tCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGACCCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGtt                                     <  1:6816260/150‑1 (MQ=255)
                                                                                                                  tCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGtt                                     <  1:4279822/150‑1 (MQ=255)
                                                                                                                    aGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTc                                                                                                                   <  2:506285/70‑1 (MQ=255)
                                                                                                                     gggAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCaa                                  <  2:7080570/150‑1 (MQ=255)
                                                                                                                      ggAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAAc                                 <  1:6630499/150‑1 (MQ=255)
                                                                                                                       gAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACg                                >  1:6544952/1‑150 (MQ=255)
                                                                                                                         aaGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGcc                              >  1:1070951/1‑150 (MQ=255)
                                                                                                                             ccccTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCa                          >  1:6180959/1‑150 (MQ=255)
                                                                                                                                     cGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGacc                  >  2:7226793/1‑150 (MQ=255)
                                                                                                                                      gtgtAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCCAGGCCTGGTTCAACGCCAGCACGTTGacca                 >  2:6627769/1‑150 (MQ=255)
                                                                                                                                       tgtAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGaccac                >  1:8728829/1‑150 (MQ=255)
                                                                                                                                       tgtAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGaccac                >  1:841461/1‑150 (MQ=255)
                                                                                                                                        gtAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGaccacc               <  1:6367212/150‑1 (MQ=255)
                                                                                                                                         tAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcg              <  2:729852/150‑1 (MQ=255)
                                                                                                                                         tAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcg              >  1:449591/1‑150 (MQ=255)
                                                                                                                                          aGGAGCGGGTTTACCCGCGAATACGTCAGTGGCGCCACCACCGCATTCGCGGGTGAACCCGCTCCGACAGGCTTCCTGCCAGCGCGTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGGTCAACGCCCGCACGTTGACCACcgg             >  1:3740285/1‑150 (MQ=38)
                                                                                                                                            tcgCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGc                                                                                    >  2:6164811/3‑77 (MQ=255)
                                                                                                                                            gAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggcg           >  2:6338744/1‑150 (MQ=255)
                                                                                                                                             aGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACcggcgg          <  1:4096764/150‑1 (MQ=255)
                                                                                                                                                     tACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGATCa  <  2:3974581/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGGAACACCTTCAGCCGGCCCCGGCCGGTGCGTGGGAGGTTGGCCAACAGCGCATCGGCGCGGGCGGAGTCACTCATGTTTCATCCTTTACTTTGTCCAACAGTCCGCCGCTTCAGGGAAAGCCCCCTTAACGTGTAGGAGCGGGTTTACCCGCGAATACGTCAGTGGCACCACCACCGCATTCGCGGGTGAACCCGCTCCTACAGGCTTCCTGCCAGCGCTTCGATTACCTGACAATGGCGAGCTTTTCCAAGGCCTGGTTCAACGCCAGCACGTTGACCACCGGCGGGCCGATCA  >  NC_002947/4705604‑4705902

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: