Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 725,370 G→A T556I (ACC→ATC)  kdpB ← K(+) transporting P‑type ATPase subunit KdpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913725,3700GA98.1% 163.3 / NA 53T556I (ACC→ATC) kdpBK(+) transporting P‑type ATPase subunit KdpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/22);  minor base C (1/0);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCAGTCAG  >  NC_000913/725129‑725611
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
tAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGtt                                                                                                                                                                                                                                            >  2:281289/1‑249 (MQ=255)
tAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGtt                                                                                                                                                                                                                                            >  1:41778/1‑249 (MQ=255)
    gCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTc                                                                                                                                                                                                                                      <  1:84810/251‑1 (MQ=255)
       aGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCgag                                                                                                                                                                                                                                   <  1:129041/251‑1 (MQ=255)
       aGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACAGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCAAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCg                                                                                                                                                                                                                                     >  1:166627/1‑249 (MQ=255)
                aaaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCat                                                                                                                                                                                                                          >  2:297892/1‑251 (MQ=255)
                 aaaaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCata                                                                                                                                                                                                                         <  1:238490/251‑1 (MQ=255)
                   aaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATAtt                                                                                                                                                                                                                       >  2:107816/1‑251 (MQ=255)
                   aaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATAtt                                                                                                                                                                                                                       <  1:352183/251‑1 (MQ=255)
                   aaGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAAGCACATGATGTTCAGCGCATTTAACTGCGGATCCGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCCATGCTGCAGGTGGTCAGCGAGCCACGGGTCATCCGCGTCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATAtt                                                                                                                                                                                                                       >  2:147125/1‑251 (MQ=255)
                    aGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTg                                                                                                                                                                                                                      >  2:344792/1‑251 (MQ=255)
                          tAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGcc                                                                                                                                                                                                                >  2:84997/1‑251 (MQ=255)
                              cAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGCCGGGCGAATGACGGCCCATGATGTTCAGCGCATTTAACTGCGGCTACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCtt                                                                                                                                                                                                            >  1:189308/1‑251 (MQ=255)
                               aaGGCTTTGAAAATCACCGCCCGGAGGATTGGGGGGTCGGGCAAATGCAGGCACATGATGTTCAGCGCATTTAACGGCGGATACGTGGCCGCGAATGCCGCCGGAATAAGGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCCCCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATGGCCCGCCTCttt                                                                                                                                                                                                           <  2:23681/251‑1 (MQ=255)
                                 ggCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCGGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTcg                                                                                                                                                                                                         <  2:373496/251‑1 (MQ=255)
                                                          aTTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGc                                                                                                                                                                                >  1:150571/1‑251 (MQ=255)
                                                                                 gCACATGATGTTCAGCGCATTAAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCGGAAGGTGGTCAGCGGGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCCCCTCGATCATCTTTATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCa                                                                                                                                                         <  1:240987/251‑1 (MQ=255)
                                                                                     atgatgTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAgcgc                                                                                                                                                     >  1:204501/1‑251 (MQ=255)
                                                                                              aGCGCATTTAACTGCGGATACGTTGCCGCGAAGGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGACCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATc                                                                                                                                            <  1:257420/251‑1 (MQ=255)
                                                                                                 gcATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCAAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGtt                                                                                                                                         <  1:138699/251‑1 (MQ=255)
                                                                                                                   gTTGCCGCTAATCCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAAGGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCCCCTCTATCAACTTGATCGGTTTAGAGTCGAGATCGACCATATTTCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCAt                                                                                                                       <  1:59877/251‑1 (MQ=255)
                                                                                                                           gAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTAcc                                                                                                                >  1:277849/1‑250 (MQ=255)
                                                                                                                             aTGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCaa                                                                                                              >  1:94087/1‑250 (MQ=255)
                                                                                                                                         taatGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGcc                                                                                                 >  1:44521/1‑251 (MQ=255)
                                                                                                                                          aatGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCt                                                                                                >  2:151913/1‑251 (MQ=255)
                                                                                                                                           atGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCCCCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTg                                                                                               <  1:129718/251‑1 (MQ=255)
                                                                                                                                             ggCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGCTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGCGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGAt                                                                                             >  2:329303/1‑251 (MQ=255)
                                                                                                                                                   gTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCGTCCGCCTGATACTGAc                                                                                       <  1:327422/251‑1 (MQ=255)
                                                                                                                                                      tttCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACgg                                                                                     >  1:129310/1‑250 (MQ=255)
                                                                                                                                                        tcgcccaaTCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTGGCCAATGTGCACCTCCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGAt                                                                                   <  2:271687/243‑1 (MQ=255)
                                                                                                                                                         cGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATc                                                                                  >  2:366184/1‑250 (MQ=255)
                                                                                                                                                                         aTGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCCGCTTGGCCTc                                                                 >  1:239630/1‑251 (MQ=255)
                                                                                                                                                                             tGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGgt                                                             >  1:79593/1‑251 (MQ=255)
                                                                                                                                                                              gAAGGGGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGgtg                                                            <  2:133017/251‑1 (MQ=255)
                                                                                                                                                                                     gtCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTc                                                     <  1:344792/251‑1 (MQ=255)
                                                                                                                                                                                          cGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCg                                                 >  2:123460/1‑250 (MQ=255)
                                                                                                                                                                                                   ggTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAatcat                                       >  2:307780/1‑251 (MQ=255)
                                                                                                                                                                                                    gTCATCGGCAGCTGTGTGCCAGTGGGCGCCACCTCGAGCAGCGTGATCGGGTTAGAGGCAAGGTCGCCCTGAGTGCTCGCCTCTTTCGCCACCTGGGTGCCGGAGTTCATCGCTACCGCGACAGCGGGCTGCGCCAGCGCCGGGGCAGCGTTTGTGCCGTCGCCGGTCATCTCTTCCAAACGACCTTCCTCCGGATACTGACGGATCAAGGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAatcatc                                      <  1:254167/251‑1 (MQ=255)
                                                                                                                                                                                                                  tttGCCAATGTGCTCCCCCTCGATCACCTTTATCGGGTTAGAGTCGAGATCGCCCTTATTCCCCTCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGTATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTcc                         <  1:163177/250‑1 (MQ=255)
                                                                                                                                                                                                                      ccAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCgcagc                    <  1:185841/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccAATGTGCACCACCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCgcagc                    >  1:263800/1‑251 (MQ=255)
                                                                                                                                                                                                                       cAATGTGCACCCCCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGCTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGCCACCCGCTTCCgcagc                    >  2:317253/1‑250 (MQ=255)
                                                                                                                                                                                                                       cAATGTGCACCACATCGATCAACTTGATCGGGATAGAGTCGAGATCGACCATATTGACCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGGGCCGTCGCCGGGCATCGCTACCAAACGACCTTCCGCAAGATACGGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGATACCCGCTTCCgcagcg                   >  2:17114/1‑250 (MQ=255)
                                                                                                                                                                                                                       cAATGTGCAACACCTCGATCAACTTGATCGGGCTAGAGTCGAGATCGACCATATTGACCGCCTCTTTCGACGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGACGGCCATCGCTACCAAACGACCTTCCGCCTGCTACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCgcagc                    >  2:17104/1‑250 (MQ=255)
                                                                                                                                                                                                                            tgCACCACCTCGATCAACTTGATCGGGTTAGAGTCGCGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCCCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGACGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCAGCTTCCGCCGCAATcgc              >  1:232879/1‑251 (MQ=255)
                                                                                                                                                                                                                              caccacCTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCgcgg            >  1:176275/1‑251 (MQ=255)
                                                                                                                                                                                                                               accaccTCGATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCgcggc           <  1:358081/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccTCGATCAACTGGATCGGGTTAGAGGCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGTCCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCa       <  1:97507/251‑1 (MQ=255)
                                                                                                                                                                                                                                     ccgtccAACTTGATCGGGTGAAAGCCGAGATCGACCATATCGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCCCCGCGACAGCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCAGt     <  2:321892/246‑1 (MQ=255)
                                                                                                                                                                                                                                       gATCAACTTGATCGGGTTCGAGTAGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTCCATCGCCACCGCGACATCTGCCTGCGCCAGAGCCGGAGCATCGTTGGTGCCGTCGACGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCAGTc  >  2:32678/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCAGTc    >  1:352375/1‑250 (MQ=255)
                                                                                                                                                                                                                                       gATCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCAGTc    >  1:311595/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aTCAACTTGATCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCAGTCAg  <  1:234543/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
TAACGCCAGGGGAATCAAAAAGACGATAATCAAGGCGTTGAAAATCACCGCACTGAGGATTGCGGAGTCGGGCGAATGCAGGCACATGATGTTCAGCGCATTTAACTGCGGATACGTTGCCGCGAATGCCGCCGGAATAATGGCGAAGTATTTCGCCACATCGTTGGCAATGCTGAAGGTGGTCAGCGAGCCACGGGTCATCAGCATCTGTTTGCCAATGTGCACCACCTCGATCAACTTGGTCGGGTTAGAGTCGAGATCGACCATATTGCCCGCCTCTTTCGCCGCCTGGGTGCCGGAGTTCATCGCCACCGCGACATCTGCCTGCGCCAGCGCCGGAGCATCGTTGGTGCCGTCGCCGGTCATCGCTACCAAACGACCTTCCGCCTGATACTGACGGATCAATGCCAGCTTGGCCTCCGGTGTCGCTTCGGCGAGAAAATCATCGACACCCGCTTCCGCAGCAATCGCGGCGGCAGTCAG  >  NC_000913/725129‑725611

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