Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,117,707	G→A	G301S (GGT→AGT)	trhO →	tRNA U34 hydroxylase TrhO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,117,707	0	G	A	97.7%	129.1 / NA	44	G301S (GGT→AGT)	trhO	tRNA U34 hydroxylase TrhO
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/23);  minor base C (1/0);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CAAAGTCTGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAGGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAG  >  NC_000913/1117463‑1117953                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       aaaaGTCTGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACTTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCGGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAgttgtt                                                                                                                                                                                                                                                  <  2:414578/250‑1 (MQ=255)     gTCTGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAg                                                                                                                                                                                                                                              <  1:261969/251‑1 (MQ=255)         ggCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTga                                                                                                                                                                                                                                           <  1:294542/250‑1 (MQ=255)           cATATCGAGGGCGGAATTATTGAATTCGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACTAAAGTTGGTGTAGTGAGAt                                                                                                                                                                                                                                        <  2:416024/251‑1 (MQ=255)           cATATCGAGGGCGGAATTATTGAATACGCCCGTCAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCACAAATTTTGTTTTTGACGAGCGGACGGGCGAACGTATATCGGATGAGCTTATCGCGCATTGCCACCAGCGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAACTACAAAAGTTGTTGTAGTg                                                                                                                                                                                                                                            >  1:221684/1‑247 (MQ=255)              atCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTtg                                                                                                                                                                                                                                     <  1:609309/251‑1 (MQ=255)                               tGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGt                                                                                                                                                                                                                    >  1:502735/1‑251 (MQ=255)                                    aCGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCGGCGAAGAAAGCGCGTTAcc                                                                                                                                                                                                                <  2:567354/250‑1 (MQ=255)                                                cgcgcgAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGCTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTg                                                                                                                                                                                                                                            >  1:346363/1‑210 (MQ=255)                                                cgcgcgAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGCTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTg                                                                                                                                                                                                                                            <  2:346363/210‑1 (MQ=255)                                                        cAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGGTTTTGACGAGCGGATGGGAGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGAGCCGTGCGACAGCCATACCAACTGGAAAAATGATGGATGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTAGCGCCAGAAGAACAGCGACGCCg                                                                                                                                                                                           >  1:33773/1‑251 (MQ=255)                                                           gttctGCCGGTGCGTTTTATGGGCAAAAATTTTGTTTTTGACGATCGTATGGGCGAACGTATATCGGATGAGATTATCGCGCTTTGCCCCCAGTGCGGTGCGCCGTGCTACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTTGTTATTCAGTGTCCAGTATGCGCGGAAAAATTCAAAAGTTGTGGTAGTGAGATGTGCTTCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCg                                                                                                                                                                                        <  2:418553/247‑1 (MQ=255)                                                               tGCCGGTGCGTTTTATTGGCAAAAATTGTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCTGGc                                                                                                                                                                                     <  2:510752/250‑1 (MQ=255)                                                                            tATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAAt                                                                                                                                                                        <  1:436781/250‑1 (MQ=255)                                                                              ttGGCAAAAATTTTGTTTTTGCCGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGTGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGc                                                                                                                                                                     >  2:629566/1‑251 (MQ=255)                                                                                 gCAAAAATTTTGTTTTGGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAAGGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAAt                                                                                                                                                                  <  2:189666/251‑1 (MQ=255)                                                                                  aaaaaaTTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATa                                                                                                                                                                 <  1:15018/250‑1 (MQ=255)                                                                                    aaaaTTTTGTTTGTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGGGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGGATGCGCGGAAAAATACAAAAGTGGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAg                                                                                                                                                               <  1:99348/251‑1 (MQ=255)                                                                                     aaaTTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGa                                                                                                                                                              <  2:16477/251‑1 (MQ=255)                                                                                          ttGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCttt                                                                                                                                                         >  1:444235/1‑251 (MQ=255)                                                                                            gtgtttGACGAGCGGATGGGCGAACGTATAGCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGTGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTtaa                                                                                                                                                       <  1:629566/248‑1 (MQ=255)                                                                                                 tGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTCCCACCAGTGCGGCGCGCCGAGCGACAGCCATACCAACTATAAAAATGCTGCCTGCCACCTGCTGTTTATTCAGTGTCCAGCATGCGCGGACAAATACAAAAGATGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGCCGGGCAGGACGTGAAACTGGCAATAAGATCTTTAATAAgc                                                                                                                                                  >  1:399070/1‑250 (MQ=255)                                                                                                  gACGAGCGGATGGGAGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGAGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCCGTATGCGCGGAACAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAATCGCGTTAACGACAGAAGAACAGCGAAGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTATTAAGtc                                                                                                                                                 >  2:14974/1‑251 (MQ=255)                                                                                                     gAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCg                                                                                                                                              >  1:154901/1‑251 (MQ=255)                                                                                                     gAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCg                                                                                                                                              >  1:566285/1‑251 (MQ=255)                                                                                                        cGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAATCTCGTgg                                                                                                                                           >  2:25359/1‑251 (MQ=255)                                                                                                            ggggCGAACGTATATCGGATGAGAGTATCGCGCATTGCCACCAGTGCGGAGCGCGGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTAGTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGt                                                                                                                                       <  2:377253/250‑1 (MQ=255)                                                                                                                         atCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAGAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACATGAAAATGGCAATAAGATCTTTAATACGTCTCGGGGACCTCTGAAAACAacac                                                                                                                          >  1:278560/1‑251 (MQ=255)                                                                                                                            ggATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACtgt                                                                                                                       >  1:324192/1‑251 (MQ=255)                                                                                                                               tGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAACGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGc                                                                                                                     >  1:296977/1‑250 (MQ=255)                                                                                                                                   aTTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTcc                                                                                                                <  1:180752/251‑1 (MQ=255)                                                                                                                                   aTTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTcc                                                                                                                >  2:44699/1‑251 (MQ=255)                                                                                                                                                gCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGaa                                                                                                   <  1:73544/251‑1 (MQ=255)                                                                                                                                                 ccaccaGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGaat                                                                                                  >  2:449366/1‑251 (MQ=255)                                                                                                                                                 ccaccaGTGCAGCGCGACGGGAGACAGACATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTAAGTGTCCAGTATGCGCGGAAAAATCCAACAGTTGTTGTAGTGAGATTTGCAGCGAAGAAAGCGCGTTACCGCCAGAAGACCAGCGACGCCGGCGGGAAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGACTGAATACAAAACTGTGCATTCCTGATCCAACAGaat                                                                                                  >  2:250994/1‑251 (MQ=255)                                                                                                                                                                  cGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGAt                                                                                 >  1:530669/1‑251 (MQ=255)                                                                                                                                                                             cATACCAACTGTAAAAATGATGGCTGCCACATGCTGTTTATTCAGTGTCCCATATGCGCGGAAAAATACAACAGTTGTTATAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGCATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATcc                                                                      >  1:574130/1‑251 (MQ=255)                                                                                                                                                                                          aaaaTGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCAc                                                         <  1:202274/251‑1 (MQ=255)                                                                                                                                                                                                        cACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTTCTGCGAAGAAAGCGCGTGACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTTGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTctgctg                                           <  1:283698/251‑1 (MQ=255)                                                                                                                                                                                                                gTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCtt                                   <  1:352091/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaaaaaACAAAAGTTTTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTGACCGGCAGAAGAACAGCGACGCCGGCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTgg                                                                                                                                           <  1:322671/121‑1 (MQ=255)                                                                                                                                                                                                                                          aaaaaTACACAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAACATGGCAATAAGATCTTTCATAAGTCTCGTgg                                                                                                                                           >  2:322671/1‑121 (MQ=255)                                                                                                                                                                                                                                           aaaaTACAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCt        >  1:431654/1‑251 (MQ=255)                                                                                                                                                                                                                                                 cAAAAGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACGCTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAg  <  2:398447/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CAAAGTCTGGCATATCGAGGGCGGAATTATTGAATACGCCCGTAAGGCGCGCGAGCAGGGCTTGCCGGTGCGTTTTATTGGCAAAAATTTTGTTTTTGACGAGCGGATGGGCGAACGTATATCGGATGAGATTATCGCGCATTGCCACCAGTGCGGCGCGCCGTGCGACAGCCATACCAACTGTAAAAATGATGGCTGCCACCTGCTGTTTATTCAGTGTCCAGTATGCGCGGAAAAATACAAAGGTTGTTGTAGTGAGATTTGCTGCGAAGAAAGCGCGTTACCGCCAGAAGAACAGCGACGCCGTCGGGCAGGACGTGAAAATGGCAATAAGATCTTTAATAAGTCTCGTGGACGTCTGAATACAACACTGTGCATTCCTGATCCAACAGAATAAATATCATTGCCGGATGCGTGCCATCCGGCAACATTTCACGCTTACTTCTGCTGTACGCCTTCCACTGAAATAATCAGATCCACTTCCTGAGAAG  >  NC_000913/1117463‑1117953

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