Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,796,473	C→T	I165I (ATC→ATT)	waaC →	ADP‑heptose:LPS heptosyltransferase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,796,473	0	C	T	100.0%	139.7 / NA	42	I165I (ATC→ATT)	waaC	ADP‑heptose:LPS heptosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/26);  total (16/26)

AGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTT  >  NC_000913/3796227‑3796719                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aGCAGAGAACTATGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcac                                                                                                                                                                                                                                                    <  2:415183/251‑1 (MQ=255)  gCAGAGAACTATGACGCAGTTAGCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGcaca                                                                                                                                                                                                                                                   <  2:376643/251‑1 (MQ=255)      agaACTATGACGCAGTTATCGATGCACAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGGAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAGCCCAGGGCGATTATGCTATTGCACAGCAt                                                                                                                                                                                                                                               <  2:416528/251‑1 (MQ=255)            aTGACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtttcct                                                                                                                                                                                                                                          <  2:287234/250‑4 (MQ=255)              gACGCAGTTATCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGGGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACTCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGa                                                                                                                                                                                                                                        <  1:187813/250‑1 (MQ=255)                       aTCGATGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGAAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGcc                                                                                                                                                                                                                              >  1:51684/1‑251 (MQ=255)                                    ggCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCAGGGCGTAAAGCATGGATTGG‑ACTGGCAACGCGTTTGCGACCCGTTAGCCTGCCTGTTTTACAATCGTAAGGATCATATTGCAAAACAGCGGCACGCCGTAGGACGCTCCAGCGAACTGTTTTCCAAAAGTCTGGGCTATATCAATCCGTAGACCCAGTGCGATTGGGCTATTTCTCAGCAGTTCCTGTGGAACCTGCGTACGGATGCTGGc                                                                                                                                                                                                                 <  2:193809/251‑1 (MQ=255)                                     gCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGGTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATGATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCg                                                                                                                                                                                                                <  1:414378/251‑1 (MQ=255)                                              aaGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCg                                                                                                                                                                                                       <  1:330714/251‑1 (MQ=255)                                                  gcggcgCTGATGACGCGTCTGGCGCATGGCGTACAGCATGGATTGG‑ACTGGCAAACCGCTCGCGCACCGTTAGCCAGCCTGTTTTACACTCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGACCTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCAGTAtt                                                                                                                                                                                                   >  2:359712/1‑251 (MQ=255)                                                   cggcgCTGGTGACGCGTCTGGCGCATGGCTGAAAGAATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTAttt                                                                                                                                                                                                  <  2:196702/251‑1 (MQ=255)                                                       gcTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATttcttc                                                                                                                                                                                              <  1:283845/251‑1 (MQ=255)                                                                   tCTGGCGCATGGCGAAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGAAAGCCTGGTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACACGCGAACTGATTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGATATTGCACAGCATTTCCTGCCGAACCTGCCTACAGAAGCTGGCGAATCTGCCGAATTTCTTCATGCGAAGAccc                                                                                                                                                                                  >  1:334565/1‑251 (MQ=255)                                                                      ggCGCATGGCGGAAAGCATGG‑ATGGGACTGGCAAACCGCTCGCGAACCGTTAGCCAGCATGTTTTACTATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACTGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATTCCGTATTTGGTTATGCGATGACCCg                                                                                                                                                                                 <  2:366772/249‑1 (MQ=255)                                                                         gcagggCGTAAAGCATGGAATGG‑ACGGGCAAACCTCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCATCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGtgatg                                                                                                                                                                            <  2:68302/247‑1 (MQ=255)                                                                                gTAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAac                                                                                                                                                                      >  1:291820/1‑250 (MQ=255)                                                                                 tAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtt                                                                                                                                                                                                                                              >  1:203011/1‑177 (MQ=255)                                                                                 tAAAGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCAtt                                                                                                                                                                                                                                              <  2:203011/177‑1 (MQ=255)                                                                                  aaaGCATGGATTGG‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACt                                                                                                                                                                   <  2:93760/251‑1 (MQ=255)                                                                                       atggattgg‑ACTGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACGCTGGCCg                                                                                                                                                              <  2:270425/245‑1 (MQ=255)                                                                                                   tGGCAAACCGCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGCCCCGTGATTATAAACACTGGCCGGAAGACCACTg                                                                                                                                                   >  2:417850/1‑251 (MQ=255)                                                                                                       aaaccccaagcgAACCGTTAGCCAGCCTGTTTTACAACCGTAAGCATCATATTGCAACACCGCAGCACGCCGTAGAACGCACCCGCGAACTGCTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCCGACGAACCTGACTTCCGATGCTGGCGAATATGACGTATTTCTTCATGCGACGACCCGTGATGATAAACAATTGCCGGAAAAaccattgggg                                                                                                                                               >  2:434817/10‑243 (MQ=255)                                                                                                            gCTCGCGAACCGTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCTAGAatt                                                                                                                                          <  2:69861/251‑1 (MQ=255)                                                                                                                             cagccTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACAGTGGCCGGAAGAGCACTGGCGAGAATTGATTGGTTTACTGGCTg                                                                                                                         <  1:359712/251‑1 (MQ=255)                                                                                                                                   gTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGACAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACTAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTTATGATAAACACTGGCCGGAATAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTc                                                                                                                     <  2:395261/249‑1 (MQ=255)                                                                                                                                   gTTTTAAAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCAAACGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCAGTATTTCTTCATGCGACGACCCGTGATGATAAACCCTGGCCGGACGAACAATGGAGAGAATTGATTGGTTTAATGGCTGCTTc                                                                                                                     >  2:416025/1‑249 (MQ=255)                                                                                                                                                aaGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAAGACGGATTa                                                                                                       <  1:98133/250‑1 (MQ=255)                                                                                                                                                 aGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTaaa                                                                                                     <  1:181613/251‑1 (MQ=255)                                                                                                                                                       atatTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCg                                                                                               <  2:318111/251‑1 (MQ=255)                                                                                                                                                                     cagcaCGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGcgccgc                                                                                  >  1:23469/1‑250 (MQ=255)                                                                                                                                                                      agcaCGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGAGTGGTTTACTGGCTGATGCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCAt                                                                                <  2:291820/251‑1 (MQ=255)                                                                                                                                                                          cGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAgg                                                                            >  2:329839/1‑251 (MQ=255)                                                                                                                                                                                aGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGCGCATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTggg                                                                                           >  1:174101/1‑230 (MQ=255)                                                                                                                                                                                aGAACGCACCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACCGATTAAACTTCCGTggg                                                                                           <  2:174101/230‑1 (MQ=255)                                                                                                                                                                                       aCCCGCGAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGa                                                                >  2:283797/1‑250 (MQ=255)                                                                                                                                                                                            cgAACTGTTTGCCAAAAGTCTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGAc                                                          <  1:329839/251‑1 (MQ=255)                                                                                                                                                                                                                   gCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCCTGAGGAAGAAc                                                                      >  2:5084/1‑216 (MQ=255)                                                                                                                                                                                                                   gCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAAc                                                                      <  1:5084/216‑1 (MQ=255)                                                                                                                                                                                                                      ataGCAAACCGCAAACCCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTg                                  <  2:90213/249‑1 (MQ=255)                                                                                                                                                                                                                               cgcaaaccCAGGGCGATTATGCTATTGCACAGCCTTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTTCTTATGTTGAAGTATTgtcg                       >  1:237938/1‑248 (MQ=255)                                                                                                                                                                                                                                  aaaccCAGGGCGATTATGCTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTTTGTTGGAGTATTGCCGa                      >  2:135470/1‑249 (MQ=255)                                                                                                                                                                                                                                                    cTATTGCACAGCATTTCCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAAATTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGtt  >  1:45828/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AGCAGAGAACTATGACGCAGTTATCGACGCTCAGGGGCTGGTAAAAAGCGCGGCGCTGGTGACGCGTCTGGCGCATGGCGTAAAGCATGGCATGG‑ACTGGCAAACCGCTCGCGAACCTTTAGCCAGCCTGTTTTACAATCGTAAGCATCATATTGCAAAACAGCAGCACGCCGTAGAACGCACCCGCGAACTGTTTGCCAAAAGTTTGGGCTATAGCAAACCGCAAACCCAGGGCGATTATGCTATCGCACAGCATTTTCTGACGAACCTGCCTACAGATGCTGGCGAATATGCCGTATTTCTTCATGCGACGACCCGTGATGATAAACACTGGCCGGAAGAACACTGGCGAGAATTGATTGGTTTACTGGCTGATTCAGGAATACGGATTAAACTTCCGTGGGGCGCGCCGCATGAGGAAGAACGGGCGAAACGACTGGCGGAAGGATTTGCTTATGTTGAAGTATTGCCGAAGATGAGTCTGGAAGGCGTT  >  NC_000913/3796227‑3796719

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: ‑