Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	421,439	T→C	F152L (TTC→CTC)	proY →	putative transporter ProY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	421,439	0	T	C	97.8%	139.8 / ‑4.0	45	F152L (TTC→CTC)	proY	putative transporter ProY
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (23/21);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

CCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGTCGTTGCAAATGGTGATGTTTGCTTAC‑GGTGGGATCGAAATTATCGGGATTACCGCCGGTGAAGCGAAAGATCCTGA  >  NC_000913/421192‑421671                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           cccGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGctct                                                                                                                                                                                                                                         >  1:182188/1‑251 (MQ=255)        gccAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCtt                                                                                                                                                                                                                                  >  1:282215/1‑251 (MQ=255)        gccAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCtt                                                                                                                                                                                                                                  <  2:156604/251‑1 (MQ=255)               ctTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGCCCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTc                                                                                                                                                                                                                           >  2:210354/1‑251 (MQ=255)                               tcAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTCCTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACcatcatcatc                                                                                                                                                                                                            >  1:50352/2‑250 (MQ=255)                                   gAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGAt                                                                                                                                                                                                       >  2:72247/1‑251 (MQ=255)                                       aCCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTc                                                                                                                                                                                                   <  1:278760/251‑1 (MQ=255)                                             gCCCGCTGGCAGGTTACCTTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGACATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGt                                                                                                                                                                                             >  2:370720/1‑251 (MQ=255)                                               ccGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGACATCCTTATTGTCGCAATCGCCGCTGTGCCCGCTTTTGGAATCTATCTGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGttt                                                                                                                                                                                           >  1:211381/1‑251 (MQ=255)                                                          ttACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGcatcatc                                                                                                                                                                                 >  1:113483/1‑250 (MQ=255)                                                                   cGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGACCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGa                                                                                                                                                                       <  1:341586/251‑1 (MQ=255)                                                                            cTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTg                                                                                                                                                                                                     >  2:143618/1‑212 (MQ=255)                                                                            cTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCGGGAATTCTGGCTCTCGTTCTTTAAAGGCGCCACCATCATCATCATGATTg                                                                                                                                                                                                     <  1:143618/212‑1 (MQ=255)                                                                                    ttGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTGTTGGTATCTATATGGGTGTCTGGTTCACGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTGCGGAATCATCATCTGGGGGATTGGCAACGGCGGGCa                                                                                                                                                       <  1:207753/250‑1 (MQ=255)                                                                                      gAAATCCTTATTGTCTCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCTCCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCCACGGCGGGCAacc                                                                                                                                                    <  1:45312/251‑1 (MQ=255)                                                                                       aactCCTTATTGTCAACATCGTCGATGTGAACGCTTTTGGTATCTATCTGGATGTCTGGATCCTGCAGGTGACGCACTGGATTTGGGTACTGAGAGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGACTTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGTTTGTCGCCGGTTTCGGCATAATCATCTGGGGGATTGGCAACGGCgggacaccc                                                                                                                                                   >  2:338160/4‑245 (MQ=255)                                                                                                            ccGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTgg                                                                                                                              <  2:59638/251‑1 (MQ=255)                                                                                                                 gtgACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACATGTGGAGCa                                                                                                                          >  1:117447/1‑250 (MQ=255)                                                                                                                     ccGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGGAAcggc                                                                                                                     <  2:342671/251‑1 (MQ=255)                                                                                                                       gCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCg                                                                                                                                                                                                                          <  2:333551/148‑1 (MQ=255)                                                                                                                       gCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCg                                                                                                                                                                                                                          >  1:333551/1‑148 (MQ=255)                                                                                                                         ttttGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACggcggcc                                                                                                                 >  1:221541/1‑250 (MQ=255)                                                                                                                          tttGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGct                                                                                                                 <  1:308104/250‑1 (MQ=255)                                                                                                                              gTATCTATATGGGTGTCTGGTACCCGACGGTGCAGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCCACCGACCCGGATTAATAACATGGGGGGCGACCGCGGCctcttc                                                                                                            >  2:345201/1‑251 (MQ=255)                                                                                                                                      atgggttgcTGGTCCCCGCCGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATCCGGTGAGTTGGAAGTCTGGCTCTCGTTCTGTAAAGTCGCCACCATCATCATCAGGATTGGCGCCGGTTCCGGCATCATCATCTGGGGTATTGGCAACTGCGGGCAACCGTCCTGTAGTCATACCCTGTGGAGCAACGGCGGCTTCTTCAGTAACg                                                                                                     <  1:338160/242‑1 (MQ=255)                                                                                                                                      atGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCAGTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGTCAACCGACCGGTATTCAGAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACgg                                                                                                    <  2:182188/251‑1 (MQ=255)                                                                                                                                        gggTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACggct                                                                                                  <  1:206447/251‑1 (MQ=255)                                                                                                                                               tGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGGATTCATAACCTGTGGGGCAACGGCGGCGTCTTCAGTAACggccgggtttt                                                                                           >  2:324888/1‑243 (MQ=255)                                                                                                                                                 gTTCCCGACGGTTCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCa                                                                                         <  2:24225/251‑1 (MQ=255)                                                                                                                                                   tCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCAt                                                                                        <  2:272351/250‑1 (MQ=255)                                                                                                                                                    cccGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCgggat                                                                                                                                           >  2:304824/1‑195 (MQ=255)                                                                                                                                                    cccGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTAt                                                                                                                                           <  1:304824/198‑1 (MQ=255)                                                                                                                                                       gACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATggga                                                                                    >  1:226493/1‑248 (MQ=255)                                                                                                                                                        aCGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AAt                                                                                 <  1:96841/251‑1 (MQ=255)                                                                                                                                                          ggTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGt                                                                               <  1:210354/251‑1 (MQ=255)                                                                                                                                                          ggTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGGATTCATAACCTGTGGAGCAACGGCGGCTTTTTCCGGAACggcttggttttgcttggtattt                                                                                >  2:215246/1‑233 (MQ=255)                                                                                                                                                               cGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGTCGTTg                                                                          <  1:320755/251‑1 (MQ=255)                                                                                                                                                                         gtGGGTACTGGGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGAGTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGTCGTTGCAAATGGTg                                                                 <  2:73538/249‑1 (MQ=255)                                                                                                                                                                                      gtggtgCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGACTTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACAGGTATTCATCACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCCTGGT‑CATGTCGGTTCAAATgggggtggtttctttcgg                                                   >  1:234591/1‑235 (MQ=255)                                                                                                                                                                                          tgCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGAATTCATAACCTGTGGAGAAAAGGAGGCTTCTTCAGTAACGGCTGGCTTTGGATGGGGATTGT‑GTTGCAAATGGGGATGTTTGCTTACCGGTgg                                               >  1:219171/1‑251 (MQ=255)                                                                                                                                                                                           gCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGGATGGG‑AATGTCGTTGCAAATGGGGATGTTTGCTTAC‑GGTgggg                                             >  1:297113/1‑250 (MQ=255)                                                                                                                                                                                            cTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGG‑AATGTCGTTGCAAATGGTGATGTTTTCTTAC‑GGTggggt                                            >  1:70751/1‑249 (MQ=255)                                                                                                                                                                                                      gcgcCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGTCGTTGCAAATGGTGATGTTTGCTTAC‑GGTGGGATCGAAATTAt                                   <  2:282215/250‑1 (MQ=255)                                                                                                                                                                                                           gTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGTCGTTGCAAATGGTGATGTTTGCTTAC‑GGTGGGATCGAAATTATCGGGAt                             <  2:365385/251‑1 (MQ=255)                                                                                                                                                                                                                                gTATTCGGTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCCTTGT‑AATGTCGTTGCAAATGGGGATGTTTGGTTAC‑GGTGGGATCGAAATTTTCGGGATTACCGCCGGGGAAGCGaa           >  1:240189/1‑248 (MQ=255)                                                                                                                                                                                                                                      ggTGAGCTGGAATTCTGGCTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGTCGTTGCAAATGGTGGTGTTTGCTTAC‑GGGGGGGTCGGAATTATCGGGATTACCGCCGGTGAAGCGGAAGATCCTgg  >  1:212777/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           CCCGGCCGCCAGCTCTTTCTCGCGTTATGCGCAGGAAAACCTCGGCCCGCTGGCAGGTTACATTACCGGCTGGACCTACTGCTTTGAAATCCTTATTGTCGCCATCGCCGATGTGACCGCTTTTGGTATCTATATGGGTGTCTGGTTCCCGACGGTGCCGCACTGGATTTGGGTACTGAGCGTGGTGCTGATCATTTGCGCCGTAAACCTGATGAGCGTGAAGGTATTCGGTGAGCTGGAATTCTGGTTCTCGTTCTTTAAAGTCGCCACCATCATCATCATGATTGTCGCCGGTTTCGGCATCATCATCTGGGGGATTGGCAACGGCGGGCAACCGACCGGTATTCATAACCTGTGGAGCAACGGCGGCTTCTTCAGTAACGGCTGGCTTGGCATGGT‑AATGTCGTTGCAAATGGTGATGTTTGCTTAC‑GGTGGGATCGAAATTATCGGGATTACCGCCGGTGAAGCGAAAGATCCTGA  >  NC_000913/421192‑421671

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