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	98.3%	204.7 / NA	59	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 (31/27);  minor base C (1/0);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGG  >  NC_000913/2732748‑2733236                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     tAATGGTCTGCGTACCAGAAAGCGATGACTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGTCGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAgg                                                                                                                                                                                                                                                >  1:12075/1‑251 (MQ=255)     ggTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCCATCAGGTCg                                                                                                                                                                                                                                             >  2:311870/1‑250 (MQ=255)      gTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGAt                                                                                                                                                                                                                                           >  2:49305/1‑251 (MQ=255)              aCCAGAAAGCGATGACTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTTTATGCTGCTTCATCAATCAGGTCGATGGCTTTAt                                                                                                                                                                                                                                   <  1:353917/251‑1 (MQ=255)                        gATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTg                                                                                                                                                                                                                         <  1:6645/251‑1 (MQ=255)                                    gCTTTCCACGCTTCTGTTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTTAGCTGGATGATACGACGATCGAGTCGGTCGATTTCTTCTGGTTGTGAGTCAATCTGCATACGAATGCTGGATGCTGCGTCATCAATCAGATCGATGGGGTTATCCGGCATCTGACGGTCAGCAAt                                                                                                                                                                                                             <  2:185183/251‑1 (MQ=255)                                     cTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGATTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTAGAGTTCTTCTGGTTTTGAGTCAATCTTCATAAGAATGCTGGATGCTGCTTCTTCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATg                                                                                                                                                                                                            >  2:359374/1‑251 (MQ=255)                                         ccACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGc                                                                                                                                                                                                        <  1:44151/251‑1 (MQ=255)                                         ccACTCTTCCTCTAAATCGGAGTACTGACGTTCTTTGTCGCTCAGGTCTGCGTTGAGCATACCCAGACGTTTTGTACTGGCTACATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTGCTTCTGGTTTTGAGGCAATCTGCATACGAACGCTGGATGCGGTTGCACCAATCAGGTCCATGGCTGTATCCGGCAGCTTACGGGCAGCAATGTAGc                                                                                                                                                                                                        <  2:413583/251‑1 (MQ=255)                                          cACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCg                                                                                                                                                                                                       >  2:293284/1‑251 (MQ=255)                                                tcttctAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATAGCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGTTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATgag                                                                                                                                                                                                  <  1:280383/250‑1 (MQ=255)                                                   tctAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAAAGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACa                                                                                                                                                                                               >  2:379206/1‑250 (MQ=255)                                                       aCTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTc                                                                                                                                                                                          <  1:227157/251‑1 (MQ=255)                                                       aCTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTc                                                                                                                                                                                          >  2:221936/1‑251 (MQ=255)                                                        cTCGGCGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGtac                                                                                                                                                                                         >  2:361307/1‑249 (MQ=255)                                                          cGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGCCGGTCAGCAATGTAGCGATGAGACAACGTcgcc                                                                                                                                                                                       >  2:195038/1‑251 (MQ=255)                                                                   gACGTTCTTTTTCGCCCAGTTCTTCGTTGAGCAGATCCAGACGTTTTTTACGGGCTTCATCAGACTCTTTCAGTAACGCCTGTTGTTCCAGTTGGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTTGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGTTCGATGGCTTTATCCGGCAGCAGACGGTCAGCAATGTAGCGATGAGACAGCGTCGCCGCTGctat                                                                                                                                                                               <  2:22623/250‑4 (MQ=255)                                                                   gACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGcaac                                                                                                                                                                               >  2:158945/1‑250 (MQ=255)                                                                          tttGTCGCTCATGTCTTCGTTGAGCATATCCAGACGTTTTGTACTGGCTTCACCAGACTCTTTCATTAGCGCCTGTTGTTCCAGTTGGAGCTGGATGATACGACGAGCGATTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAGTTCTGGATGCTGCTTCAGCAAGCAGGTCGATGGCTTTAGCCGGCAGCTGACGGGCAGCAATGTAGCGATGAGACAACGTCGTCGCGGCAACAATTGcc                                                                                                                                                                        <  2:260967/250‑1 (MQ=255)                                                                          tttGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCg                                                                                                                                                                       >  2:364309/1‑251 (MQ=255)                                                                                    aGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAAt                                                                                                                                                             <  2:301423/251‑1 (MQ=255)                                                                                       tcttcGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAAttt                                                                                                                                                           <  1:273777/250‑1 (MQ=255)                                                                                            gTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACAt                                                                                                                                                     <  2:12261/251‑1 (MQ=255)                                                                                              tGAGCCTATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACAtgg                                                                                                                                                   >  2:182083/1‑251 (MQ=255)                                                                                                       ccAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGGTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGGTGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCCACAATTGCCGGGTCAGTAATTTGCCCATGGTGGTGCAAt                                                                                                                                          >  1:370845/1‑251 (MQ=255)                                                                                                         aGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTc                                                                                                                                        >  1:375181/1‑251 (MQ=255)                                                                                                                       gtctgcatcaGACTCTTTCATTAACGCCTGTTGGTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGGCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTTCAATTCGTAACGTTCttt                                                                                                                            <  2:379834/244‑1 (MQ=255)                                                                                                                                aGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGcc                                                                                                                 >  2:120363/1‑250 (MQ=255)                                                                                                                                   ctctTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGCGTCAATCTGCATACGAATGCTGGATGCTGCTTCATAAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGa                                                                                                               >  1:253527/1‑250 (MQ=255)                                                                                                                                               cGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTAt                                                                                                  >  2:185381/1‑251 (MQ=255)                                                                                                                                                    gttgttCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCtt                                                                                               <  1:379206/249‑1 (MQ=255)                                                                                                                                                       gttCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAAc                                                                                           >  2:65121/1‑250 (MQ=255)                                                                                                                                                           cAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAg                                                                                      <  1:49305/251‑1 (MQ=255)                                                                                                                                                                tGAGCGGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGGTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCAGCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGGCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTc                                                                                  <  2:211800/250‑1 (MQ=255)                                                                                                                                                                 gAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCgg                                                                                <  1:108471/251‑1 (MQ=255)                                                                                                                                                                     tGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaac                                                                            >  2:207807/1‑251 (MQ=255)                                                                                                                                                                     tGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGAAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATAGCCATGGTCTCTTCACCAGAAGGCTCGGcaac                                                                            >  2:281377/1‑251 (MQ=255)                                                                                                                                                                      ggaggatACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAAGCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGcaa                                                                             <  2:370845/245‑1 (MQ=255)                                                                                                                                                                       gatgatCCGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGGATCTTCAACAGAAGGCTCGGcaaca                                                                           >  1:326927/1‑250 (MQ=255)                                                                                                                                                                        atgatACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAATCAATCTGCATACGAAGGCTGGATGCTGCTTCATCAATCAGGTCTATGGCGTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTTGTGCAAGTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCaacaaa                                                                         <  2:139079/251‑1 (MQ=255)                                                                                                                                                                             acgacgAGCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACtt                                                                    <  1:215337/251‑1 (MQ=255)                                                                                                                                                                               gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGTATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACGCTTTc                                                                  <  2:253527/251‑1 (MQ=255)                                                                                                                                                                               gacgaTCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGGCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCCGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCCACAAACACTTTc                                                                  >  2:243374/1‑251 (MQ=255)                                                                                                                                                                                       aGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGc                                                                               >  1:161808/1‑230 (MQ=255)                                                                                                                                                                                       aGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTGTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacg                                                          <  1:251091/251‑1 (MQ=255)                                                                                                                                                                                       aGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCGGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAAGCGCAATGGTATCTTCAACAGAAGGCTCGGc                                                                               <  2:161808/230‑1 (MQ=255)                                                                                                                                                                                        gtcggtcgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacg                                                          >  2:319081/1‑250 (MQ=255)                                                                                                                                                                                         tcggtcgAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAacgac                                                        >  2:83045/1‑251 (MQ=255)                                                                                                                                                                                                aGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTcc                                                   >  2:335157/1‑249 (MQ=255)                                                                                                                                                                                                            tttGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCttt                                     >  2:347948/1‑251 (MQ=255)                                                                                                                                                                                                              tGAGTCAATCTGCATCCGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgcag                                            >  1:121821/1‑242 (MQ=255)                                                                                                                                                                                                              tGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCgcag                                            <  2:121821/242‑1 (MQ=255)                                                                                                                                                                                                                   cAATCTGCATACGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATg                              <  1:227641/251‑1 (MQ=255)                                                                                                                                                                                                                             aCGAATGCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCACTGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGat                    >  2:54876/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCTGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTGTTCAATGTACTGGCGATATTCg               <  2:177410/250‑1 (MQ=255)                                                                                                                                                                                                                                     tGGATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTc             >  1:352836/1‑250 (MQ=255)                                                                                                                                                                                                                                       gATGCTGCTTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCa            >  2:10232/1‑249 (MQ=255)                                                                                                                                                                                                                                          gctgctTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCTAGcg        <  2:359041/250‑1 (MQ=255)                                                                                                                                                                                                                                          gctgctTCATCAATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgt       >  1:276519/1‑251 (MQ=255)                                                                                                                                                                                                                                               atcatcaATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAAACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGcgtctcgg  >  2:272914/2‑247 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     TAATGGTCTGCGTACCAGAAAGCGATGCCTTCTCTGCTTTCCACTCTTCTTCTAACTCGGAGTACTGACGTTCTTTGTCGCTCAGTTCTTCGTTGAGCATATCCAGACGTTTTTTACTGGCTTCATCAGACTCTTTCATTAACGCCTGTTGTTCCAGTTTGAGCTGGATGATACGACGATCGAGTCGGTCGAGTTCTTCTGGTTTTGAGTCAATCTGCATACGAATGCTGGATGCTGCTTCATCGATCAGGTCGATGGCTTTATCCGGCAGCTGACGGTCAGCAATGTAGCGATGAGACAACGTCGCCGCTGCAACAATTGCCGGGTCAGTAATTTGCACATGGTGGTGCAATTCGTAACGTTCTTTCAGGCCACGCAGAATCGCAATGGTATCTTCAACAGAAGGCTCGGCAACAAACACTTTCTGGAAACGACGTTCCAGCGCAGCATCTTTTTCAATGTACTGGCGATATTCGTCAAGCGTCGTGG  >  NC_000913/2732748‑2733236

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