Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,710,332 T→G 11.0% intergenic (‑75/+102) 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,710,3320TG11.0% 625.2 / 13.1 223intergenic (‑75/+102)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (121/69);  new base G (2/29);  total (124/99)
Fisher's exact test for biased strand distribution p-value = 7.66e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 233 total aligned reads displayed.

GCGCAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAG  >  NC_002947/4710183‑4710481
                                                                                                                                                      |                                                                                                                                                      
gcgcAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGACATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAt                                                                                                                                                        >  2:6974629/1‑150 (MQ=255)
gcgcAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAt                                                                                                                                                        <  1:7681183/150‑1 (MQ=255)
gcgcAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAt                                                                                                                                                        >  1:6093868/1‑150 (MQ=255)
gcgcAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAt                                                                                                                                                        >  1:6029777/1‑150 (MQ=255)
gcgcAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAt                                                                                                                                                        >  1:5828417/1‑150 (MQ=255)
    agcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggg                                                                                                                                                   <  2:8538404/150‑1 (MQ=255)
    agcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggg                                                                                                                                                   <  2:3577022/150‑1 (MQ=255)
    agcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggg                                                                                                                                                   <  2:3656572/150‑1 (MQ=255)
    agcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggg                                                                                                                                                   <  2:5891082/150‑1 (MQ=255)
     gcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCgggg                                                                                                                                                  <  2:5528754/150‑1 (MQ=255)
     gcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCgggg                                                                                                                                                  <  1:1813228/150‑1 (MQ=255)
     gcagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCgggg                                                                                                                                                  <  2:8450397/150‑1 (MQ=255)
      cagcGCCACCCGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggggg                                                                                                                                                 <  1:3646675/150‑1 (MQ=255)
      cagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggggg                                                                                                                                                 <  2:3409663/150‑1 (MQ=255)
      cagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggggg                                                                                                                                                 <  1:6234726/150‑1 (MQ=255)
      cagcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCggggg                                                                                                                                                 <  1:5132208/150‑1 (MQ=255)
       agcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCgggggg                                                                                                                                                <  1:4557676/150‑1 (MQ=255)
       agcGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCGGTGGGGAAATTGTCCGAGGGAGAGGCGTGAAGGGGCGAGATCgggggg                                                                                                                                                >  2:691632/1‑150 (MQ=255)
        gcgcCACCAGCAGATAAACCGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggt                                                                                                                                               <  1:1863407/150‑1 (MQ=255)
        gcgcCACCAGCAGATAAAACGCCAATACCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggt                                                                                                                                               <  1:7613875/150‑1 (MQ=255)
         cgcCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtg                                                                                                                                              <  1:4091450/150‑1 (MQ=255)
          gcCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtgg                                                                                                                                             <  2:3365811/150‑1 (MQ=255)
           ccaccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtggg                                                                                                                                            <  2:2160074/150‑1 (MQ=255)
           ccaccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtggg                                                                                                                                            >  1:4194196/1‑150 (MQ=255)
           ccaccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtggg                                                                                                                                            >  1:8597041/1‑150 (MQ=255)
            caccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGATCTTGTTCCTGTTTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTTTCCGAGGGAGAGTCGTAAAGGGGCGAGATCGgggggtgggg                                                                                                                                           <  1:8455914/150‑1 (MQ=255)
            caccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtgggg                                                                                                                                           <  2:788145/150‑1 (MQ=255)
            caccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtgggg                                                                                                                                           <  2:8117048/150‑1 (MQ=255)
             accaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtggggg                                                                                                                                          <  1:7132320/150‑1 (MQ=255)
              ccaGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGgggggtgggggt                                                                                                                                         <  1:2653707/150‑1 (MQ=255)
                   agATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGta                                                                                                                                    >  1:544270/1‑150 (MQ=255)
                    gATAACACGCAAATGCCACTGCCAGCCGCTTTGACAGCCGGTCGCTCTTGTTCTTGTGTGGTTCTCCCGGTGTTAGGGCTGGGGCCGGTTTGGGGAATTTGTCGGCGGGCGCGGGGTTCAGGGGGGTGTGGGGGGGGTGGGGGTGGGtat                                                                                                                                   <  2:2826121/150‑1 (MQ=255)
                    gATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGtat                                                                                                                                   >  2:7076276/1‑150 (MQ=255)
                     aTAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGtata                                                                                                                                  >  1:6158500/1‑150 (MQ=255)
                       aaaaCGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATATCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATaag                                                                                                                                >  2:1713355/1‑149 (MQ=255)
                       aaaaCGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCAGGTGCAGGGGGGGATCTCCCCGGGTTACGGCGGGGGCCGCTTGGGGGCACTTGTCCGAGGGAGAGGCGGAAAGGGGCGAGCTCGGGGGGTGGGGGTGGGGCTaaa                                                                                                                                >  2:4357208/1‑150 (MQ=255)
                        aaaGGCCTATTCCTCTGCCTGCCGCTGTGTCAGCCCGTCTTGCTTGTTCTTGTGTGGATCTCGGCGTGTTTCGGCTTGGGCGGCTTTGGGGAATTTGTCGGTGGGAGTGGGTTAATGGGGCGTGTTCGGGGGGTGGGGGTGGGTATAAAg                                                                                                                               <  2:714931/150‑1 (MQ=255)
                         aaCGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGa                                                                                                                              >  1:7233460/1‑150 (MQ=255)
                         aaCGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGa                                                                                                                              >  2:1239349/1‑150 (MQ=255)
                         aaCGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGa                                                                                                                              >  2:1120260/1‑150 (MQ=255)
                          aCGCCAATGCCAATGCCAGCTGCTGTTTCAGCCCGTCTTGCTTGTTCTTGTGTGGTTCTCCCGTTTTTAGGGGTGGGGCCGTTTTGGTGTTTTTGTCGGTGGTTGAGGCTTAATGGGGCTTGTTGGGGGGGTGGGGGTGGGTATAAAGaa                                                                                                                             <  2:5621550/150‑1 (MQ=255)
                           cGCCATTGCCACTGCCAGCTGCTGTGACCGCCGTTCGTGCTTTTTCTTGTGTGGATCGCCCGTTGTTCGGGGTTGGGCGGGTTTGGGTAATTTGTCGGTGGGTGAGGCGTAATGGGGCGTGATGGGGGGGTGGGGGTGGGTATAAAGAAt                                                                                                                            <  2:6282278/150‑1 (MQ=255)
                            gCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATg                                                                                                                           >  1:1766868/1‑150 (MQ=255)
                            gCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATg                                                                                                                           >  1:3489542/1‑150 (MQ=255)
                             ccAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGc                                                                                                                          >  2:8415588/1‑150 (MQ=255)
                               aaTGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGt                                                                                                                        >  1:142503/1‑150 (MQ=255)
                                aTGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTa                                                                                                                       >  2:7379483/1‑150 (MQ=255)
                                aTGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGGGGGGGGGGGTCTAAAGAATGCGTa                                                                                                                       >  2:8499968/1‑150 (MQ=255)
                                     aCTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGGGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAAt                                                                                                                  >  2:4546890/1‑150 (MQ=255)
                                       tGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGATTGGGCAGATTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCa                                                                                                                >  2:875218/1‑150 (MQ=255)
                                         ccAGCAGCAGTGACAGCCCGTCGAGCATGTTCAGGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGc                                                                                                              >  1:1292498/1‑150 (MQ=255)
                                           agcagcagTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGCGAGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCgg                                                                                                            >  2:3685189/1‑149 (MQ=255)
                                            gcttgagTGACTGCCGGTTGTGTTTGTTTTTGTTTGTATGTCCGGGTGTTAGGGCTTGGGCGTGTTTGGGTATTTTGTCGGTGGTTGTGGGGTTCTGGGGGGTGTTGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGtt                                                                                                           <  2:2058706/145‑1 (MQ=255)
                                               gcatTGACTGCCCGTCGTTATTTTTCTTGTGTGGATCTGCGGGTGTTTGGGCTTGGGCGGCTTTGGGTATTTTGTCGGTGGGTGTGGGGTATTGGGGCGCGGGGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTtgtg                                                                                                        <  2:2458778/150‑1 (MQ=255)
                                                 agTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGt                                                                                                      >  1:5179288/1‑150 (MQ=255)
                                                  gTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTc                                                                                                     >  2:7291189/1‑150 (MQ=255)
                                                          ccGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTa                                                                                             >  1:7787410/1‑150 (MQ=255)
                                                                       tcttgtgtgGATCGCCCGTTGTTCGGGCTGGGGCCGGGTGGGGGAATTGTCCCGCGGGAGTGGCGCAACGGGGCGAGAGGGGGGGGGGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTcgc                                                                                <  2:6727621/147‑1 (MQ=255)
                                                                       tCATGGGTGGACCTCCCCGTGTTACGGCGTGGGCCGCTTTGGGGAATTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGAGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTcgc                                                                                <  1:6895587/150‑1 (MQ=255)
                                                                          gggggggATCACCCCTTGTTCGGGCTTGGGCCGCTTTGGGGAATTTGGCCGTGGGAGCGGCGTAACGAGGCGAGAGCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCggg                                                                             <  2:5162876/145‑1 (MQ=255)
                                                                            tttgGTTCTCCCCGTGGTACGGCTTGGGCGGCGTTGGTGAACGTGGCCGAGGGCGCGGCGTAACGGGGCGAGCCGGGGGGGGGGGGGCGCGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCt                                                                           <  2:5788227/148‑1 (MQ=255)
                                                                              tgGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCACTTCGCGGGCTTg                                                                         >  1:5751159/1‑150 (MQ=255)
                                                                               gtaTCGCCCCTTTTTAGGGCTTGGGCCGCTTTGGGTAATTTGTCCGCGGGTGAGGCTTAATGGGGCGTTTTGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGc                                                                        <  2:3008793/148‑1 (MQ=255)
                                                                                gATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGcc                                                                       >  1:4915005/1‑150 (MQ=255)
                                                                                gATCGCCCCGTGTGAGGGCTGGGGCGGCTTGGGGGAATGGGCCCGCGGGAGCGGCGTCACGGGGCGCGCGCGGGGGGTGGGGGCGCCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGcc                                                                       <  1:6290439/150‑1 (MQ=255)
                                                                                 ttcgcccgGTGTGCCGGCTGGGGCGGGTTGGGGCAATGTGTCCGGGGGCGAGGCGTACGGGGGGGCGAGCGGGGGGGGGGGGGGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGccc                                                                      <  1:5969261/146‑1 (MQ=255)
                                                                                 aTCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGccc                                                                      >  2:7645695/1‑150 (MQ=255)
                                                                                    tcCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCcgcg                                                                   >  2:4383824/1‑150 (MQ=255)
                                                                                    gccctGTGGTCCGGGTTGGGCCGCTGTGGGTAAATTGTCCGCGGGAGCGGGTTTATGGGGCGCGATGGGGGGGGGGGGGGGGCCATAAAGAATGCGTAAATATCAGCGGTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCAAGCCcgcg                                                                   <  2:752186/149‑1 (MQ=255)
                                                                                     ccccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGa                                                                  >  2:5936646/1‑150 (MQ=255)
                                                                                      cccGTGTTCCGGCTTGGGGCGCTTGGGGGAATTTGTCCGCGGGCGCGGGGTACCGGGGCGAGCGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                 <  2:6047278/150‑1 (MQ=255)
                                                                                      cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                 >  2:4837904/1‑150 (MQ=255)
                                                                                      cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                 >  2:6511455/1‑150 (MQ=255)
                                                                                       ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                >  1:5413529/1‑150 (MQ=255)
                                                                                       ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                >  2:1535237/1‑150 (MQ=255)
                                                                                        cGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGGTTGCCCGCGaaga                                                               >  1:8274569/1‑148 (MQ=255)
                                                                                         gtgtTCCGGCTTGGGCCGCTTGGGGCAAATTGTCCGCGGCTGTGGCGTAACGGGGCGAGAGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaagag                                                              <  1:2484824/150‑4 (MQ=25)
                                                                                         gtgtTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAgg                                                              >  1:3067308/1‑150 (MQ=255)
                                                                                           gtTCCGGCTTGGGCCGCTTGCGGGAATTTGTCCGTGGGAGCGGCGTACCGGGGCGTGAGGGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGc                                                            <  2:6168779/150‑1 (MQ=255)
                                                                                             tACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCg                                                          >  1:6690296/1‑150 (MQ=255)
                                                                                              aCGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGt                                                         >  1:4716312/1‑150 (MQ=255)
                                                                                                     gggCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAATGGCCGTTACAGGc                                                  >  2:7353792/1‑150 (MQ=255)
                                                                                                     gggCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                  >  2:1374153/1‑150 (MQ=255)
                                                                                                        ccGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTa                                               >  2:4709310/1‑150 (MQ=255)
                                                                                                                  gaattggTCCGCGGGCGGGGCGTCCCGGGGGGCGATGGGGGGGGGGGGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggc                                     <  1:2265796/144‑1 (MQ=255)
                                                                                                                  gAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTGATTGTTGTggc                                     >  2:2333612/1‑150 (MQ=255)
                                                                                                                     aTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggcggc                                  >  2:7339869/1‑150 (MQ=255)
                                                                                                                            cGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTATTGTTGTGGCGGCggctgc                            >  1:15192/1‑150 (MQ=255)
                                                                                                                             gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTgcg                           >  1:8731616/1‑150 (MQ=255)
                                                                                                                             gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTgcg                           >  2:3957798/1‑150 (MQ=255)
                                                                                                                             gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTgcg                           >  1:2397894/1‑150 (MQ=255)
                                                                                                                                  agagGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCcgcg                      >  1:8229268/1‑150 (MQ=255)
                                                                                                                                   gagGCTTACAGGGGCGAGATGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGt                     <  2:2681026/150‑1 (MQ=255)
                                                                                                                                      gcgtcccggggCGCGCGCGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTAc                  <  2:8674623/143‑1 (MQ=255)
                                                                                                                                      ccgtaccggggCGCGCGCGGGGGGGGGGGGCGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTAc                  <  2:8674709/143‑1 (MQ=255)
                                                                                                                                        gtcccggggCGCGCTCGGGGGGGGGGGGGGGCCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACgg                <  1:5023303/145‑1 (MQ=255)
                                                                                                                                        gTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACgg                >  1:8066334/1‑150 (MQ=255)
                                                                                                                                        gTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCCAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCCGCTGCGCCGCGTTACgg                >  1:5323370/1‑150 (MQ=255)
                                                                                                                                        gTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTAGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACgg                >  1:1580356/1‑150 (MQ=255)
                                                                                                                                         tAAAGGGGCGAGATCGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTAGGGGTCTGCCTGTAACGGCCACTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTAATAGTAGTGGCGGCTGCAGCGCCGCGTTACgg                >  2:710282/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCGCAGCAGCGCCACCAGCAGATAAAACGCCAATGCCACTGCCAGCAGCAGTGACAGCCCGTCGAGCATGTTCATGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAG  >  NC_002947/4710183‑4710481

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: