Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,732,992	G→A	I394I (ATC→ATT)  I394I (ATC→ATT)	clpB ← clpB ←	chaperone protein ClpB chaperone protein ClpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,732,992	0	G	A	97.4%	128.2 / NA	39	I394I (ATC→ATT)  I394I (ATC→ATT)	clpB clpB	chaperone protein ClpB chaperone protein ClpB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/16);  minor base C (1/0);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCG  >  NC_000913/2732746‑2733230                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               tttAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTAAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGCGTTCTTCTGGTTTTGAGTCAATCTGCATACGCATGCTGGATGCTGCTTCATCAATc                                                                                                                                                                                                                                             >  1:419152/1‑250 (MQ=255)            tcGTACCAGAAAGCGATGCCTTCTCTGCGTTCCACTCTTCTTCTAACTCGTAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACGGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCTGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATg                                                                                                                                                                                                                                    <  1:251335/247‑1 (MQ=255)                         cGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCt                                                                                                                                                                                                                    <  1:125106/251‑1 (MQ=255)                          gATGCCTTCTCTGCTTTCCCCTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATAAGAATGCTGGATGCTGCTTCATCAATCAGGTCGCTGGCTTTATCCGGCAGCTg                                                                                                                                                                                                                   >  1:277887/1‑251 (MQ=255)                          gATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTg                                                                                                                                                                                                                   >  1:227645/1‑251 (MQ=255)                                       cTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATg                                                                                                                                                                                                      >  1:305714/1‑251 (MQ=255)                                                         aCTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTc                                                                                                                                                                                    >  1:408834/1‑251 (MQ=255)                                                            cGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTcgcc                                                                                                                                                                                 <  1:198295/251‑1 (MQ=255)                                                                        gTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAAtt                                                                                                                                                                     <  1:340288/251‑1 (MQ=255)                                                                                     cAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGt                                                                                                                                                          <  2:408834/249‑1 (MQ=255)                                                                                              gTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACAt                                                                                                                                               <  2:150243/251‑1 (MQ=255)                                                                                                 gAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACAtggt                                                                                                                                            <  1:14530/251‑1 (MQ=255)                                                                                                   gCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCTGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACACTTGACGGGTCAGGAATTGGCACAtggtgg                                                                                                                                          >  2:120053/1‑251 (MQ=255)                                                                                                     atatCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGCCGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTTCCGGGTCAGTAATTTGCACATGgtggtg                                                                                                                                        >  2:17859/1‑251 (MQ=255)                                                                                                       atCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGGCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAAGTTGCACATGGTGGTGCa                                                                                                                                      <  2:84212/251‑1 (MQ=255)                                                                                                           aGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGCTGGCTTTATCCGGCAGCTGACGGTCAGCACTGTAGCGATGAGACACCGTCGCCGCTGCAACAATTTCCGGGTCAGTACTTTGCACATGGTGGTGCAATTc                                                                                                                                  >  2:194803/1‑251 (MQ=255)                                                                                                                   tttACTGGCTTCCTAAGACTCTTTCCTTCACGCCTGTTGTTCCAGTTTGGGCTGGATGATACGACGCTCGCGTCGGTCGAGATCTTCTGGTTTTGAGTCAATATGCATACGAATGCTGGATGCTGCTTCATCCATCAGGTCGCTGGCTTTATCCGGCAGCTGCCGGTCCGCAATGTATCGATCACACACCGTCGCCGCTGCCACAATTGCCGGGCCAGTAATTTGCACATGGTGGTGCAATTCGTaccggt                                                                                                                          >  2:187618/1‑246 (MQ=255)                                                                                                                      aCTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGCGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACCATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGttatt                                                                                                                       >  1:399389/1‑248 (MQ=255)                                                                                                                      aCTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCtt                                                                                                                       >  1:252864/1‑251 (MQ=255)                                                                                                                      aCTGGCTTCATCAGAATCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTTGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGACGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGGGGTGCAATTCGTAACGTTCtt                                                                                                                       >  2:88549/1‑251 (MQ=255)                                                                                                                          gCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGcaac                                                                                                                                                                         >  2:90856/1‑197 (MQ=255)                                                                                                                          gCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGcaac                                                                                                                                                                         <  1:90856/197‑1 (MQ=255)                                                                                                                          gCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAg                                                                                                                   <  1:416426/251‑1 (MQ=255)                                                                                                                            ttCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCACTCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCACTGTAGCGATGAGACAACGTCGCCGCTGCCACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGc                                                                                                                 >  1:409775/1‑251 (MQ=255)                                                                                                                                tcaGACTCTTTAATTAACGCCTGTTGTTCCAGCTTGCGCTGGATGATACGACGCTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACCTGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACg                                                                                                             >  2:244001/1‑251 (MQ=255)                                                                                                                                      tctTTCATTAAAGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAAt                                                                                                       >  1:319814/1‑251 (MQ=255)                                                                                                                                          tCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGCATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGc                                                                                                    >  1:218533/1‑250 (MQ=255)                                                                                                                                                  gCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATc                                                                                           >  1:312953/1‑251 (MQ=255)                                                                                                                                                  gCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATc                                                                                           >  1:274559/1‑251 (MQ=255)                                                                                                                                                                                     aTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGa                                                        <  2:171906/251‑1 (MQ=255)                                                                                                                                                                                          gtcggtcgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGACTCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacg                                                    >  2:362887/1‑250 (MQ=255)                                                                                                                                                                                            cagtcgAGTTCTTCTGGTTTTGAGTCCATCTGCATACGAATGCTAGAGGCTGCTTCATCAATCAGGTCGATGGCTTTACCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCACTAATTTACACATGGTGGTGCAATTCGTAACTTTCTTTCAGGCCACGCCGAATCGCAATGGTAGCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacgac                                                  <  2:99742/248‑1 (MQ=255)                                                                                                                                                                                             ggtcgCGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGt                                                >  2:62045/1‑251 (MQ=255)                                                                                                                                                                                                cgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCGCGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTcc                                             <  1:199977/251‑1 (MQ=255)                                                                                                                                                                                                           gtttttTAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAAACAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATc                                  <  1:271823/249‑1 (MQ=255)                                                                                                                                                                                                                                 gAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGatat            >  1:338932/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaTGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGGCGTTCCAGCGCAGCATCTTTTTCAATTTTCTGGCGATAtt           <  2:149459/251‑1 (MQ=255)                                                                                                                                                                                                                                   aTGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTc          <  1:405195/251‑1 (MQ=255)                                                                                                                                                                                                                                           tgctgcTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcg  >  2:419312/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               TTTAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCG  >  NC_000913/2732746‑2733230

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