Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,321 A→G intergenic (+289/‑313) cra → / → mraZ DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,3210AG88.2% 138.2 / 7.0 52intergenic (+289/‑313)cra/mraZDNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (24/21);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCC  >  NC_000913/89076‑89567
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
gTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACt                                                                                                                                                                                                                                                   <  1:818029/251‑1 (MQ=255)
                     ccatatCCTTAAAGTAAATAAGACTAGTCCTGACTATTATTGATATGTTCTTTTAAACCCGCCCGTTAAGTAACTCACCAGATGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCCCTGTTCTTTGCTTTTGCGCTAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTc                                                                                                                                                                                                                              <  2:226494/249‑1 (MQ=255)
                       atatCCTTCAAGTAAATAAGACTATTCCTGCCTATTATTGATAAGTTCTTTTAAACCCGACCGTTAATTAAATCACCAGCTGCAATTAACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGAGATGAATAACCGATTTAACAAGCACTTTTCTTTCCTTTTGCGCAAACCCGCTAGCATCAAGCGCCACACAGACGCAACATGGACTGTTAACCGGGAAAGatatttgct                                                                                                                                                                                                                            >  2:80093/1‑246 (MQ=255)
                         atccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCTCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTa                                                                                                                                                                                                                           <  1:690524/247‑1 (MQ=255)
                                     aaTAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTcgc                                                                                                                                                                                                              >  1:871222/1‑251 (MQ=255)
                                                tCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTACCAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACATACGCAACAGGGCCTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTcgcaggaacgga                                                                                                                                                                                                     >  2:873789/1‑240 (MQ=255)
                                                tCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTcgcg                                                                                                                                                                                                             <  1:873789/241‑1 (MQ=255)
                                                   tGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGac                                                                                                                                                                                                 >  1:714039/1‑250 (MQ=255)
                                                          tattGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTtat                                                                                                                                                                                          <  1:193337/250‑1 (MQ=255)
                                                           attGATAAATTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACTCTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTgtat                                                                                                                                                                                        <  2:397238/251‑1 (MQ=255)
                                                               ataagttcTTTTAAACCCGCCCGTTAAGTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAATATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGc                                                                                                                                                                                    <  2:359102/244‑1 (MQ=255)
                                                                taagttcTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCt                                                                                                                                                                                   >  1:763115/7‑251 (MQ=255)
                                                                tAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCCTTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTACCCAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGCCACTTTATATTTGCt                                                                                                                                                                                   >  2:594036/1‑251 (MQ=255)
                                                                 aagttcTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCtg                                                                                                                                                                                  >  1:385195/6‑251 (MQ=255)
                                                                   atgcgttgaaccccgcctggttatgaaCTTACCAGCTGAAATTCACAATAATTAAGTGATATCGACATCGCGTTTTTGCATTATTTTGTTAAATGCGGCGATGAATTTCCTATTTATCCAACTCTTTCCTGTGCTTTTGCGCAAACCCGCTGGCATCATGCGCCTCACAGACGTACCAATGACTGTTAACCGTGGAATATATGTCCGAAAAGGACGCGCGCGTCGCAATCTGTCACTTGATAGTTGCtgtg                                                                                                                                                                                <  2:751569/226‑1 (MQ=255)
                                                                                cgcccgTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGCCTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGATGTGGAAAATAGTGAGt                                                                                                                                                                   >  2:742098/1‑251 (MQ=255)
                                                                                      ttaattaaCTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCAtttt                                                                                                                                                             >  1:316127/1‑251 (MQ=255)
                                                                                       taattaaCTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCAGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCaaa                                                                                                                                                                                                      >  2:344282/1‑209 (MQ=255)
                                                                                       taattaaCTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCAGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCaaa                                                                                                                                                                                                      <  1:344282/209‑1 (MQ=255)
                                                                                       taattaaCTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGt                                                                                                                                                                       >  1:546951/1‑240 (MQ=255)
                                                                                       taattaaCCCACCAGCTGAAATTCACAATAATTAAGTGATATCGACATCGTGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTGTCTTTGCTTTTGCGCAATCCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAATATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGCCACTTTATATTTGCTGTGGAAAATAGt                                                                                                                                                                       <  2:546951/240‑1 (MQ=255)
                                                                                                     gCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTTCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATTTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACg                                                                                                                                              <  1:604355/251‑1 (MQ=255)
                                                                                                     gCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGATTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATAAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACg                                                                                                                                              >  1:818669/1‑251 (MQ=255)
                                                                                                        gAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATGTTGTTACATGCAGAGATTAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATg                                                                                                                                           <  2:774132/251‑1 (MQ=255)
                                                                                                              cacaATAATTAAGTGCTATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTCACACGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGAGAAGATATGTCATAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAATTAGTGAGTCATTTTAAAACGTAGATGACTATGAggatt                                                                                                                                     >  2:829848/1‑248 (MQ=255)
                                                                                                              cacaATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGAt                                                                                                                                     <  1:80093/251‑1 (MQ=255)
                                                                                                                       tAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCtt                                                                                                                             >  2:50818/1‑250 (MQ=255)
                                                                                                                            gATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGc                                                                                                                        <  1:521198/250‑1 (MQ=255)
                                                                                                                                  gACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTCACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAACATAGTGAGTCATTTTAACACGGAGATGACGATGAGGGATTTTTTCTTACAGCTATTCAt                                                                                                                 >  2:804025/1‑251 (MQ=255)
                                                                                                                                   aCAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCGTTTGCGCAAACCCGCTGGCATCAAGCTCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTc                                                                                                                   <  2:465036/248‑1 (MQ=255)
                                                                                                                                     aGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCAACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCGACTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGGGAGTCATTTTAAAACGGTTATGAAGATGAGGGATTTTTTCTTACCGCTATTCATAAc                                                                                                              >  2:36874/1‑251 (MQ=255)
                                                                                                                                      gcgcGATTTTGCAATATTTTGTTACATGCATAGATGAATTGCCGATTTAACAAGCACTTTTCGTTGCTGTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACg                                                                                                             <  2:249674/251‑1 (MQ=255)
                                                                                                                                                 cATTATTTTGTTACCTGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGAGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGTTTTTTTCTTACAGCTATTCATAACGTTAATTTGCtg                                                                                                  >  2:556613/1‑250 (MQ=255)
                                                                                                                                                                     agATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATa                                                                              <  1:661965/250‑1 (MQ=255)
                                                                                                                                                                                cGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTg                                                                    >  2:608709/1‑250 (MQ=255)
                                                                                                                                                                                  aTTTAACAAGCACTTTTCTTTGATTCTGCGCAAACCCGCTGGCATCAAGCGCCACACAGTCTCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTAACAGCTATACATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGata                                                                 <  2:639485/251‑1 (MQ=255)
                                                                                                                                                                                            cacTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTaa                                                       <  1:233832/251‑1 (MQ=255)
                                                                                                                                                                                                       gCTTTTGCGCAAACCCGCTGGCATCAATCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGttt                                            >  1:657176/1‑251 (MQ=255)
                                                                                                                                                                                                          tttGCGCAAACACGCTGGCATCAAGAGCCACTCAGACGTCACAAGGACTGTCAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGAGATGACGATGAGGGATTTTTTCTTACCGCTATTCATAACGTTACTTTGCTTCGCACGTTGGCCGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGtttttt                                         >  2:126034/1‑251 (MQ=255)
                                                                                                                                                                                                              cgcAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGAAATTAGCCGTAAACATCGGGTTTTTTACCt                                     >  2:553308/1‑251 (MQ=255)
                                                                                                                                                                                                                   aCCCGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGt                                 >  1:487235/1‑250 (MQ=255)
                                                                                                                                                                                                                     ccGCTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATg                              >  1:144617/1‑251 (MQ=255)
                                                                                                                                                                                                                        cTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAAttt                                                                                                      >  2:784404/1‑176 (MQ=255)
                                                                                                                                                                                                                        cTGGCATCAAGCGCCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAAttt                                                                                                      <  1:784404/176‑1 (MQ=255)
                                                                                                                                                                                                                                   cgcCACACAGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCt                <  1:183955/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gcCACACAGACGCAACAGGGACTGTTAACAGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCtt               >  1:257400/1‑251 (MQ=255)
                                                                                                                                                                                                                                          caGACGCAACAGGGACTGTTATCCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACGCGTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCTCACTTTGGACCTAAAATAAACAACGCGGATATTAGCCGTAAACATCGGGTTTTGTACCTCGGTATGCCTTTTGACTGGCTTGACAAg         <  2:818669/251‑1 (MQ=255)
                                                                                                                                                                                                                                          caGACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAg         <  1:22558/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gACGCAACAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGCGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCACGTGATGCCTTGTTCCTGGCATGACCAGCt       >  2:275766/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cgcaaCAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCttt     <  1:853206/248‑1 (MQ=255)
                                                                                                                                                                                                                                               gcaaCAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCtttt    >  1:541498/3‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaCAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTcc  >  1:635300/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      
GTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGGACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCC  >  NC_000913/89076‑89567

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: