Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,707,381 T→A 11.3% intergenic (‑151/+26) kdpF ← / ← PP_4162 potassium‑transporting ATPase subunit F/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,707,3810TA11.3% 272.0 / 22.5 105intergenic (‑151/+26)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (44/46);  new base A (7/5);  total (54/51)
Fisher's exact test for biased strand distribution p-value = 7.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

Only 100 of 108 total aligned reads displayed.

GTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAAC  >  NC_002947/4707233‑4707530
                                                                                                                                                    |                                                                                                                                                     
gtgGATCTCCGCGGTTTCCGGCTGGGGCCGCTTTGGGGATATTGTCCGAGGGTGAGGCGTAATGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCtt                                                                                                                                                      <  2:4311526/150‑1 (MQ=255)
    aTCTCCCCGTGGTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTGGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGccc                                                                                                                                                  >  2:1637072/1‑150 (MQ=255)
      cgccccGTGTCCCGGCGGGGGGCGCTGGGGGGACCTGGGCCGCGGGCGCGGCGCACCGGGGGGCGCGCGGGGGGGGGGGGCGGGTATAAAGAATGCGGAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCcgc                                                                                                                                                <  2:166565/148‑1 (MQ=255)
       tcCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGCTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCcgcg                                                                                                                                               >  2:2033345/1‑150 (MQ=255)
        ccccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGa                                                                                                                                              >  2:550668/1‑150 (MQ=255)
         cccGTTTTTCGGCTTGGGCGTCTTGGGGGATATTGTGCGCGGGAGTGGGTTACTGGGGCTTGATGGGGGGGTGGGGGTGGTTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCAAGCCCGCGaa                                                                                                                                             <  2:883219/150‑1 (MQ=255)
         cccGTGTTCCGGCTTGGGCCGCTGGGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGGATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             >  1:1144784/1‑150 (MQ=255)
         cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             >  1:5119206/1‑150 (MQ=255)
           cttgtTCCGGCTTGGGGCGCTTGGGGGATTTGGTCGGGGGGCGCGGCGTCCCGGGGCGCGGTGGGGGGGTGGGGGGGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaaga                                                                                                                                           <  2:4101919/148‑3 (MQ=16)
           cGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                           >  2:3229535/1‑150 (MQ=255)
            gtgtTCCGGCTGGGGCCGCGTGGGGGACCGTGCCCGCGGGCGAGGCGTCCCGGGGCGCGCCCGGGGGGGGGGCCCGCCCCTAAAGAAGGCGGAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCAAGCCCGCGaagag                                                                                                                                          <  2:397333/150‑4 (MQ=255)
             tgttcgggCTGGGGCCGGTTTGGGTCACGGGTCCGCGGGGGCGGGGTCCGGGGGCGCGCGCGGGGGGGGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAggg                                                                                                                                         <  1:4959507/144‑1 (MQ=255)
              gtgtcGGCTTGGGGCGCGGGGGGGCACTGGGCCGAGGGCGCGGCGTCCCGGGGGGAGACGGGGGGGGGGGGGGGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGc                                                                                                                                        <  2:4151694/146‑1 (MQ=255)
               ttACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                                                                                                       >  2:2172314/1‑150 (MQ=255)
                 cgggcgggggCCGGTTGGGGGTAATTGTCGGCGGGAGAGCCGTACGGGGGGGATCGGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGt                                                                                                                                     <  2:1691000/143‑1 (MQ=255)
                   ggCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTa                                                                                                                                   >  1:3361825/1‑150 (MQ=255)
                   ggCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTa                                                                                                                                   >  1:3361852/1‑150 (MQ=255)
                    gCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTAc                                                                                                                                  >  2:1594025/1‑150 (MQ=255)
                        gggCCGCTTTGCGGACTTTGTCGGTGGGCGTGGCGGCAGGGGGCGCGCTGGGGGGGGGGGGGCGCCAATAAACAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                                                                                              <  2:326185/150‑1 (MQ=255)
                        gggCCGCTTGGGTTCATTTGTCCTCGGGCGAGGCGTAACGGGGCGAGCGGGGGGGGTGGGGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                                                                                              <  2:455773/150‑1 (MQ=255)
                           ccGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGTTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTa                                                                                                                           >  1:2219821/1‑150 (MQ=255)
                            cGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTAt                                                                                                                          >  1:3197815/1‑150 (MQ=255)
                              cTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttg                                                                                                                        >  1:391296/1‑150 (MQ=255)
                               tttGGGGTCCTTGTCGGCGGGTGCGGCGCACGGGGGCGGGCGGGGGGGGTGGGGCCGCGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgt                                                                                                                       <  2:391296/150‑1 (MQ=255)
                                    gtcacgtGTCGGCGGGCGAGGCGCACGGGGGCGAGCGCGGGGGGGGGGGGTGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                                                                                                  <  1:4458559/144‑1 (MQ=255)
                                    ggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                                                                                                  >  1:4531287/1‑150 (MQ=255)
                                              cGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCtgctgc                                                                                                        >  2:2526089/1‑150 (MQ=255)
                                              cGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTAATTGTTGTGGCGGCtgctgc                                                                                                        >  2:728386/1‑150 (MQ=255)
                                              cGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCACTACGCGGGCTTGCCAGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCtgctgc                                                                                                        >  2:728373/1‑150 (MQ=255)
                                                          ggaaaGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTACGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACgg                                                                                            >  2:269405/3‑150 (MQ=255)
                                                             ccggggCGCGCCCGGGGGGTGGCGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc                                                                                         <  1:2561975/148‑1 (MQ=255)
                                                              cggggCGCGTTCGGGGGGTGGGGGAGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcc                                                                                        <  2:2436370/149‑1 (MQ=255)
                                                               ggggCGCGCGCGGGGGGGGGGGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                       <  2:3856839/150‑1 (MQ=255)
                                                               ggggCGCGCCCGGGGGGGGGGGGCGGCAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                       <  1:3179126/150‑1 (MQ=255)
                                                                gggCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccag                                                                                      >  2:1451291/1‑150 (MQ=255)
                                                                  gCGAGATCGGGGGGAGGGGGTGGGACTAAAGAATGCGTAAAAATCAGCGGGGTGGTCTGCCTGTAACGGCCTCTGCGCGCGCTTGCCCGCGAAAGGGCCGGTACAGGCGGCGGGGTGAGGCGGCTGCGGCGCCGCGGGCCGGCtccagac                                                                                    >  2:3268585/1‑148 (MQ=255)
                                                                    gcgcggggggggTGGGGGGGGTAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCGGCCGGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  <  1:2274552/144‑1 (MQ=255)
                                                                    gcgcgcGGGGGGGGGGGGGGGGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                  <  1:3849044/145‑1 (MQ=255)
                                                                    gagaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTTTGGCGGCTGCTGCGCCGCGTTACGGctccactcca                                                                                  >  1:614119/1‑150 (MQ=255)
                                                                     cgctCGGGGGGGGGGGGCGGGTCGAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtccag                                                                                 <  2:1596950/147‑1 (MQ=255)
                                                                      gcggggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagc                                                                                <  2:2455522/146‑1 (MQ=255)
                                                                      gaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagc                                                                                >  2:573899/1‑150 (MQ=255)
                                                                       aTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAAAGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGGCGTTACAGGCTTATTTTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTgcagca                                                                               >  1:2667283/1‑150 (MQ=255)
                                                                          ggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCGTGCCCCCGAAAGGGCCGTTACACGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGACTCCAGTCcagcagca                                                                            >  1:1899628/1‑150 (MQ=255)
                                                                          ggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGATGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagcagca                                                                            >  2:3927849/1‑150 (MQ=255)
                                                                                   gggcggctATAAAGAAGGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGc                                                                   <  2:3143091/150‑1 (MQ=255)
                                                                                    ggtgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGcc                                                                  >  2:1828008/1‑150 (MQ=255)
                                                                                      tgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTa                                                                >  1:2674004/1‑150 (MQ=255)
                                                                                        ggtATAAAGAACGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAcc                                                              <  2:494524/150‑1 (MQ=255)
                                                                                        ggaagaaaGAACGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAcc                                                              <  1:1410472/145‑1 (MQ=255)
                                                                                        gctataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAcc                                                              <  2:2984401/148‑1 (MQ=255)
                                                                                        gccataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAcc                                                              <  2:494578/147‑1 (MQ=255)
                                                                                          tataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGc                                                            <  2:5119206/150‑1 (MQ=255)
                                                                                          cagaaaGCATGCGTAACAATCAGCGTGGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGc                                                            <  2:804154/147‑1 (MQ=255)
                                                                                              aaGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGc                                                        >  2:2020716/1‑150 (MQ=255)
                                                                                              aaGAATGCGTAAAAATCAGCGTTGTGGTCTGCCAGTAACCGCCACTTCGCGGGCTTGCCCGCGCAAGGGCCGTTACAGGGGTATTGTTGTGGCGGCTGCGGCGCCGCCTTACCGCGCCACTCCAGCCGCAGGCGGTAGCCAACccccccc                                                        >  1:1850933/1‑147 (MQ=255)
                                                                                                gAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGCAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTAAAGGCGTATTGGTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCac                                                      >  1:200634/1‑149 (MQ=255)
                                                                                                    gCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTTGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTc                                                  <  1:4480070/150‑1 (MQ=255)
                                                                                                    gCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTc                                                  <  1:1538541/150‑1 (MQ=255)
                                                                                                      gTAAAACTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCgg                                                <  1:3229535/150‑1 (MQ=255)
                                                                                                        aaaaaTCAGCGTTGGGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGcc                                              <  1:3871912/150‑1 (MQ=255)
                                                                                                         aaaaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGccc                                             <  1:4192772/150‑1 (MQ=255)
                                                                                                           aaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                           <  2:2928177/150‑1 (MQ=255)
                                                                                                            aTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                          <  2:2039750/150‑1 (MQ=255)
                                                                                                             tCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAATGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCCCCAGCAAGGTCGGCCCGATg                                         >  1:1144603/1‑150 (MQ=255)
                                                                                                                 cGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                     <  2:1413571/150‑1 (MQ=255)
                                                                                                                 cGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                     <  2:1144784/150‑1 (MQ=255)
                                                                                                                    tgtgGTCTGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                                  <  1:4664128/150‑1 (MQ=255)
                                                                                                                    tgtgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAGGGc                                  <  2:1888132/150‑1 (MQ=255)
                                                                                                                    tgtgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                                  <  1:1234217/150‑1 (MQ=255)
                                                                                                                      tgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCa                                <  1:2526089/150‑1 (MQ=255)
                                                                                                                        gTCTGCCTGTAACGGCCTCTTCGCGGGCGTGCCCGCGACAGGGCCGTTCCCGGCGTAGGGGTGTGGCGGCTGCGGCGCCGCGGTACGGCTCCAGTCCCGCAGCAGGCGGGAGCCTCCCGCCCGCCCGGGCGGCCCGATGAAAAGGCCCAt                              >  2:2170890/1‑150 (MQ=255)
                                                                                                                           tGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                           <  1:1594025/150‑1 (MQ=255)
                                                                                                                           tGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                           >  2:1596161/1‑150 (MQ=255)
                                                                                                                            gCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAg                          <  1:1637072/150‑1 (MQ=255)
                                                                                                                              cTGTAACGGCCTCTTCGCGGGCTTTCCAGCGAAAGGGCCGTTACAGGCGTATTGTTGTGGCGGATGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCCATGAAGcc                        >  1:3928812/1‑150 (MQ=255)
                                                                                                                              cTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                        >  1:1662773/1‑150 (MQ=255)
                                                                                                                              cTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                        >  2:4265478/1‑150 (MQ=38)
                                                                                                                              cTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAGGAAAAGGCCAATGAAGcc                        >  2:60125/1‑150 (MQ=38)
                                                                                                                              cTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCGTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                        >  1:1689910/1‑150 (MQ=38)
                                                                                                                                 tAACGGCCTCTTCGCGGGCTGGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCaaa                     <  2:126324/150‑1 (MQ=255)
                                                                                                                                 tAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCaaa                     >  2:4946638/1‑150 (MQ=38)
                                                                                                                                   aCGGCCTCTTCTCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTATTGTTGTTGCGGCTGCTGCGCCGCGTTACAGCTCCCGTCCAGCAGCAGGCAGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAATGCCAATGAAGCCAAAgg                   >  1:858289/1‑150 (MQ=255)
                                                                                                                                    cGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgc                  >  1:3369861/1‑150 (MQ=255)
                                                                                                                                     ggccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgca                 >  2:1689440/7‑150 (MQ=38)
                                                                                                                                     ggCCTCTTCGCGGGCGTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgca                 >  2:2887626/1‑150 (MQ=255)
                                                                                                                                      gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                >  2:2723196/6‑150 (MQ=39)
                                                                                                                                      gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                >  1:1604545/6‑150 (MQ=39)
                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTa                                                                >  1:1344799/1‑102 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTa                                                                <  2:1344799/102‑1 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                <  2:2219821/150‑1 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTCGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                >  2:1374178/1‑150 (MQ=255)
                                                                                                                                        ccctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc              >  1:4826140/3‑150 (MQ=255)
                                                                                                                                         tctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagca             >  1:627839/1‑150 (MQ=255)
                                                                                                                                          ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagcaa            <  1:2172314/150‑1 (MQ=255)
                                                                                                                                          ctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAACCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagcaa            <  2:4779295/150‑1 (MQ=255)
                                                                                                                                               cgGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccgc       <  2:4531287/150‑1 (MQ=255)
                                                                                                                                                gggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccgcc      >  2:1637016/1‑150 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAac  >  1:4676116/1‑150 (MQ=255)
                                                                                                                                                    ttGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAac  >  1:2525000/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAAC  >  NC_002947/4707233‑4707530

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