Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,707,392 A→G 13.7% intergenic (‑162/+15) 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,3920AG13.7% 403.1 / 43.6 154intergenic (‑162/+15)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (71/61);  new base G (14/7);  total (85/69)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 167 total aligned reads displayed.

CCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAG  >  NC_002947/4707243‑4707541
                                                                                                                                                     |                                                                                                                                                     
ccgggctACGGCTTGGGCCGCTTTGGTTACCTGGTCCGTGGGAGCGGCGTCTCGGGGGGCGCTGGGGGGGGGGGGGTGGCAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       <  2:1216191/144‑1 (MQ=255)
ccGTGTTACGGCTTGGGCGGGTTTGGGGAATTTGTCCGAGGGTGAGGGGTTATGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       <  2:91/150‑1 (MQ=255)
ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       >  1:3291691/1‑150 (MQ=255)
ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       >  2:4395612/1‑150 (MQ=255)
ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       >  1:4897856/1‑150 (MQ=255)
ccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       >  1:566302/1‑150 (MQ=255)
 cGTTTTCCGGCTGGGGCGTCTTGGGGGACTGTGTCCGCGGGTGCGGCGTACGGGGGGGCGTGGGGGGGGTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                                      <  2:4558508/150‑1 (MQ=255)
 cGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                                      >  1:4625151/1‑150 (MQ=255)
  gggtcccGGCTTGGGCCGCTTTGGGCAACTGGCCCGCGGGCGCGGCGCCCCGGGGCGAGCGCGGGGGGTGGGGGGGGCCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAgg                                                                                                                                                     <  1:389721/144‑1 (MQ=255)
  ggggtccGGCTGGGGCCGCTTGGGGGTATTTGTCCGCGGGCGAGGCGTAAGGGGGCGAGCTCGGGGGGTGGGGGTGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAgg                                                                                                                                                     <  2:4194344/144‑1 (MQ=255)
     ttcgggCTTGGGCGGGTTGGGGTAATGTGCCGGTGGGCGCGGCGTACCGGGGCGCCCTGGGGGGGGGGGGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                                                  <  1:709235/146‑1 (MQ=12)
     ttcgggCTGGGCCCGCTTTGGGCAACTTGCCCGCGGCCGCGGCGCACCGGGGCGCGCCGGGGGGGGGGGGGCGGCAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                                                                                                                  <  1:2173257/146‑1 (MQ=255)
     ttACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                                                                                                                  >  1:1297472/1‑150 (MQ=255)
     ttACGGCTCGGGCCGCTTTGGGGAATGTGTCCGCGGGCGCGGCTTCCGGGGGCGAGCGCGGGGGGGGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGcc                                                                                                                                                  <  2:242927/150‑1 (MQ=21)
          gCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTAc                                                                                                                                             >  1:5372802/1‑150 (MQ=255)
           ctggggCCGCTGTGGGCACCTGGCCCGCGGGCGCGGCGTCACGGGGCGCGCTCGGGGGGGGCGGGCGCCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACa                                                                                                                                            <  1:1928668/147‑1 (MQ=255)
                     tttGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgt                                                                                                                                  >  1:3259587/1‑150 (MQ=255)
                     cttGGGTAACTGGGCCGCGGGCGCGGCGCAACGGGGCGCGAGCGGGGGGGGCGGGCCGGTATAAAGAAGGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgt                                                                                                                                  <  1:3996248/149‑1 (MQ=255)
                        ggggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATtgttgt                                                                                                                               >  2:1006460/1‑150 (MQ=255)
                        ggggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATtgttgt                                                                                                                               >  1:4657178/1‑150 (MQ=255)
                         gggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGTCGTTACAGGCTTATTGTtgtg                                                                                                                              >  2:719449/1‑150 (MQ=255)
                             aaTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTGGCGGGCTTGGCCGCGGAAGGGCCGTTTGAGGGTTTTTGTTGTggggg                                                                                                                          >  1:250171/1‑150 (MQ=255)
                             aaTTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTggcgg                                                                                                                          >  2:3241302/1‑150 (MQ=255)
                                     gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAAAGGCCACATCGCGGGCTTTCCCCCGAAAGGGCCGTTACAGGCATATTTTTGTTGCGGCTGCTgcg                                                                                                                  >  1:1071483/1‑150 (MQ=255)
                                     gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGGATAAAGAATGCGTCAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTCTTGTTGTGGCGGCTGCTgcg                                                                                                                  >  2:1437474/1‑150 (MQ=255)
                                           gagGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATCCGTAAAAATCAGCGGTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCATGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCAGCGCCGCGt                                                                                                            >  2:2811537/1‑150 (MQ=255)
                                              gCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTTTGGCGGCTGCTGCGCCGCGTTAc                                                                                                         >  1:2972467/1‑150 (MQ=255)
                                                gTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTTACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTAAAGGCTTTTTGTTTTTGCGGCTGCTTGGCCGCGTTACgg                                                                                                       >  1:2101030/1‑150 (MQ=255)
                                                 tAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGc                                                                                                      >  2:4691537/1‑150 (MQ=255)
                                                 tAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTATTTTTGTGGCGGCTGCTGCGCCCCGTTACGGc                                                                                                      >  1:5030263/1‑150 (MQ=255)
                                                   ccggggCGCGCGCGGGGGGGGGGGGTGGCAACAAAGAAGGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc                                                                                                    <  1:2077353/148‑1 (MQ=255)
                                                    cggggCGTGCGCGGGGGGGGGGGGTGGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGACCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcc                                                                                                   <  1:4736490/149‑1 (MQ=255)
                                                     ggggCGCGCGCGGGGGGGGGGGGTGGGAATAAAGAAGGCGGAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtcca                                                                                                  <  2:726133/150‑1 (MQ=255)
                                                        gCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCGTGCCCGCGAAAGGGCCGTTACAGGCTGATTGGTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtc                                                                                               >  1:88097/1‑150 (MQ=255)
                                                         cgcccccGGGCGGCGGGGGTCGCTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCGGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcc                                                                                              <  1:1053105/144‑1 (MQ=255)
                                                         cGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACAGCtccagtcc                                                                                              >  2:1357984/1‑150 (MQ=255)
                                                          gcgcgcGGGGGGTGGGGGGGGGCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                             <  2:579123/145‑1 (MQ=255)
                                                          gagcgcGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                             <  2:3680294/145‑1 (MQ=255)
                                                          gagaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                             >  2:1625498/1‑150 (MQ=255)
                                                            gaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       >  1:271861/1‑90 (MQ=255)
                                                            gaTCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                                       <  2:271861/90‑1 (MQ=255)
                                                             ctgggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagca                                                                                          <  2:3935678/147‑1 (MQ=255)
                                                                ggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagcagca                                                                                       >  1:3816832/1‑150 (MQ=255)
                                                                    ggtgggggtgggTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCt                                                                                   >  2:4075098/1‑150 (MQ=255)
                                                                           gggggaATAAAGAATGCGTAAAAATCACCGTTGGGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCt                                                                            <  2:357094/148‑1 (MQ=255)
                                                                             gggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGACGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAc                                                                          <  2:4350664/150‑1 (MQ=255)
                                                                              ggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTAcc                                                                         <  2:3816832/150‑1 (MQ=255)
                                                                                 ataAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTTCCCGCGAAAGGGCCGTTAAAGGCGTAATGTTGTGGCGGCTGCGGCGCCGCGGTACGGCTCCAGTCCACCCGCAGGCGGTAGCCAACCGcc                                                                      >  1:1420930/1‑150 (MQ=255)
                                                                                  taAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagca                                                                                          <  2:333475/129‑1 (MQ=255)
                                                                                   aaaGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTATTGTTTTTGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCAACCCCCac                                                                    >  1:2550650/1‑149 (MQ=39)
                                                                                      gAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCaa                                                                 >  2:4794396/1‑150 (MQ=255)
                                                                                      gAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCaa                                                                 >  2:1419402/1‑150 (MQ=255)
                                                                                      gAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTAAAGGCTTATTGTTGTTGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCGGTAGCCTACCCCCAGcca                                                                 >  1:1993975/1‑148 (MQ=255)
                                                                                       aaTGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                             <  1:2648501/121‑1 (MQ=255)
                                                                                             tAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGc                                                          >  2:365110/1‑150 (MQ=255)
                                                                                                 aaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                                      <  2:635000/150‑1 (MQ=255)
                                                                                                  aTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAACAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                                     <  1:145687/150‑1 (MQ=255)
                                                                                                  aTCAGCGTGGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                                     <  2:4919819/150‑1 (MQ=255)
                                                                                                   tCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATg                                                    >  2:4656391/1‑150 (MQ=255)
                                                                                                    cAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGa                                                   <  1:1625498/150‑1 (MQ=255)
                                                                                                      gCGTTGTGGTCTGCCTGTAACGGCCTCTTAGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTTTGGCGGCTGCTGCGCCGCGTTAAGGCTCCAGTCCACCAGCAGGCTGTAGCCTAACGCCACCAAGGTAGGCCCCATGaaa                                                 >  1:3857071/1‑150 (MQ=255)
                                                                                                       cGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                                <  2:4859404/150‑1 (MQ=255)
                                                                                                       cGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                                <  2:318444/150‑1 (MQ=255)
                                                                                                          tgtgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                                             <  1:1439315/150‑1 (MQ=255)
                                                                                                            caataagCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGc                                                          <  2:2076558/129‑1 (MQ=25)
                                                                                                                  gCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAGGGCCAATGaa                                      >  2:1677148/1‑149 (MQ=38)
                                                                                                                    cTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                                   <  1:2844300/150‑1 (MQ=255)
                                                                                                                    cTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                                   >  1:3359690/1‑150 (MQ=255)
                                                                                                                     tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCt                                                                            >  1:4396554/1‑108 (MQ=255)
                                                                                                                     tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                                >  2:5376641/1‑136 (MQ=255)
                                                                                                                     tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                                <  1:5376641/136‑1 (MQ=255)
                                                                                                                     tGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCa                                  >  1:3179490/1‑150 (MQ=255)
                                                                                                                         aCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAgg                              >  1:1736027/1‑150 (MQ=255)
                                                                                                                         aCGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAg                                                                >  2:2193448/1‑116 (MQ=37)
                                                                                                                           ggccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                                     <  2:3008254/119‑1 (MQ=38)
                                                                                                                           ggccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgca                            >  1:1196236/7‑150 (MQ=38)
                                                                                                                           ggccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgca                            >  2:2053822/7‑150 (MQ=38)
                                                                                                                           ggCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgca                            >  1:2984009/1‑150 (MQ=255)
                                                                                                                            gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa                           >  1:768710/6‑150 (MQ=39)
                                                                                                                             ccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                             <  2:5402756/79‑1 (MQ=21)
                                                                                                                             ccTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcacg                          >  2:1881840/1‑148 (MQ=255)
                                                                                                                             ccTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaag                          >  2:2499237/1‑150 (MQ=255)
                                                                                                                              ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc                         >  1:4344186/1‑150 (MQ=255)
                                                                                                                              ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc                         >  1:4951097/1‑150 (MQ=255)
                                                                                                                              cctttCGCGGGCAAGCCCGCGAAGGGGCCGTTACAGGCTTATTGTTGTGTCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                                             <  1:4165099/79‑1 (MQ=25)
                                                                                                                              cctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc                         >  1:844691/4‑150 (MQ=39)
                                                                                                                              cctttCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaa                                                 >  2:1257470/4‑126 (MQ=255)
                                                                                                                              cctttCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaa                                                 <  1:1257470/123‑1 (MQ=255)
                                                                                                                               tctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagca                        >  1:4018916/1‑150 (MQ=255)
                                                                                                                                    gcgGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccg                   >  2:2569754/1‑150 (MQ=255)
                                                                                                                                      gggcaagCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccgcc                 >  1:2257089/7‑150 (MQ=255)
                                                                                                                                       ggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGccgccg                <  2:881965/150‑1 (MQ=255)
                                                                                                                                        gcaagCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGa               >  2:628029/5‑150 (MQ=255)
                                                                                                                                          ttGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAgg                                                                                     >  1:612224/1‑78 (MQ=255)
                                                                                                                                          aagCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAg                                     <  1:4713513/124‑1 (MQ=255)
                                                                                                                                             cccGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACg          >  2:661197/1‑150 (MQ=255)
                                                                                                                                               cgcgAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGcc        >  2:275304/1‑150 (MQ=255)
                                                                                                                                                  gAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCcag     <  2:3291691/150‑1 (MQ=255)
                                                                                                                                                   aaaGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCcagc    >  2:4643454/1‑150 (MQ=255)
                                                                                                                                                    aaGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCcagca   <  1:628029/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGCCGAACACGCCCAGCAG  >  NC_002947/4707243‑4707541

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