Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,357,202	C→T	intergenic (‑92/+221)	sapA ← / ← ymjA	putative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,357,202	0	C	T	100.0%	150.5 / NA	44	intergenic (‑92/+221)	sapA/ymjA	putative periplasmic binding protein SapA/DUF2543 domain‑containing protein YmjA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/21);  total (23/21)

TTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTCA  >  NC_000913/1356959‑1357445                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    ttGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGATATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGtggt                                                                                                                                                                                                                                              >  2:30970/1‑251 (MQ=255)          ttAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTg                                                                                                                                                                                                                                     >  1:509939/1‑251 (MQ=255)           tAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTg                                                                                                                                                                                                                                     <  2:206567/250‑1 (MQ=255)                   ttGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTACTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAgggg                                                                                                                                                                                                                             >  2:649057/1‑248 (MQ=255)                    tGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAACGCGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTAATAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGCGAc                                                                                                                                                                                                                           >  2:578679/1‑251 (MQ=255)                       cgttcccGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTc                                                                                                                                                                                                                        <  1:556011/247‑1 (MQ=255)                                         aGACAAACCCGCTGTCGCGGATAGCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGTGCTGTGTGGTTGGTCGTTGAACAGGGGACGTCGCCATCTGTTCAAc                                                                                                                                                                                                          <  2:454259/247‑1 (MQ=255)                                         aGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                                                                                                                                          >  1:454259/1‑247 (MQ=255)                                             aaaaCCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAc                                                                                                                                                                                                  <  2:377325/251‑1 (MQ=255)                                                      gTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAttt                                                                                                                                                                                         <  1:448926/251‑1 (MQ=255)                                                        cgcgGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtc                                                                                                                                                                                        <  1:615317/250‑1 (MQ=255)                                                        cgcgGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACACACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtct                                                                                                                                                                                       >  1:546930/1‑251 (MQ=255)                                                            gATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctat                                                                                                                                                                                   >  2:214583/1‑251 (MQ=255)                                                             atatCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctat                                                                                                                                                                                   >  1:100462/1‑250 (MQ=255)                                                             atatCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctat                                                                                                                                                                                   <  2:100462/250‑1 (MQ=255)                                                                       gCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                                                                                                                                        <  2:691574/251‑1 (MQ=255)                                                                            ggAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGa                                                                                                                                                                   <  1:617513/251‑1 (MQ=255)                                                                             gagaTTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                                                                                                                                  <  1:611203/251‑1 (MQ=255)                                                                               gaTTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaaca                                                                                                                                                                <  1:30970/251‑1 (MQ=255)                                                                                aTTCAGGCGCGGCGATTGCCTGCCCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaaca                                                                                                                                                                >  2:504769/1‑250 (MQ=255)                                                                                 ttCAGACGCGGCGCTTGCGTGACCACTCACAAGTACAGCAATCACCAAAAGAGACGATAATACCTGGCGCATACTATTAAGGGATTTTATGTACAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTATGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCAAAGaacaac                                                                                                                                                              >  2:675482/1‑251 (MQ=255)                                                                                       cgcgGCGATTGCCTGACCACTCACAAGTCCAGCAATCAACAAAAGAGACGATAATACCTCGCGCCTAATATTAAGGGCTTTTATGTAACGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCg                                                                                                                                                        >  2:47340/1‑251 (MQ=255)                                                                                           gcgATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTTTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACATGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagga                                                                                                                                                    <  1:371140/251‑1 (MQ=255)                                                                                            cgATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgaggag                                                                                                                                                   <  1:588416/251‑1 (MQ=255)                                                                                                 gCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCagag                                                                                                                                              <  2:482607/251‑1 (MQ=255)                                                                                                                  tCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTgaga                                                                                                                              >  2:642222/1‑250 (MQ=255)                                                                                                                   ccAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACACACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGGGAGaa                                                                                                                             <  2:385221/250‑1 (MQ=255)                                                                                                                                   gagaCGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCt                                                                                                            >  1:450722/1‑251 (MQ=255)                                                                                                                                       cgctactaCCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGCGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAAATTAGCGAGGAGCAGAGCGGTACGATCGTGAGACACCCATTAAAGCACGCTGc                                                                                                          >  2:506874/7‑249 (MQ=255)                                                                                                                                           aataCCTGGCGCATAATATTAAGGGAGTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAATCACGCTGCAGTGCg                                                                                                    <  1:553343/251‑1 (MQ=255)                                                                                                                                                              tAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGTCAGTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAAAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGAt                                                                                  <  1:649057/250‑1 (MQ=255)                                                                                                                                                              tAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGAt                                                                                  <  1:232279/250‑1 (MQ=255)                                                                                                                                                                     ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTg                                                                          >  1:375023/1‑251 (MQ=255)                                                                                                                                                                     ttttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTAAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTg                                                                          >  2:573171/1‑251 (MQ=255)                                                                                                                                                                                          tACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACaa                                                        <  1:274221/248‑1 (MQ=255)                                                                                                                                                                                                  tttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCCTCTGTTCAACATTCGTACCCGTGCTTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGAcc                                             >  2:627893/1‑251 (MQ=255)                                                                                                                                                                                                           aTTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTg                                    >  2:531820/1‑251 (MQ=255)                                                                                                                                                                                                                ccATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGCTACTGACCGATCGTCTGCAAtt                               >  2:557820/1‑251 (MQ=255)                                                                                                                                                                                                                ccATTACCGTTTGTGTTCAGGGGACGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAAtt                               >  1:209976/1‑251 (MQ=255)                                                                                                                                                                                                                        gTTTGTGTTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAATTAAAtcatt                       >  2:613274/1‑251 (MQ=255)                                                                                                                                                                                                                            gtgtTCAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGCAAAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGATTTATCCGGACAACGATACTGACCGCTCGTCTGCAATTAAAtcattactc                   >  1:602220/1‑251 (MQ=255)                                                                                                                                                                                                                                 cAGGGGTCGTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTGAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTTCAGTGCGTGGCATTGTTAGCCAGATACGGCGTGAACGCTTTATCCAGACAACGATACTGCCCGATCGTCTGCAATTAAAtcattacgcatta               <  2:456039/250‑1 (MQ=255)                                                                                                                                                                                                                                         gTAATGAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTg       >  1:115416/1‑250 (MQ=255)                                                                                                                                                                                                                                              gAGGGCTGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCGCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACTATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTCa  <  1:506874/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    TTGGATGGGTTAAAGGTGTTGACTTGCCCGCTGACGCAATAGACAAAACCGCTGTCGCGGATATCAGCATGCGGGGGAGATTCAGGCGCGGCGATTGCCTGACCACTCACAAGTCCAGCAATCACCAAAAGAGACGATAATACCTGGCGCATAATATTAAGGGATTTTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAAGCACGCTGCAGTGCGTGTCATTGTTAGCCAGATGCGGCGTGAACGCTTTATCCGGACAACGATACTGACCGATCGTCTGCAATTAAATCATTACTCATTACCCCATTGATTCA  >  NC_000913/1356959‑1357445

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