Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,189,980 A→G T156T (ACT→ACC)  phoP ← DNA‑binding transcriptional dual regulator PhoP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,189,9800AG98.0% 159.8 / NA 49T156T (ACT→ACC) phoPDNA‑binding transcriptional dual regulator PhoP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (23/25);  minor base T (0/1);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATC  >  NC_000913/1189735‑1190226
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
ccGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGgtgt                                                                                                                                                                                                                                                    >  1:758874/1‑250 (MQ=255)
 cGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTAt                                                                                                                                                                                                                                                  <  2:299278/251‑1 (MQ=255)
                aaaaaaaGACGCAGTAATTTTGTCATCAGCGCAATTCGAACAGATAGCCCTGGCTGCGAACGGTGGTAATCTCTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGTCTTGCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGACTCTTTGCTGTCCACTTTGCCATTATGGCGTATCAACGTTTCCATAATGGTGTACTCGAACGCGGTCAGt                                                                                                                                                                                                                                   <  2:806082/251‑1 (MQ=255)
                  aaaaaGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGttt                                                                                                                                                                                                                                 <  1:747194/251‑1 (MQ=255)
                    aaaGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGa                                                                                                                                                                                                                               <  1:651148/251‑1 (MQ=255)
                     aaGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCACCGTTTCCATAACGGTGTATTCGAACGCGGTAAGTTTGAt                                                                                                                                                                                                                              >  2:589974/1‑251 (MQ=255)
                        aCGCAGTCATTTTTTCATCAGCGCAATTCGCACAGATAGCCCTGGCCGCGAACGGTGGTACTCACTTCTTGGGGATATTGCGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCAGCATCCGGCTAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATACACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCAc                                                                                                                                                                                                                           >  2:161339/1‑251 (MQ=255)
                            aGTAACTTTTACATCAGCGCATTTCTAACAGATAGCCGTTACCGCGAACGTTGGTAATAACTTCTTGTGGATATTGTGCCTGAATTTGTTTGCGCAGACGTCCCAGCAGTACATCAATGGTTTGGCTTTCCCGCTGCGCCGCATCCGGATAGAGTTGGATCATTACCGAATCTTTGCTGTCCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGGCAGGTGGATCACTTCg                                                                                                                                                                                                                       <  2:59811/251‑1 (MQ=255)
                             gTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATTGTGGATTCGAACGCGGTCAGTTTGATCACTTCGt                                                                                                                                                                                                                      <  1:636620/251‑1 (MQ=255)
                               aaTTTTTTCATCAGCGCAATTCGAACAGATATCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTc                                                                                                                                                                                                                     >  1:856213/1‑250 (MQ=255)
                                aTTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACAGCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCAt                                                                                                                                                                                                                   <  1:101413/251‑1 (MQ=255)
                                                               cccTGTCCGCGAACGGTGGTAATCACTTCTTTGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAAc                                                                                                                                                                                    >  1:608152/1‑251 (MQ=255)
                                                                         gAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGACATTATTGCGTATCAACGTTTCCATAATGGTGAATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCCACCTGCAACggg                                                                                                                                                                          >  1:141732/1‑251 (MQ=255)
                                                                          aaCGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACgggg                                                                                                                                                                         <  1:481495/251‑1 (MQ=255)
                                                                                       aCTTCTTGTGGATATTGTGCCTGAATTTTTTTGCGGAGACGTCTCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCACCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGGCCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGTTGTATTCGACCGCGGt                                                                                                                                                                                                                                       >  2:62601/1‑176 (MQ=255)
                                                                                       aCTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGAATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGAAGCGAAATGAc                                                                                                                                                            >  1:820872/1‑251 (MQ=255)
                                                                                       aCTTCTTGGGGATATTGTGCCTGAATGTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCGCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGt                                                                                                                                                                                                                                       <  1:62601/176‑1 (MQ=255)
                                                                                            ttGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCtgtg                                                                                                                                                       >  1:439557/1‑251 (MQ=255)
                                                                                                  atatTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCa                                                                                                                                                 <  1:161339/251‑1 (MQ=255)
                                                                                                         gCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGGCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCt                                                                                                                                          >  1:143643/1‑251 (MQ=255)
                                                                                                                 ttttttGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTAcgc                                                                                                                                   <  1:257110/250‑1 (MQ=255)
                                                                                                                  tttttGAGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCAGCAGCTCCGCCTCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTCATAGATAATTCACGGCGAGCGAGATCAACATGAAACGGGGGGAGCGAAATGACATGTGCCGCCAGACCGCTATTAcgccg                                                                                                                                 >  2:816686/1‑251 (MQ=255)
                                                                                                                     ttGCGCAGACGTCACATCAGTACAACTATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAAAGGTGTATTCGAACGCGGTCAGTTTGATCACTGCGTCATTAATAGATAAATCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCa                                                                                                                               <  2:591241/250‑1 (MQ=255)
                                                                                                                                  ccATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTTACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCAtt                                                                                                                 >  1:856437/1‑251 (MQ=255)
                                                                                                                                  ccATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCAtt                                                                                                                 <  1:571201/251‑1 (MQ=255)
                                                                                                                                        gTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCTAGAGAGATCAACCTGAAACGGGGGGCGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCTTTAATGCCTGCATTCGCGcc                                                                                                           >  1:32846/1‑251 (MQ=255)
                                                                                                                                              cAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCAcc                                                                                                     >  2:206814/1‑251 (MQ=255)
                                                                                                                                                  ggTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCtt                                                                                                 <  1:458981/251‑1 (MQ=255)
                                                                                                                                                     aTGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCaa                                                                                              <  1:596410/251‑1 (MQ=255)
                                                                                                                                                     aTGGCTTGCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGGCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGTAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATGAATGCCTGCATTCGCGCCATCACCTCTTCa                                                                                               <  2:112966/250‑1 (MQ=255)
                                                                                                                                                        gcgtttccGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTAGCTGACCACATTGCCATTATTGCGTATCAACGTTTCAATAAAGGTGTATTCGAACGCGGTCAGGTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTTAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTCCAatat                                                                                           <  2:613689/245‑1 (MQ=255)
                                                                                                                                                           ttctcGCAGCTCCGCATCCGGATAGATTTGGAGCATTAACGAATCTTTGCTGACCTCTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTGTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCCCCTCTTCAATATg                                                                                          <  2:75132/245‑1 (MQ=255)
                                                                                                                                                            tcacGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCCCTTCAATATGaa                                                                                        <  2:480546/247‑1 (MQ=255)
                                                                                                                                                               cGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTGGCTGACCACTTTGCCATTATTGCGTATCACCGTTGCCATAATGGTGTATTCGAACTCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGTGATCAACCTGAAACGGGGGTTGCGAAATTACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGTGCCATCACCTCTTCAATATGAAACgg                                                                                    <  1:572415/251‑1 (MQ=255)
                                                                                                                                                                gCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGCCCACTTTGCAATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCCTTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGCAATGACCTGTGAAGCCAGACCGCTATTCCGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATTTGACACGGt                                                                                   >  2:369420/1‑251 (MQ=255)
                                                                                                                                                                   gCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGttt                                                                                 >  1:601066/1‑250 (MQ=255)
                                                                                                                                                                       cGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGCCCGCTATTACGCCGCATTAATGCCTGCCTTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTc                                                                             >  2:689314/1‑250 (MQ=255)
                                                                                                                                                                        gCATCCGGATAGAGTTGGAGCATTAACGAATATTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACAGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACATCTTCAATATGAAACGGTTTAGTcac                                                                           >  2:501198/1‑251 (MQ=255)
                                                                                                                                                                           tCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACAt                                                                         >  1:514285/1‑250 (MQ=255)
                                                                                                                                                                              ggATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATactc                                                                     <  2:450823/251‑4 (MQ=255)
                                                                                                                                                                                        ggAGCATTAACGAATCTTTGCTGACCACGTTGCCATTATTGCGTATCAACGTTGCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACGTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCg                                                           <  2:299155/251‑1 (MQ=255)
                                                                                                                                                                                              ttAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACtt                                                     >  1:59815/1‑251 (MQ=255)
                                                                                                                                                                                                aaCGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTaa                                                   <  1:589974/251‑1 (MQ=255)
                                                                                                                                                                                                              gACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTg                                      <  2:28185/250‑1 (MQ=255)
                                                                                                                                                                                                                ccACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGACCGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTcc                                   <  1:653808/251‑1 (MQ=255)
                                                                                                                                                                                                                 cACTTTGCCATTATTGCGTATCAACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCt                                  >  1:550767/1‑251 (MQ=255)
                                                                                                                                                                                                                             attGCGTATCAACGTTTCCATAATGGTGTATGCGAACGCGGTCAGTTTGATCACGTCGTCATTAATAGATAAATCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCCTCAATATGAATCGGTTTAGCCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTCTCAc                      <  2:608152/251‑1 (MQ=255)
                                                                                                                                                                                                                              ttGCGTATCCACGTTTCCATAATGGTGTATTCGAACGCGGTCAGTATGATAACTTCGTAAAAACTAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTACACg                     >  1:690519/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ttCCATAATGGTGAATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATAcc        >  1:546865/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 taatgGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATc  >  1:563101/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
CCGCAGCGAGAGCGGGAAAAAAAGACGCAGTAATTTTTTCATCAGCGCAATTCGAACAGATAGCCCTGGCCGCGAACGGTGGTAATCACTTCTTGGGGATATTGTGCCTGAATTTTTTTGCGCAGACGTCCCATCAGTACATCAATGGTATGGCTTTCCCGCAGCTCCGCATCCGGATAGAGTTGGAGCATTAACGAATCTTTGCTGACCACTTTGCCATTATTGCGTATCAACGTTTCCATAATAGTGTATTCGAACGCGGTCAGTTTGATCACTTCGTCATTAATAGATAATTCACGGCGAGAGAGATCAACCTGAAACGGGGGGAGCGAAATGACCTGTGAAGCCAGACCGCTATTACGCCGCATTAATGCCTGCATTCGCGCCATCACCTCTTCAATATGAAACGGTTTAGTCACATAATCATCAGCACCGGCACTTAATACTTCGACTTTGTCCTGCCAGCTTTCACGGGCGGTTAATACCAGAATC  >  NC_000913/1189735‑1190226

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