Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,808,866	T→C	D109D (GAT→GAC)	waaA →	KDO transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,808,866	0	T	C	100.0%	178.0 / NA	54	D109D (GAT→GAC)	waaA	KDO transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/27);  total (27/27)

AGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTT  >  NC_000913/3808619‑3809112                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       aGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTg                                                                                                                                                                                                                                                     >  2:300529/1‑251 (MQ=255)       cGGCCTATCGAAAACGCTGGGGGGAACGTTACGGTTTTGACCGCCACCCGCTAAAACCAGGCGGCATTATGCTGCTCTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGAGTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCTTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTAGCTGCCGTATGACCTGCCCGAt                                                                                                                                                                                                                                               <  2:438018/251‑1 (MQ=255)                          gggTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTg                                                                                                                                                                                                                             >  1:85052/1‑250 (MQ=255)                            gTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAAt                                                                                                                                                                                                                          <  2:61138/251‑1 (MQ=255)                                  gTTACGGTTTTTACCGCCATCCGCTTAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTc                                                                                                                                                                                                                    >  2:182821/1‑251 (MQ=255)                                      cGGTTTTTACCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGGTCAGCACGTTTATCTTCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGAcc                                                                                                                                                                                                                <  1:414658/251‑1 (MQ=255)                                            ttACCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGCTGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTacac                                                                                                                                                                                                          >  1:230792/1‑248 (MQ=255)                                              aCCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACt                                                                                                                                                                                                         <  1:2411/250‑1 (MQ=255)                                                cGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGgt                                                                                                                                                                                                      >  1:47186/1‑251 (MQ=255)                                                          tAAAACCAGGCGGCATTATGCTGCAATCAGTCTCCGTCGGTGAAACTATGGCGGCGATCACGTTGGTGCGCGCGCTGCGTCATAGTTATCCTGATTTACCGCTTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCAGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCACCCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATg                                                                                                                                                                                            >  2:262609/1‑251 (MQ=255)                                                           aaaaCCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTAt                                                                                                                                                                                             <  1:300529/249‑1 (MQ=255)                                                           aaaaCCAGGCGGCATTATGCTGCACTCCGTCTACGTCGGGGAAACTCTGGCGGCGATCCCGTTGGAGCGCGCGCTGCGTCATCGTTATCCTGTTTTACCGATTACCGTTACAACCACGACGCCAACCGGTTCGGATCGCGTACAATCGGCTTTCGGGAAGCATGTTCAGCACGTCTATCAGCCGTATGACCTGCCCGATGCGCTCAACAGTATCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATgg                                                                                                                                                                                           >  2:242195/1‑251 (MQ=255)                                                              aCCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTgg                                                                                                                                                                                                       >  2:80440/1‑236 (MQ=255)                                                              aCCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTgg                                                                                                                                                                                                       <  1:80440/236‑1 (MQ=255)                                                                          tATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGc                                                                                                                                                                            >  1:124556/1‑251 (MQ=255)                                                                           aTGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGcc                                                                                                                                                                           >  2:223410/1‑251 (MQ=255)                                                                                 cACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCt                                                                                                                                                                     <  1:298650/251‑1 (MQ=255)                                                                                 cACTCCGTCTCCGTCGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGAGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTAGGGCCTAACCt                                                                                                                                                                     <  2:60762/251‑1 (MQ=255)                                                                                             gtcGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCt                                                                                                                                                         <  2:259753/251‑1 (MQ=255)                                                                                               cGGTGAAACTCTGGCGGCGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTAc                                                                                                                                                       <  2:300270/251‑1 (MQ=255)                                                                                                ggTGAAAATATGGCGGCGCTCCCGTTGGTGAGCCCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTAc                                                                                                                                                       >  1:121099/1‑250 (MQ=255)                                                                                                         ctGGCGGCGATCTAGTTGGTGCGCGCGCTGCGTCATCGTTATCCGGATTTACCGATGACCGTTATAACCATGACGCCAACCGGTTTGGAGCGCGTACAATCGGCTTTCGGGAAGGATGGTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCGCAAACGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGTGGCGCTACATAAACGt                                                                                                                                               <  1:438734/249‑1 (MQ=255)                                                                                                            gcggcTATCCCGTTTGTGCGCGAGCTGCGTCATCGTTATCCGGATTTACCGATGTCCGTGACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCGGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTaatat                                                                                                                                          <  2:347177/251‑4 (MQ=255)                                                                                                             cggcGATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAAtt                                                                                                                                         >  1:173414/1‑251 (MQ=255)                                                                                                                  aTCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCt                                                                                                                                    >  1:32600/1‑251 (MQ=255)                                                                                                                       gTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATaa                                                                                                                                                                                                                        <  1:362476/162‑1 (MQ=255)                                                                                                                       gTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATaa                                                                                                                                                                                                                        >  2:362476/1‑162 (MQ=255)                                                                                                                                cgcgcTGCGTCATCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAAcg                                                                                                                      <  1:155923/251‑1 (MQ=255)                                                                                                                                           aTCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGGGATCGCTAACGCGCGActctct                                                                                                           >  1:243210/1‑251 (MQ=255)                                                                                                                                           aTCGTTATCCTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGCAAAATTCCGCTGGTGATCGCTAACGCGCGActctct                                                                                                           <  2:47186/251‑1 (MQ=255)                                                                                                                                                tagccTGAATTACCGATTACCGGTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGTCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGATAACGCGCGACTCTCTGccc                                                                                                       <  2:59721/247‑1 (MQ=255)                                                                                                                                                   ccTGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTc                                                                                                   <  1:266945/251‑1 (MQ=255)                                                                                                                                                     tGATTTACCGATTACCGTTACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCgg                                                                                                 >  1:27979/1‑251 (MQ=255)                                                                                                                                                                      ttacaacCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCaa                                                                                 >  1:373452/3‑250 (MQ=255)                                                                                                                                                                        acaacCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACt                                                                              >  1:279318/1‑251 (MQ=255)                                                                                                                                                                         caacCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACTCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTg                                                                             <  2:27979/251‑1 (MQ=255)                                                                                                                                                                          aacCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGACGCAGGTTATGCCAAACTgg                                                                            >  1:416088/1‑251 (MQ=255)                                                                                                                                                                                    ccAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTGGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTcgt                                                                  <  1:59154/251‑1 (MQ=255)                                                                                                                                                                                         cGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTcgtccgtc                                                             >  2:371931/1‑251 (MQ=255)                                                                                                                                                                                                     cgTACAATCGGCTGTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGcgtc                                                 <  2:11578/251‑1 (MQ=255)                                                                                                                                                                                                         cAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCGCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTAt                                             <  1:3465/251‑1 (MQ=255)                                                                                                                                                                                                           aTCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTAtt                                            >  1:207185/1‑250 (MQ=255)                                                                                                                                                                                                                  ttCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGAt                                    >  2:89128/1‑251 (MQ=255)                                                                                                                                                                                                                     gggAAGGATGTTCAGCACGTTTATCTGCCGTATGCCCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAACCGTAAAATTCCGCTGGTTATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCTTCTCTTGCTGCGTCGTATTACGCTGATtgc                                 <  2:265978/251‑1 (MQ=255)                                                                                                                                                                                                                       gAAGGATGTTCAGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATtgccg                               <  2:85052/251‑3 (MQ=255)                                                                                                                                                                                                                                  aGCACGTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAATgaa                    <  1:124840/251‑1 (MQ=255)                                                                                                                                                                                                                                       gTTTATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAATgaaga                  >  1:429876/1‑248 (MQ=255)                                                                                                                                                                                                                                          tATCTGCCGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAATGAAGAAGATg             <  1:89128/250‑1 (MQ=255)                                                                                                                                                                                                                                               gCCGTATGACCTGCCCGAGGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCGAACCTGATTGCGGCGTTACATAAACGTAAAATTCCGCTGGTGAGCGCTAACGCGCGACGCTCTGCCCGCTCGGCCGCAGGTGATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAATGAAGAAGATGGTg          <  2:209947/248‑1 (MQ=255)                                                                                                                                                                                                                                                ccGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAAt                       >  2:251609/1‑234 (MQ=255)                                                                                                                                                                                                                                                ccGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAAt                       <  1:251609/234‑1 (MQ=255)                                                                                                                                                                                                                                                ccGTATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAATGAAGAAGATGGTGCACg      >  2:408724/1‑251 (MQ=255)                                                                                                                                                                                                                                                   tATGACCTGCCCGATGCGCTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAATGAAGAAGATGGTGCACGttt   >  2:311668/1‑251 (MQ=255)                                                                                                                                                                                                                                                    aTGACCTGCCCGATGCGCTCAACCGTCTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGGTATGCCAAACTGGGAAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCCGCCCAAAATGAAGAAGATGGTGCACGtttt  >  2:177895/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       AGGCTCCGGCCTATCGAAAACGCTGGGGTGAACGTTACGGTTTTTACCGCCATCCGCTAAAACCAGGCGGCATTATGCTGCACTCCGTCTCCGTCGGTGAAACTCTGGCGGCAATCCCGTTGGTGCGCGCGCTGCGTCATCGTTATCCTGATTTACCGATTACCGTAACAACCATGACGCCAACCGGTTCGGAGCGCGTACAATCGGCTTTCGGGAAGGATGTTCAGCACGTTTATCTGCCGTATGATCTGCCCGATGCACTCAACCGTTTCCTGAATAAAGTCGACCCTAAACTGGTGTTGATTATGGAAACCGAACTATGGCCTAACCTGATTGCGGCGCTACATAAACGTAAAATTCCGCTGGTGATCGCTAACGCGCGACTCTCTGCCCGCTCGGCCGCAGGTTATGCCAAACTGGGTAAATTCGTCCGTCGCTTGCTGCGTCGTATTACGCTGATTGCTGCGCAAAATGAAGAAGATGGTGCACGTTTT  >  NC_000913/3808619‑3809112

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