Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,707,302 A→C 11.9% intergenic (‑72/+105) 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,3020AC11.9% 602.9 / 4.9 252intergenic (‑72/+105)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (130/67);  new base C (2/35);  total (133/119)
Fisher's exact test for biased strand distribution p-value = 1.38e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

GGCGCGCAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAG  >  NC_002947/4707153‑4707451
                                                                                                                                                     |                                                                                                                                                     
ggCGCGCAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCga                                                                                                                                                       <  1:7600324/150‑1 (MQ=255)
ggCGCGCAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCga                                                                                                                                                       <  2:4255575/150‑1 (MQ=255)
ggCGCGCAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCga                                                                                                                                                       <  2:3522044/150‑1 (MQ=255)
 gcgcgcAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCgag                                                                                                                                                      <  2:552633/150‑1 (MQ=255)
   gcgcAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAt                                                                                                                                                    >  1:8279917/1‑150 (MQ=255)
      cagcagCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCgg                                                                                                                                                 >  2:2134047/1‑150 (MQ=255)
        gcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCgggg                                                                                                                                               <  1:593968/150‑1 (MQ=255)
              ccaccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtggg                                                                                                                                         <  1:2371049/150‑1 (MQ=255)
               caccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtgggg                                                                                                                                        <  2:5750856/150‑1 (MQ=255)
               caccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtgggg                                                                                                                                        >  2:663202/1‑150 (MQ=255)
                  cagcagATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGggggtgggggtg                                                                                                                                     <  1:1629177/150‑1 (MQ=255)
                   agcagATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGgggtgggggtgg                                                                                                                                    <  2:303895/150‑1 (MQ=255)
                    gcagATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGCGCTTGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGggtgggggtggg                                                                                                                                   <  2:194995/150‑1 (MQ=255)
                    gcagATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGggggggggtggg                                                                                                                                   >  2:1536656/1‑150 (MQ=255)
                         tAACACGCCAAGGCCACTGCCAGCAGCAGTGACCGCCCGTCGAGCATGTTCATGTGTGGATCTCCCGGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCGGCGGGAGCGGGGTACAGGGGCGTGATCGGGGGGTGGGGGTGGGTATaa                                                                                                                              <  2:320469/150‑1 (MQ=255)
                              cGCCTATGCCACTGCCAGCAGCAGTGACAGCCCTTGTAGCTTGTTCTTGTGTGGATGTCCCCGTGTTAGGGCTTGGGCGGGTTTGGTGAAATTGTCGGTGGGAGAGGCGTATAGGGGCGTGTTCGGGGGGTGGGGGTGGGTATAAAGAAt                                                                                                                         <  2:5476171/150‑1 (MQ=255)
                              cGCAAAGGCCACTGCCAGCTGCTGTGACAGCCGGTCGTGCTTGTTGTTTTGTGGTTCTCCCGGTTTTAGGGGTTGGGCGGCTTTGGGTAATTTGTCCGTGGGTGTGGCTTAACGGGGGGTGATGGGGGGGTGGGGGTGGGTATAAAGAAt                                                                                                                         <  2:8612770/150‑1 (MQ=255)
                               gCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATg                                                                                                                        >  2:6977934/1‑150 (MQ=255)
                                ccAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGc                                                                                                                       >  2:824629/1‑150 (MQ=255)
                                ccAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGGGGGGGGGGGTATAAAGGATGc                                                                                                                       >  2:3810030/1‑150 (MQ=255)
                                  aaTGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGt                                                                                                                     >  1:3625107/1‑150 (MQ=255)
                                        aCTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAAt                                                                                                               >  1:1367757/1‑150 (MQ=255)
                                        aCTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAAt                                                                                                               >  2:1630291/1‑150 (MQ=255)
                                         cTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATc                                                                                                              >  1:2166682/1‑150 (MQ=255)
                                          tGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATTGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCa                                                                                                             >  2:5659130/1‑150 (MQ=255)
                                              agcagcagTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGt                                                                                                         >  1:507477/1‑150 (MQ=255)
                                              agcagcagTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGt                                                                                                         >  1:8384043/1‑150 (MQ=255)
                                                 aggtgTGACAGGCCGTCGTGGATGTTCTTGTGTGGTTCTCGGCGTGTTAGGGCTTGGGCGGGTTTGGGGAATTTGTCGGAGGGAGAGGGGTAAAGGGGGGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTtgt                                                                                                      <  2:7456107/146‑1 (MQ=255)
                                                  gttgTGACTGCCGGTGGCGCCTGTTTTTGTGTGTTTCGCCCTTTGTTTGGGGTTGGGGGGGTTTGGGTACATTGTCGGTGGTTGTGGCTTACCGGGGCGCGTTGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTtgtg                                                                                                     <  2:4091/147‑1 (MQ=255)
                                                    agTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGt                                                                                                   >  1:5472708/1‑150 (MQ=255)
                                                    agTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGGTGTggg                                                                                                   >  2:2232748/1‑149 (MQ=255)
                                                      tGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGTGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCt                                                                                                 >  1:4557480/1‑150 (MQ=255)
                                                       gACCGCCCTTCGAGTTTGTTCTTGTGTGGATCTCGCGGTTTTTGGGTTTGGGCCGCTTTGGGGAATTTTTCCTTGGGAGTGGCTTATCGGGGCGATCTGGGGGGGTGGGGGTGTGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTg                                                                                                <  2:3265265/150‑1 (MQ=255)
                                                        cccgCCCGTCTAGCATGTTCTGGTGTGGATCTCCCCGTGTTAGGGCTTGGGCCGCTTTGGGGAATTTGTCCGAGGGAGAGGCGTAACGGGGCGAGCTGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGc                                                                                               <  2:6731915/147‑1 (MQ=255)
                                                        ccAGCCCGTCGAGCTTGTTCTTGTGTGGATCTCCCCGTGTTAGGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGCGTGGCGTACCGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGc                                                                                               <  2:5999482/149‑1 (MQ=255)
                                                          cgCCCGTCGTGCTTTTTCATGTGTGTATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGTATATTGTCCGAGGGAGAGGCGTATCGGGGCGAGATCGGGGGGTGGGGGTGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCt                                                                                             <  2:3670031/149‑1 (MQ=255)
                                                          aGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCt                                                                                             >  1:1031664/1‑150 (MQ=255)
                                                            cccGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGt                                                                                           >  1:7675852/1‑150 (MQ=255)
                                                                  gcgcttGTTCATGTGTGTATCTCCCGGTGTTCGGGCTTGGGCGGCTTTGGGGAATTTGCCGGTGGGTGAGGCGTAACGGGGCGAGATGGGGGGGTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGc                                                                                     <  1:6310281/145‑1 (MQ=255)
                                                                  gAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTATAACGGc                                                                                     >  1:43496/1‑150 (MQ=255)
                                                                    gCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCct                                                                                   >  2:1189808/1‑150 (MQ=255)
                                                                       tGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCtt                                                                                >  2:2927510/1‑150 (MQ=255)
                                                                            ttgtgtgGCTCTCCCGGTGTTCCGGGTTGGGCCGCGTTGGGGATATTGTCCGTGGGCGGGGCGTCCGGGGGCGCGCCGGGGGGGTGGGGGGGGGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCgg                                                                           <  1:8848142/149‑1 (MQ=255)
                                                                            tggtgtgtATCCCCCGGGGTTCGGGCTGGGGCGGGTTTGGGGAAGGGGCCGGCGGGCGCGGCGTCCCGGGGCGCGCGGGGGGGGGGGGGGGGGCAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCgg                                                                           <  1:4611771/148‑1 (MQ=255)
                                                                            aTGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCgg                                                                           >  1:3778766/1‑150 (MQ=255)
                                                                              gtgtgGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGc                                                                         >  1:2975647/1‑150 (MQ=255)
                                                                                gtgGATCTCCGGGTGTTCGGGCGTGGGCCGCTTTGTGGTATTTGGCCGAGGGTGCGGCTTATTGGGGCGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCtt                                                                       <  2:5424575/150‑1 (MQ=255)
                                                                                   gATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGcc                                                                    >  2:3109388/1‑150 (MQ=255)
                                                                                   gATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGcc                                                                    >  1:6908633/1‑150 (MQ=255)
                                                                                    aTCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGccc                                                                   >  1:5195302/1‑150 (MQ=255)
                                                                                     tctcCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCcg                                                                  >  2:4265959/1‑150 (MQ=255)
                                                                                      ctcCCCGTGTTCCGGCTTGGGCGGCTTTGGGGAAATTGTCCGCGGGAGAGGCGTAACGGGGCGCGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCcgc                                                                 <  2:8020834/150‑1 (MQ=255)
                                                                                      ctcCCCGTGTTCCGGCTTGGGCCGCTTTGGGGACCTGGTCCGTGGGCGCGGCGTAACGGGGCGCGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCcgc                                                                 <  2:3406466/150‑1 (MQ=255)
                                                                                         cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                              >  2:3748628/1‑150 (MQ=255)
                                                                                         cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                              >  1:6087660/1‑150 (MQ=255)
                                                                                         cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCGTGCCCCCGaa                                                              >  1:3831467/1‑150 (MQ=255)
                                                                                           cGTGTTCGGGCTTGGGCGGCTTGGGGTAATTGGCCCGCGGGCGCGGCGTAGGGGGGCGCGAGGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaaga                                                            <  1:3243667/150‑3 (MQ=19)
                                                                                           cGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                            >  2:7089711/1‑150 (MQ=255)
                                                                                            gggtTCCGGCTGGGGCCGGTGTGGGGAACTTGTCCGCGGGAGCGGCGTAACGGGGCGAGCGGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaagag                                                           <  1:7304671/148‑4 (MQ=25)
                                                                                              gtTCCGGCTTGGGCGGCTTTGGGGAAATTGTCCGCGGGAGCGGCGTAACGGGGCGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGc                                                         <  2:1789893/150‑1 (MQ=255)
                                                                                               ttccGGCGTGGGCCGCTTTGGGGAATTGGTCGGCGGCCCGGGCGTCAGGGGGCGAGCTCGGGGGGGGGGGGGGGCAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                        <  1:7696988/147‑1 (MQ=255)
                                                                                               ttACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                        >  1:2724709/1‑150 (MQ=255)
                                                                                               gcccGGCTTGGGCCGGTTGGGGGAATTGGTCCGCGGGCGCGGCGTCCCGGGGCGCGCGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                        <  1:6585306/147‑1 (MQ=255)
                                                                                                 aCGGCTTGGGCGGCTTTGGGGAAAGTGCCCGCGGGAGAGGCGTACCGGGGCGCGCGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGt                                                      <  2:2368813/150‑1 (MQ=25)
                                                                                                   ggCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGTGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTa                                                    >  1:2619140/1‑150 (MQ=255)
                                                                                                     gtggggCCCCTTGGGGTTAATTGCCCGTGGGTGCGCCTTACGGGGGCGAGCCGGGGGGGTGGGGGTGGGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACa                                                  <  1:2779832/147‑1 (MQ=255)
                                                                                                     ctggggCCGCTTTGGGGCAATGGGCCGCGGGCGCGGCGTCCCGGGGCGCGAGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACa                                                  <  2:523529/147‑1 (MQ=21)
                                                                                                     cctGGGCCGCTGTGGGGAACTGGGCCGCGGGCGCGGCGTACCGGGGCGAGAGCGGGGGGTGGGGGGGGGTATACAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACa                                                  <  2:636449/148‑1 (MQ=255)
                                                                                                        gggCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAAAGGCCACTTAGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                               >  1:3905803/1‑150 (MQ=255)
                                                                                                        gggCCGCTTGGGCGAACTTGTCGGCGGGAGAGGCGTACCGGGGCGCGGGGGGGGGGTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                               <  2:346676/150‑1 (MQ=255)
                                                                                                          tccGCTCTGGGGAAATTGGCCGCGGGTGAGGCGTACCGGGGCGAGATCGGGGGGTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCtt                                             <  2:6459752/149‑1 (MQ=255)
                                                                                                            cGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAt                                           >  1:2276746/1‑150 (MQ=255)
                                                                                                              gtttGGGGAACTCGGCCGCGGGAGTGGCGTCCTGGGGCGCGATGGGGGGGTGGGGGTGGCCATAAAGAATGCGTAAAAATCCGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttg                                         <  1:5524889/149‑1 (MQ=255)
                                                                                                              cTTTGGGGAACTGGTCCGAGGGAGAGGCGTAACGGGGCGCGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttg                                         <  2:8384043/150‑1 (MQ=255)
                                                                                                                ttGGGGAACGTGGCCGCGGGAGAGGCGTACCGGGGCGAGCTCGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgtt                                       <  2:3172190/150‑1 (MQ=255)
                                                                                                                  ggggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATtgttgt                                     >  2:2728448/1‑150 (MQ=255)
                                                                                                                  ggggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGCGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATtgttgt                                     >  2:5342784/1‑150 (MQ=255)
                                                                                                                   gggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGc                                                                         <  2:7669511/113‑1 (MQ=255)
                                                                                                                    ggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTGTAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                   <  2:2709965/150‑1 (MQ=255)
                                                                                                                    ggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                   >  2:8697325/1‑150 (MQ=255)
                                                                                                                    ggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTgg                                   >  1:5956524/1‑150 (MQ=255)
                                                                                                                         ttGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgct                              >  1:2013623/1‑150 (MQ=255)
                                                                                                                          tggccGCGGGTGTGGCGTACCGGGGCGAGTTGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCGTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctg                             <  2:849515/147‑1 (MQ=255)
                                                                                                                          tGTCCGTGGGAGAGGCGTACAGGGGCGTGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctg                             <  2:5807056/150‑1 (MQ=255)
                                                                                                                            tCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctgct                           >  2:7561455/1‑150 (MQ=255)
                                                                                                                              cgtgggAGCGGCGTTTTGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCtgctgc                         <  2:1427597/147‑1 (MQ=255)
                                                                                                                              cGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGc                                                                                                                       <  2:5083065/56‑1 (MQ=255)
                                                                                                                               gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGCAAGGGCCGGTACAGGCTTATTGTTGTGGCGGCTGCTgcg                        >  1:8076837/1‑150 (MQ=255)
                                                                                                                               gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTgcg                        >  2:4510004/1‑150 (MQ=255)
                                                                                                                                  ggAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                        <  2:8821021/115‑1 (MQ=255)
                                                                                                                                         cggaaaGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTGGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGTACAGGCTTCTGGTTGTGGCGGCTGCTGCGCCGCGTTACg              >  2:5300943/4‑150 (MQ=39)
                                                                                                                                          ttaccggggCGCGAGGGGGGGGTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTGGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACgg             <  1:7115226/145‑1 (MQ=255)
                                                                                                                                          gTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTTTGGCGGCTGCTGCGCCGCGTTACgg             >  1:1580766/1‑150 (MQ=255)
                                                                                                                                             aaGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc          >  2:1817421/1‑150 (MQ=255)
                                                                                                                                             aaGGCGCGAGATCGGGGGGTGGGGGTGCGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc          <  2:7924167/150‑1 (MQ=255)
                                                                                                                                               ggggCGCGCGCGGGGGGGGGGGGGGGCAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca        <  1:6905658/150‑1 (MQ=255)
                                                                                                                                                  gCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCCGAAAAATCAGCGTGGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGCAAGGGCCGGTACAGGCTGATTGTTGTGGCGGCTGCTGCGCCGCGTTACGCCtccaggc     >  2:6895565/1‑148 (MQ=255)
                                                                                                                                                    gagaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca   >  2:9042128/1‑150 (MQ=255)
                                                                                                                                                    gagaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca   >  2:4649011/1‑150 (MQ=255)
                                                                                                                                                     agaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCtt                                                                       >  1:7497825/1‑81 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGCGCAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAG  >  NC_002947/4707153‑4707451

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: