Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,547,958	G→T	M186I (ATG→ATT)	fdnG →	formate dehydrogenase N subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,547,958	0	G	T	92.9%	170.5 / ‑0.4	56	M186I (ATG→ATT)	fdnG	formate dehydrogenase N subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base T (28/24);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATG‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGC  >  NC_000913/1547716‑1548190                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          aaCCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACTGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGt                                                                                                                                                                                                                                  >  1:395/1‑251 (MQ=255)  aCCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTGCCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTa                                                                                                                                                                                                                                 <  1:33278/251‑1 (MQ=255)     gTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAACTTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGAc                                                                                                                                                                                                                              >  2:157476/1‑251 (MQ=255)         gcgcTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAAcc                                                                                                                                                                                                                          <  2:190412/251‑1 (MQ=255)          cgcTACCCGGAATAGCGTGCGCCAGGTTCTGACAAATTGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGGATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACGTTATTGAAAAGAACGAGCAGGTCGTAACGGTAAACCGTTGGCTTTATACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGAGGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCa                                                                                                                                                                                                                         <  2:47047/251‑1 (MQ=255)               ccaggAAGATCGTGCGCCAGGTTTCGACAAAGTGCAGCGCAGTAGCTGGGAAGAAGTATTCTCCCGTGTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTAGTGAAAAGAACGAGCAGGGCGTACCGGTAAACCGTTGGCTGTCTACCGGTAGGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAACGTTGACCGCGCCCTCGGTATT‑CTGGCGGTAGACAACCAggtg                                                                                                                                                                                                                     <  2:130035/247‑3 (MQ=255)                   gAATATCGCGCGCCCGGTTCTGACAAATGGCCGCGCATTAGCTGGGCAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTACCGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCCACGAAACCGGGATCCTCACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGcgcgcgc                                                                                                                                                                                                                >  2:33516/1‑250 (MQ=255)                         cGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGaca                                                                                                                                                                                                          <  1:287340/251‑1 (MQ=255)                             cgcCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAACTTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACCCgg                                                                                                                                                                                                       >  2:43835/1‑250 (MQ=255)                               ccAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGGAGACAACCAGGCGCGCGTCTGACACGGAcc                                                                                                                                                                                                    >  1:191169/1‑251 (MQ=255)                                               ggCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCtt                                                                                                                                                                                    >  1:272141/1‑251 (MQ=255)                                                               ggAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTc                                                                                                                                                                    >  1:124170/1‑251 (MQ=255)                                                                                  aTTGCGAAGCTGATGAACGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAACAACGAAACCGGGATGCTGACCCAGAACTTTGCACGCTCCCTCGGGATT‑CTGGCGGTCGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCATGCTCCAACATTTGGTCGCGGGGCGATGaccaaccc                                                                                                                                                 >  2:334413/1‑250 (MQ=255)                                                                                  aTTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGaccaacc                                                                                                                                                  >  1:324171/1‑250 (MQ=255)                                                                                       gAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACCCGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTggg                                                                                                                                            >  2:8406/1‑251 (MQ=255)                                                                                                           gagACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAAAATTTGGTCGCGGTGCGATGACCAACCACTGGGGGGATATCAAAAACGCTAAc                                                                                                                        >  2:108567/3‑251 (MQ=255)                                                                                                                   aaCTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCTGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCgtgat                                                                                                                <  2:168863/251‑1 (MQ=255)                                                                                                                         aTTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCgtgatggtgat                                                                                                          <  2:329624/251‑1 (MQ=255)                                                                                                                           tGAAAAGAACGAGCAGGGCGTACCGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCCACCACTGGGTGGATATCAAAAAAGCTAACGTCgtgatggagctgg                                                                                                        >  2:224640/1‑245 (MQ=255)                                                                                                                            gAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATggg                                                                                                       <  1:165684/251‑1 (MQ=255)                                                                                                                             aaaaGAACGAGCAGTGCGGAACGGGAAACCGTGGGAGTTCCACCGGTATGCTGTGTGCCGCCGGAGCGAGCAACGGAACCGGTATTCCGGCCCAGAAATGTTTCCGCCTCATCGGG‑TGTCTGGCGGTAGACAACCAGGCGCGCGGCTGACACGGGCCAACGGTAGCAAGTCTTGCTTAAACATTTGGTAGCGGTGCGATGACCAGCCGCTGGGTGGAGATCAAAAACGCTAACGTCGTGATGGTGATGGGc                                                                                                      <  2:88941/251‑1 (MQ=255)                                                                                                                             aaaaGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAACTTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGATAACGTCGTGATGGTGATGGGc                                                                                                      >  2:327364/1‑251 (MQ=255)                                                                                                                             aaaaGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCCCTGGGGGGATATCCAAAACGCTAACGTCGGGATGGTGATGGGc                                                                                                      >  2:87218/1‑251 (MQ=255)                                                                                                                                        cAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACgctgc                                                                                           >  1:232691/1‑251 (MQ=255)                                                                                                                                         aGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCaaaaa                                                                                                                               >  2:285155/1‑214 (MQ=255)                                                                                                                                         aGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCaaaaa                                                                                                                               <  1:285155/214‑1 (MQ=255)                                                                                                                                         aGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACgctgct                                                                                          <  1:150978/251‑1 (MQ=255)                                                                                                                                           ggCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCCACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAACCGCTAAAGTCGTGATGGTGATGGGCGGTAACGctgctg                                                                                         >  2:170534/1‑250 (MQ=255)                                                                                                                                                   ggTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATc                                                                                >  1:65731/1‑251 (MQ=255)                                                                                                                                                       aaCCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGt                                                                            <  1:241006/251‑1 (MQ=255)                                                                                                                                                        aCCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTc                                                                           <  1:275134/251‑1 (MQ=255)                                                                                                                                                          cGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCgg                                                                         <  1:157476/251‑1 (MQ=255)                                                                                                                                                           gTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTCGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGGAACGCTGCTGAAGCGCATCCCGTc                                                                           >  1:339389/1‑248 (MQ=255)                                                                                                                                                           gTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGt                                                                        >  1:290299/1‑251 (MQ=255)                                                                                                                                                           gTTGGCTTTCTACCGGTATGCTGTGGGCCTCCGGGGCCAGCAACGAAAGCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATG‑TGGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGGCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCTTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTc                                                                           <  2:339389/248‑1 (MQ=37)                                                                                                                                                                   tCTACCGGTATGCTGTGTGCCTCCGGTGCCAGAAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGCCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGACTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCGCTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGTGCGTCCCGTCGGTTTCCGCTg                                                                >  1:233462/1‑251 (MQ=37)                                                                                                                                                                     tACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTAATGGTCGGGAACGCTGCTGAAGCGCATCCCGTAGGCTTCCGCTggg                                                              >  2:60418/1‑251 (MQ=255)                                                                                                                                                                                 gtgtgCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGc                                                    <  1:110974/249‑1 (MQ=255)                                                                                                                                                                                 gtgtgCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGGAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAACAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTACGCTGGGCGATgggag                                                     >  2:251033/1‑245 (MQ=255)                                                                                                                                                                                    tgCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGTTGGTGCTGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGCTGGAAGCGaaa                                                >  2:35169/1‑250 (MQ=255)                                                                                                                                                                                     gCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGaaaaa                                              <  2:335128/251‑1 (MQ=255)                                                                                                                                                                                     gCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGaaaa                                               >  1:144161/1‑250 (MQ=255)                                                                                                                                                                                               ccAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCCACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGt                                                                                                                                                               >  1:156812/1‑128 (MQ=255)                                                                                                                                                                                               ccAGCAACGAAACCGGGATGCTGACCCAGAAAGTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACGGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAacgac                                    <  1:224640/251‑1 (MQ=255)                                                                                                                                                                                                    aaCGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAAc                               >  1:271050/1‑251 (MQ=255)                                                                                                                                                                                                        aaaCCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCACCCTTg                           <  1:312912/251‑1 (MQ=255)                                                                                                                                                                                                        aaaCCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGGCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGTTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTg                           <  1:305204/251‑1 (MQ=255)                                                                                                                                                                                                          aCCGGGGTGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTTGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGAt                         <  2:269399/251‑1 (MQ=255)                                                                                                                                                                                                           ccGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCTGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATc                        <  1:327364/251‑1 (MQ=255)                                                                                                                                                                                                             gggATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGACT‑CTGGCGGTAGACAACCAGGCGCGCGGCTGCCACGGACCAACGGTAGCAAGTCGTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGGCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCATCCTTGATc                        <  2:88075/249‑1 (MQ=25)                                                                                                                                                                                                                    tGACCCAGAAATTTGCCCGCTCCCTCGGGATG‑CTGGCGGTAGACCACCAAGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGGTGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAaacaaacacga                                     >  2:163105/1‑223 (MQ=255)                                                                                                                                                                                                                    tGACCCAGAAATTTGCCCGCTCCCTCGGGATG‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAacga                                     <  1:163105/229‑1 (MQ=255)                                                                                                                                                                                                                    ggACCCAGTAATTTGCCCGATATCGCTGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCTGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGAt               <  2:339456/250‑1 (MQ=25)                                                                                                                                                                                                                     gACCCAGAAATTTGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATc              <  1:43835/251‑1 (MQ=255)                                                                                                                                                                                                                                 tGCCCGCTCCCTCGGGATT‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCAATGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCAGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAAAAACGACGCAACATTGATCGTTGCCGATCCCCGTTTTAcgc  >  1:237856/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCCCGCTCCCTGGGGATT‑CTGGGCGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTAcgc  <  2:114598/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                          AACCGTCTGCGCTACCCGGAATATCGTGCGCCAGGTTCTGACAAATGGCAGCGCATTAGCTGGGAAGAAGCATTCTCCCGTATTGCGAAGCTGATGAAAGCTGACCGTGACGCTAACTTTATTGAAAAGAACGAGCAGGGCGTAACGGTAAACCGTTGGCTTTCTACCGGTATGCTGTGTGCCTCCGGTGCCAGCAACGAAACCGGGATGCTGACCCAGAAATTTGCCCGCTCCCTCGGGATG‑CTGGCGGTAGACAACCAGGCGCGCGTCTGACACGGACCAACGGTAGCAAGTCTTGCTCCAACATTTGGTCGCGGTGCGATGACCAACCACTGGGTGGATATCAAAAACGCTAACGTCGTGATGGTGATGGGCGGTAACGCTGCTGAAGCGCATCCCGTCGGTTTCCGCTGGGCGATGGAAGCGAAAAACAACAACGACGCAACCTTGATCGTTGTCGATCCCCGTTTTACGC  >  NC_000913/1547716‑1548190

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