Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,466,729	G→A	15.9%	pseudogene (1338/2513 nt)	ydbA →	pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,466,729	0	G	A	15.9%	89.2 / 4.9	36	pseudogene (1338/2513 nt)	ydbA	pseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (21/8);  new base A (2/4);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCAC  >  NC_000913/1466467‑1466992                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                         gcggcggTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATACAGGTTGACGGAGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TACACCGGTACGGAAATTGCCGGCAAAAACGGGAAGGTTATTCAGGAACGCGACCTGGaggc                                                                                                                                                                                                                                                >  2:312983/1‑286 (MQ=17)         cACGATATTGATATCACCGGCGACAGCGCAACTGTGGATAACAAGGGCACCATGCCCGTCACCGCTCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGTGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCTTTGCCGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAgcggcg                                                                                                                                                                                                                                        <  1:312983/289‑1 (MQ=255)               tttGAGATCACCGGCGAAAGCGAAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAGCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGAACCGGCACGCAGATTAACGGGGATGACGCCACGGAAAACAACAAGGGCAAAAGCAGCGTCGGAGGTAAGGAA‑GACACCGGTACGGAAATTACTGGCAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑Gt                                                                                                                                                                                                                                  <  2:73482/287‑1 (MQ=17)                            gCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAg                                                                                                                                                                                                                                             <  2:409884/265‑1 (MQ=255)                            gCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCAATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAg                                                                                                                                                                                                                                             >  1:409884/1‑265 (MQ=255)                                        cGGTGGAGAACAAGGGCACCAGGACCGTAACCGATCCGGAGTCCATCGGTATCAAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCa                                                                                                                                                                                                                                <  2:138542/265‑1 (MQ=255)                                                  cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGat                                                                                                                                                                                                                      >  1:138542/1‑265 (MQ=255)                                                     gggCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCTATCACCAATGGTGGCACCGGCACGCAGATTAACGGTTACGACGCCACGGCGAACAACAACGGAAAAACCACCGTTGACGGAAAGGAT‑TACACCGGTACGGAAATCGCGGGCACTAACGGGAATGTGATTCTGGATGGCGACCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATAACATCGGCGACAGCGCAACGGTTGac                                                                                                                                                                                           >  2:228555/1‑288 (MQ=17)                                                                   tCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTCAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGACGACGCCACCGCAAATAACAGCGGCAAAACCACCGTTGACGGCAAAGAT‑TCCACGGGTACCGAAATTGCTGGCAATAACGGGAAGGTAAATCAGGACGCCTATATACATGTAAGC‑GTGGCGAGGCACGGTATTGTTAAACCCGGCGAGAGCGCCGTGGTGTATAACAAaggaatcaa                                                                                                                                                                             >  2:336426/1‑283 (MQ=255)                                                                              gAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCAGCCTGACCGTCAcc                                                                                                                                                                   >  1:364916/1‑288 (MQ=255)                                                                                   cATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGACGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCCACGGTGGATAAAAAGGGCACCATGACCGTCACAGATCCg                                                                                                                                                             >  1:132227/1‑289 (MQ=255)                                                                                    aTCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTAATTCAGGACGGGGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAAAGGTGGATAACCACGGAACGATGACCGTCACCGATCCgg                                                                                                                                                            >  2:411166/1‑289 (MQ=255)                                                                                           tCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATCGCTGGCAATAACGGGAAGGTGATTCAGGACGGAGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACCGCGCAACGGTGGCTAACAAGGACATCATGACCGTCACAGAACCGGAGTcc                                                                                                                                                       >  1:219691/1‑287 (MQ=255)                                                                                                   gACGGCTACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCa                                                                                                                                                                                                                                <  2:393450/206‑1 (MQ=255)                                                                                                   gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGCTCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATGACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTAccc                                                                                                                                             >  1:213560/1‑289 (MQ=255)                                                                                                   gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGACGACGCCACGGCAAATAACAGCGGCAAAACCACTGTTGACGGAAAGGAT‑TCCACCGGTACCGAAATCGCTGGCAATAACGGGAAGGTGAATGTGGACGGCGACCTGAATGTTACCGGCGGTG‑AGCATGGTATTGAGAGCATTGGTGACAGCGGAACGGTGGACAACAAAGGAACTATGGCTGTTACCGATCCGCAgcctattggcttgc                                                                                                                                             >  2:202533/1‑276 (MQ=255)                                                                                                   gACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCAACGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCa                                                                                                                                                                                                                                >  1:393450/1‑206 (MQ=255)                                                                                                     cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGt                                                                                                                                                                                                                            >  1:357452/1‑208 (MQ=255)                                                                                                     cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGt                                                                                                                                                                                                                            <  2:357452/208‑1 (MQ=255)                                                                                                     cGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGAAAATTGCGGGCAATAACGGGAGGGTGATTCAGGACGGCGACCTGGATGTCAGCGGCGGCG‑CGCACGGTATTGAGAACACCGGCGACAGCGCCACGGTGGATAACAAGGGCACCATTACCTGCACCGATCCGGAGTCCATTGGGATcccg                                                                                                                                           >  2:129363/1‑287 (MQ=17)                                                                                                                             cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAg                                                                                                                                                                                     >  1:271351/1‑223 (MQ=255)                                                                                                                             cGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAg                                                                                                                                                                                     <  2:271351/223‑1 (MQ=255)                                                                                                                                         cgcCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTAAATAACGAAGGc                                                                                                       >  2:356491/1‑289 (MQ=255)                                                                                                                                                    aCGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACGAGGCCATCCTCAATACACAAGGCGCGCGCACTa                                                                                             >  1:387768/1‑288 (MQ=255)                                                                                                                                                                                  atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGCTCa                                                               >  1:322768/1‑288 (MQ=255)                                                                                                                                                                                  atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCa                                                               >  2:27951/1‑288 (MQ=255)                                                                                                                                                                                  atgaCGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGAGTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCa                                                               >  1:91957/1‑288 (MQ=255)                                                                                                                                                                                                  agacaacGGAAAAACCACCGTGTACGGC‑AGTATGTCCACCGGTACGGAAATTGCGTGCAATAACGGGAAGGTGAGTCAGGACGGCGATCTGCATTTCTGCGGCGTCG‑GTCGCGGTATTGCGAGCACCGGCAAAAGCGCAACGGCGGCAACCTAGGGTACCTTTAAGGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCGATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcc                                              <  2:132227/287‑1 (MQ=255)                                                                                                                                                                                                     gaagtataaggtcaaggTTGACGGGAAAGAT‑TCCACCGGTACGGAAATCACCGGTCATAAGGTGAAAGTGATTCAGGACGGTGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTAATATCACCGGTGACAGCGCAACGGTGGGTAACAAGGGCACCATGACCGTCACCGATCCTGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCa                                                <  2:111872/268‑1 (MQ=11)                                                                                                                                                                                                                    cGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGCTTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACCCAGATCAACGGTAACGACGCCCCCGCGAATAAAAGTGGaaaa                            >  2:61322/1‑289 (MQ=255)                                                                                                                                                                                                                       tatccgCAAGGAT‑TCAACCGGTACGGAAATAAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTGACGGTATTGATATCACCGGGGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAAc                          <  2:159854/283‑1 (MQ=14)                                                                                                                                                                                                                                ggAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGTTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCa                                                <  1:298417/257‑1 (MQ=255)                                                                                                                                                                                                                                ggAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCa                                                >  2:298417/1‑257 (MQ=255)                                                                                                                                                                                                                                   t‑tCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac                                                                             <  2:165953/225‑1 (MQ=11)                                                                                                                                                                                                                                   t‑tCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGcggcac                                                                             >  1:165953/1‑225 (MQ=11)                                                                                                                                                                                                                                               cggaaatcaacggtaaTAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCAc  >  2:180109/15‑289 (MQ=11)                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                         GCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAT‑TCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCG‑GTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCACCGCGAATAACAGTGGAAAAACCACTGTTGATGGAAAAGATTCCAC  >  NC_000913/1466467‑1466992

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