Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	389,867	T→C	intergenic (‑140/‑384)	hemB ← / → yaiT	porphobilinogen synthase/putative autotransporter YaiT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	389,867	0	T	C	97.6%	151.7 / ‑6.3	43	intergenic (‑140/‑384)	hemB/yaiT	porphobilinogen synthase/putative autotransporter YaiT
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/27);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATG  >  NC_000913/389619‑390110                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    caaAGATCTGCCACACCAGGCCGTTAAGGCTAAGTGGTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGCGGGCGTTGACGTAAGTCTGACATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGGATAACTGAAGCGGGGTAGGGAGGTTTTCAGAAAGCTGTCGCGATGTTGAAAAGAAGAGACTGGAAGAGAGGCCAGGACATCTTGTTAATAGAACGAGCAATACttt                                                                                                                                                                                                                                                   >  2:187792/1‑251 (MQ=255)    aGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTtat                                                                                                                                                                                                                                                >  1:242988/1‑251 (MQ=255)     gATCGGCAACAACAGGTCGTTAAGGATAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGCGGGCGGTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGAGATGTTGACAAGCAGAGAATGGAAAAGAGGCCAGGACATCTTTTTAGTCGAAACAGCAATAACTTtat                                                                                                                                                                                                                                                >  2:158391/1‑250 (MQ=255)      aTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTatat                                                                                                                                                                                                                                              <  2:370246/251‑1 (MQ=255)       tCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACCTNTTTTTAATAGAAACAGCAATAACTTTATAtt                                                                                                                                                                                                                                             >  2:242177/1‑251 (MQ=255)           cAACACCAGGTCGTTACGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGCTGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAACGTTGTCGCGATGTTGACAAGAAGAGAATGGACGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACt                                                                                                                                                                                                                                         >  1:254912/1‑251 (MQ=255)             acacCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGa                                                                                                                                                                                                                                       <  1:2625/251‑1 (MQ=255)             acacCAGGTCGTTAAGGCTAAGTGTTGTCTCNTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGa                                                                                                                                                                                                                                       <  2:240601/251‑1 (MQ=255)               acCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGaaa                                                                                                                                                                                                                                     <  1:520264/251‑1 (MQ=255)                     tCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGGAGGGATGTTTTACGAAAGATGTCGCGATGTTGCCAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCACTAACTTTATATTAAATGAAACTAttt                                                                                                                                                                                                                               >  2:359307/1‑251 (MQ=255)                          aaGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTaa                                                                                                                                                                                                                           <  1:388197/250‑1 (MQ=255)                             gCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGCTGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtt                                                                                                                                                                                                                       >  1:222344/1‑251 (MQ=255)                                  gtgtTGTCTCTTCAAACAGAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttattt                                                                                                                                                                                                                  <  1:38584/251‑1 (MQ=255)                                    gttgtCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCtttattta                                                                                                                                                                                                                 <  1:168264/250‑1 (MQ=255)                                                 aCATAGCGCGCAGCGCAGGAGATTTGCGCAGGCAACGAGGGCGTTGGATTAAGGATGTCATGGTCTGCCTGATGTTTGTGTAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAATTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTATTCTTTATTTATCGCTGTTGGttat                                                                                                                                                                                                   >  1:455/1‑251 (MQ=255)                                                   aTAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGttatta                                                                                                                                                                                                 <  1:115356/251‑1 (MQ=255)                                                         ggcagcgcggGAGGTTTGCGCATGCGACGAGNGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGGTTGTGGAAGCGAGGGGCATAGTAGACCTGAAGCAGGGTAGGGATGTTTTACGAAATTTTTCGCGATGTTGAAAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAAGCTTTATTTATAGCTGTTGGTTATTAtttttt                                                                                                                                                                                           <  2:254912/250‑1 (MQ=255)                                                         cgcagcgcagGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTAtttttt                                                                                                                                                                                           <  1:181674/251‑1 (MQ=255)                                                                     gaTTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTAATGGTCTGCCTGACGTTTGTGGAATAGAGGGGCATAGTATACCTGATGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCAAGGCCATCTTTTTATTAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGttt                                                                                                                                                                               >  2:435099/1‑250 (MQ=255)                                                                                  gacgaGGGAGTTGGATTAAGTCGGTCATGGTCTGCATGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGAAAAGAAGAGAATGGAAGAGAGTCCAGGTCATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTGTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAAt                                                                                                                                                                  <  2:319961/251‑1 (MQ=255)                                                                                                  tAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACt                                                                                                                                                   <  1:217768/250‑1 (MQ=255)                                                                                                            aTGGTCTGCCTGATGTTTGTGGAATCGAGGGNCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCatat                                                                                                                                        <  2:240884/251‑1 (MQ=255)                                                                                                              ggTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATCTTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCgct                                                                                                                                                                           >  2:452549/1‑214 (MQ=255)                                                                                                              ggTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCgct                                                                                                                                                                           <  1:452549/214‑1 (MQ=255)                                                                                                                                ggAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTgg                                                                                                                     <  1:413160/250‑1 (MQ=255)                                                                                                                                          ggCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGgtttgtgtttgt                                                                                                          <  1:242865/251‑1 (MQ=255)                                                                                                                                                               ggTAGGGATGTTTTACGAAAGGTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCt                                                                                     <  2:42617/251‑1 (MQ=255)                                                                                                                                                                gcaGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATCACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATCGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATATTTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATCTTTATAATGCGTATGCtt                                                                                    >  2:72536/3‑251 (MQ=255)                                                                                                                                                                        gTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAAc                                                                             >  1:456408/1‑250 (MQ=255)                                                                                                                                                                          tttACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGCATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACaa                                                                           >  1:505765/1‑250 (MQ=255)                                                                                                                                                                                aaaGTTGTCTCGATGTTGACAAGAAGAGAATGGAAGAGAGGCAAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGATATTATTTTTTGGAGTTTGTTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAACATTATc                                                                    >  2:83214/1‑251 (MQ=255)                                                                                                                                                                                 aaGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTTAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCt                                                                   <  2:541356/251‑1 (MQ=255)                                                                                                                                                                                      gTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCaa                                                               <  1:444561/250‑1 (MQ=255)                                                                                                                                                                                        cgcgATGTTGACAAGAAGAGAATGGAAGAGNGGCCAGTACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACa                                                             <  2:236566/250‑1 (MQ=255)                                                                                                                                                                                                 gaaaagaagAGACTGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTATTGGTTATTATTTTTTGGCGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGCTATTTATCTGGTTTGTGTTTGTGATATTTATTATGAGTAGGCTTAAAAAACAAAATTATCTGCAACATGTTTAAAt                                                    >  2:357863/4‑250 (MQ=255)                                                                                                                                                                                                 gACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGAAACATGTTGAAAt                                                    <  1:374081/250‑1 (MQ=255)                                                                                                                                                                                                      gaagaGAATGGAAGAGAGGCCAGGACATCTTTGTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATGTTTAATCTGTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCTCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTGTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACtgtg                                              <  2:504622/251‑1 (MQ=255)                                                                                                                                                                                                        agagAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCt                                            <  1:1422/251‑1 (MQ=255)                                                                                                                                                                                                                            ggACATCTTTTTAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGAAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATTCTGTGCTTTTATGAATTTGATGCgtgt                        <  1:303665/251‑1 (MQ=255)                                                                                                                                                                                                                                    ttttAATAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTGGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTGTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGAGGCGTGttttt                    <  1:468040/247‑1 (MQ=255)                                                                                                                                                                                                                                         aTAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGa                                                                                                        >  1:484010/1‑158 (MQ=255)                                                                                                                                                                                                                                         aTAGAAACAGCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGa                                                                                                        <  2:484010/158‑1 (MQ=255)                                                                                                                                                                                                                                          tAGAAACAGCAATAAATTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGCAGTTTGGTTGCGCTGCTATAAATTGAATAATTTAAACTTTGTTGCATATTCATAGGACATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACACAATTATCTGCAACTTGTTGAAATACTGTGATTTTATGAATTTGATTCGTGTTTTTCTCCATAAAt          >  2:35756/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gCAATAACTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTTCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATg  <  1:499600/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    CAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAGGGCGTTGGATTAAGTCTGTCATGGTCTGCCTGATGTTTGTGGAATCGAGGGGCATAGTATACCTGAAGCAGGGTAGGGATGTTTTACGAAAGTTGTCGCGATGTTGACAAGAAGAGAATGGAAGAGAGGCCAGGACATCTTTTTAATAGAAACAGCAATAATTTTATATTCACTGAAAATATTTTTAATCTTTATTTATAGCTGTTGGTTATTATTTTTTGGAGTTTGGTTGCGCTGCTATAAATTGAATAATTAAAACTTTGTTGCATATTCATAGGATATTTATCTGGTTTGTGTTTGTGATATTTATTATGCGTATGCTTCAAAAACAAAATTATCTGCAACATGTTGAAATACTGTGCTTTTATGAATTTGATGCGTGTTTTTCTCCATAAATTATATATG  >  NC_000913/389619‑390110

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