Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,793,898	Δ1 bp	coding (212/552 nt)	btuE ←	glutathione peroxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,793,897	0	T	.	98.9%	430.6 / ‑3.1	94	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (44/49);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793754‑1794040                                                                                                                                                |                                                                                                                                                aCTGCGGTCCGCCCCGCGGCAACCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTcc                                                                                                                                           <  1:100904/149‑1 (MQ=255)    gCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACGCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGa                                                                                                                                        <  2:1148239/149‑1 (MQ=255)     cGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaa                                                                                                                                       >  2:512585/1‑149 (MQ=255)        tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACt                                                                                                                                    >  1:1091737/1‑149 (MQ=255)        tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGACGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACt                                                                                                                                    >  2:936837/1‑149 (MQ=255)          cgcgcTGCGGCCACCATTTTTTGATAGCGCGGATGGCGTGCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTATAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTgg                                                                                                                                  <  1:56219/148‑1 (MQ=255)            cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGtt                                                                                                                                <  2:1463217/149‑1 (MQ=255)                  ggCAATCAATTTTTGATACAGCGGAGGGCGGCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACgg                                                                                                                          >  1:802306/1‑149 (MQ=255)                  ggCAATCAATTTTCGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACgg                                                                                                                          <  2:139006/149‑1 (MQ=255)                   gcgatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACggg                                                                                                                         >  2:413787/4‑149 (MQ=255)                    caatcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGa                                                                                                                        <  2:1175867/149‑1 (MQ=255)                      atcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCCCGGCAAt                                                                                                                      >  2:386923/1‑149 (MQ=255)                        caatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTAAAGCTTTTTGTGCCTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATcc                                                                                                                    <  1:1062473/149‑1 (MQ=255)                         aatTTTTGCTACAGCGGATGGCGTCCGTCGCCATTAACTTCAAGCTTACTGAACCTCGGGAACGTCACCCCCCTGGTGGTGGTACAGTAAGGTTTAATCTCTGCATCGCTGCCCGGTTC‑TGTTCCCGAAACTGGTTGCACGGGAATccc                                                                                                                   <  1:365688/149‑1 (MQ=255)                         aatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCAAGAAACTGGTTGCAAGGGAATccc                                                                                                                   <  1:150854/149‑1 (MQ=255)                         aatTTTCGATACCGCGGATGGCGTCCTTCGCCATTAACGTCCATCTTACTGAACATCCCGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTTTTCATCGCTTCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATccc                                                                                                                   <  1:295448/149‑1 (MQ=255)                         aatTTTCGATACAGCGGATGGCGTCCTTCGCCATTAACTGCAACCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATccc                                                                                                                   <  2:802306/149‑1 (MQ=255)                           tttttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGCTGCACGTGAATCCCac                                                                                                                 >  2:33947/1‑148 (MQ=255)                             tttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGcg                                                                                                               >  1:926685/1‑148 (MQ=255)                             tttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCACCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCa                                                                                                               <  1:934991/149‑1 (MQ=255)                                 aTACAGCGGATGGCGTCCTCTGCCATTAACTTGAATCTAACCGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCAt                                                                                                           <  1:810445/149‑1 (MQ=255)                                   aCAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                                                                         >  1:426792/1‑149 (MQ=255)                                    cAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAACCTGGTTGCACGGGCATCCCAGCACCATaaa                                                                                                        >  1:76710/1‑149 (MQ=255)                                     aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaac                                                                                                       >  1:405712/1‑148 (MQ=255)                                     aGCGGATGGCGTCCTTCGCCACTAACTTCAATCTTCCTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                       <  2:1556707/149‑1 (MQ=255)                                        ggATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACCTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCt                                                                                                    <  1:1336200/149‑1 (MQ=255)                                        ggATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCTCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCt                                                                                                    <  1:1555186/149‑1 (MQ=255)                                            ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                                <  1:512585/149‑1 (MQ=255)                                                  cTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAccc                                                                                          <  1:140336/149‑1 (MQ=255)                                                  cTTCGCCATTAACTTCAACCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAccc                                                                                          <  1:168948/149‑1 (MQ=255)                                                   ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGGTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                         >  1:539581/1‑149 (MQ=255)                                                   tccGCCATTAACTGCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                         >  2:1179911/3‑149 (MQ=255)                                                    tCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAg                                                                                        <  1:585221/149‑1 (MQ=255)                                                        cATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCtt                                                                                    >  2:63512/1‑149 (MQ=255)                                                          ttAACTTCCATCTTCCGTAACCTCGGGAACGTCACCCCCCCTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTc                                                                                  <  2:1299107/149‑1 (MQ=255)                                                           tAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCt                                                                                 >  1:816462/1‑149 (MQ=255)                                                            aaCTTCAATCTTACTGAACATCGGGATCCTCACCCCCCTTGTGGTGGGACAGTAAGTTTGAATGTCGGCATCGCTGCCCGGTTCTT‑TTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTg                                                                                <  2:524337/149‑1 (MQ=255)                                                               ttCAATCTTCCCGAACATCGGGAACCTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAATTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                             <  2:1007747/149‑1 (MQ=255)                                                               ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGATACGGGAGGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                             >  2:1171845/1‑149 (MQ=255)                                                                 cAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGActat                                                                           >  1:61782/1‑149 (MQ=255)                                                                   aTCTTCCTGAACCTCGTGAACGTCACCTCCCATGTGGTGGTACAGTAAGTTTTACTCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                                                         <  2:671171/149‑1 (MQ=255)                                                                    tCTTACTGAACACCGGGACCGTCACCCCCCATGTGGTGGTACAGTAAGATTTAATCTCTTCATCGCTGCCCGGTTC‑CGTTCCAAAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtct                                                                        >  2:412080/1‑149 (MQ=255)                                                                     cTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtctc                                                                       <  2:802637/149‑1 (MQ=255)                                                                     cTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtctc                                                                       >  1:950612/1‑149 (MQ=255)                                                                                 cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGGAACTGGTTGCACGGGCATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCata                                                           >  1:27959/1‑149 (MQ=255)                                                                                 cGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCCCata                                                           >  2:1277827/1‑149 (MQ=255)                                                                                     aaCTTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCGGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgc                                                       <  1:1169442/149‑1 (MQ=255)                                                                                     aaCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgc                                                       <  1:884946/149‑1 (MQ=255)                                                                                     aaCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgc                                                       <  1:76562/149‑1 (MQ=255)                                                                                     aaCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGAGCGACCCAGGCTGTCTGAATATTCTCCAACTGCTCCTATTgc                                                       >  1:347771/1‑149 (MQ=255)                                                                                      aCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                                                      <  2:654376/149‑1 (MQ=255)                                                                                      aCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                                                      <  1:936837/149‑1 (MQ=255)                                                                                        gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCGGCCCCATAAAACCTCGATCGACCCCGGCGTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                    >  2:504846/1‑149 (MQ=255)                                                                                        gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                    >  1:990066/1‑149 (MQ=255)                                                                                        gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                    >  1:966207/1‑149 (MQ=255)                                                                                         tCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggcg                                                   >  1:221083/1‑149 (MQ=255)                                                                                         tCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCCCCAACTGCTCATCTTgcggcg                                                   >  2:359740/1‑149 (MQ=255)                                                                                         tCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGAATTCTGAATAGACTCCAACTGATAATATTcctgcg                                                   >  1:354669/1‑149 (MQ=255)                                                                                             cccccATGTGGTGGTACAGTAATTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                               <  2:1291534/149‑1 (MQ=255)                                                                                             cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGGATCCCAGCACCATAAAACCTCGATCGACCCCGGCTTTCTGAATATTCTCCAACTGCCCATATTacggcgcgaa                                               >  1:646327/1‑145 (MQ=255)                                                                                             cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGAGCCACGCTGTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTaa                                               <  1:988870/149‑1 (MQ=255)                                                                                              ccccATGTGGGGGTACAGTACGTTTTCATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCGTATTGCGGCGTTAAg                                              <  1:183378/149‑1 (MQ=255)                                                                                               cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGc                                             <  1:1364205/149‑1 (MQ=255)                                                                                                 cATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGGGTTAAGCca                                           <  2:408135/149‑1 (MQ=255)                                                                                                 cATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca                                           <  2:990066/149‑1 (MQ=255)                                                                                                      ggtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACtt                                      >  1:1153685/1‑149 (MQ=255)                                                                                                       gtggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                                     <  2:1153685/149‑1 (MQ=255)                                                                                                         ggtACAGTAAGTTTTAATCTCTTCCTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                   <  2:1091737/149‑1 (MQ=255)                                                                                                          gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                                  >  1:343951/1‑149 (MQ=255)                                                                                                          ggaCAGTAAGTTTTAATCTCGTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                                  <  1:599908/147‑1 (MQ=255)                                                                                                            aCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑GGGTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGc                                <  2:82155/149‑1 (MQ=255)                                                                                                              aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCGACACTTTGAGGCgc                              >  1:22595/1‑148 (MQ=255)                                                                                                              aGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGAGAAACCCATCACCATAAAACCTCGATCGACCAACGCTCCCCAAACATCCGCCAACTTCTCATATTGCGGCGTTAAGCCACACCTTGAGGCGa                              >  2:1098427/1‑149 (MQ=255)                                                                                                                tAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAc                                                                                                              >  2:1404543/1‑67 (MQ=255)                                                                                                                tAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAc                                                                                                              <  1:1404543/67‑1 (MQ=255)                                                                                                                   gTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTc                                                              >  2:680250/1‑112 (MQ=255)                                                                                                                   gTTTTAATCTCTTCAACGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTc                                                              <  1:680250/112‑1 (MQ=255)                                                                                                                    ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                        <  1:352815/149‑1 (MQ=255)                                                                                                                    ttttAATCTCTTCATCGCTGCCCGGATC‑CGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                        >  1:871781/1‑149 (MQ=255)                                                                                                                       tAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACaa                     >  1:913452/1‑149 (MQ=255)                                                                                                                         aTCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGCAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATc                   <  1:1274173/149‑1 (MQ=255)                                                                                                                          tctcTTCATCGCTGCCCGGTTC‑TGTTCCCGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCa                  >  1:455613/1‑149 (MQ=255)                                                                                                                           ctctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                 <  1:182938/149‑1 (MQ=255)                                                                                                                            tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc                <  2:816462/149‑1 (MQ=255)                                                                                                                             ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa               >  2:207902/1‑149 (MQ=255)                                                                                                                             ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa               <  2:495284/149‑1 (MQ=255)                                                                                                                             ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTATAAGTCGACCTTGCCCATCcaca               >  2:891699/1‑149 (MQ=255)                                                                                                                             cgtCATCGCTGCCCGGTTC‑TGTTCCAGAGACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGACTATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa               <  1:634207/147‑1 (MQ=255)                                                                                                                              gtactCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTCCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAg              <  2:1154818/145‑1 (MQ=255)                                                                                                                                 aTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGcac           >  1:1233469/1‑149 (MQ=255)                                                                                                                                         ccGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACCCTTTGAGGCGACATTGACAATAAACAGCACATTACCgg   >  2:281493/1‑149 (MQ=255)                                                                                                                                          cGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGCCAATCAACACCACATTACCGGc  >  2:70161/1‑149 (MQ=255)                                                                                                                                          cGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc  >  2:199003/1‑149 (MQ=255)                                                                                                                                          cGGTTC‑TGTCCCAGAAACGGGTTGCACGGGAATCCCAGCAACGTAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCGCATATTTCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc  <  2:137083/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                ACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793754‑1794040

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