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.5%	303.0 / ‑3.3	67	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (33/33);  minor base C (1/0);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793751‑1794040                                                                                                                                                   |                                                                                                                                                gcgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgt                                                                                                                                               >  1:211723/1‑147 (MQ=255)   gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTGCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                             >  2:719726/1‑145 (MQ=255)    accgCGGTCGGCGCTGCGGCAATCAATTTTTGATCCAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTcc                                                                                                                                           <  2:420346/146‑1 (MQ=255)     ctacGGCCGCCGCTGCGGCAACCAATTTTTGATACAGCGCATGGCGTCCTTCGCCATTAACTTCAATCTCACCAAACCCCCGCAACGCCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTGCATCGCTCCCCCGTTC‑TGTTCCa                                                                                                                                          <  2:438412/146‑1 (MQ=255)      tGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAg                                                                                                                                         >  1:848267/1‑149 (MQ=255)        gcagtcgcgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaa                                                                                                                                       >  1:1034557/7‑149 (MQ=255)         ggTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaaa                                                                                                                                      <  1:93343/149‑1 (MQ=255)         ggTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaaa                                                                                                                                      <  2:48424/149‑1 (MQ=255)          gTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCGTTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAGGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAAc                                                                                                                                     >  1:93029/1‑149 (MQ=255)           tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACt                                                                                                                                    >  1:754647/1‑149 (MQ=255)             ggCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTgg                                                                                                                                  <  2:74527/149‑1 (MQ=255)               cgcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGtt                                                                                                                                <  1:907393/149‑1 (MQ=255)                gcTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCCTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTg                                                                                                                               >  2:54696/1‑149 (MQ=255)                    cGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCCTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCTCGGTTC‑TGTTCCAGAAACTGGTTGCACg                                                                                                                           <  1:614154/149‑1 (MQ=255)                     ggCAATCAATTTTTGACACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACgg                                                                                                                          <  2:167725/149‑1 (MQ=255)                         atcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTCCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAAt                                                                                                                      >  2:574188/1‑149 (MQ=255)                         atcaatTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCACAGGTGGGTGTACAGTAAGGTTTAATCTCTTCATCGCTGCCCGGTTC‑GGTTCCAGAAACTGGTTGCACGGGagt                                                                                                                      >  2:819327/1‑147 (MQ=255)                          tcaatTTTTGATACAGCGGATGGCGTCCGTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATc                                                                                                                     >  1:46479/1‑149 (MQ=255)                                tttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCa                                                                                                               >  1:859092/1‑149 (MQ=255)                                 ttGATACAGCGGATGGCCGCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCAc                                                                                                              <  2:708900/149‑1 (MQ=255)                                      aCAGCGGATGGCGTCCTTCGCCGTTAACTTCAATCTTAATGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaa                                                                                                         >  2:879895/1‑149 (MQ=255)                                        aGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCAGGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                                                       >  2:1080888/1‑149 (MQ=255)                                          cGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAAcc                                                                                                     >  1:34622/1‑149 (MQ=255)                                            gATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtc                                                                                                   <  1:1063785/149‑1 (MQ=255)                                               ggCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACGCCCCCTGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGG‑TCCTGTTCCGGACACTGGTTCCCCAGGAATCCCAGCCCCATAAAAcaccgac                                                                                                >  2:936812/1‑143 (MQ=255)                                                      ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCa                                                                                         >  1:345113/1‑149 (MQ=255)                                                      ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                         >  1:629068/1‑149 (MQ=255)                                                                  ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                             >  2:679896/1‑149 (MQ=255)                                                                  ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                             >  2:923413/1‑149 (MQ=255)                                                                   tCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAata                                                                            <  2:1034557/149‑1 (MQ=255)                                                                    cAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTACGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGACTCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAatat                                                                           >  1:893830/1‑149 (MQ=255)                                                                     aaTCTTACTGAACATCGGGAACGTCACCCCCCACGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATAtt                                                                          <  2:341939/149‑1 (MQ=255)                                                                        cTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtctc                                                                       >  2:608987/1‑149 (MQ=255)                                                                          tACTGAACATCGTGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCACAGCCCCATAAAACCTCGAACGACCCAGGCTTTCGGACTCTTCTCCa                                                                     >  2:16093/1‑149 (MQ=255)                                                                               aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGc                                                                >  1:1002314/1‑149 (MQ=255)                                                                                ccATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCGAGAAACTGGTTGCACGGGAATCCCAGCCCCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                                                               <  2:386828/148‑1 (MQ=255)                                                                                 cATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCCTCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTc                                                              <  2:102365/149‑1 (MQ=255)                                                                                      ggAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTCCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATAtt                                                         >  1:369480/1‑149 (MQ=255)                                                                                        aaCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATCTTCACCAACTGCTCATAtttc                                                       >  2:517708/1‑147 (MQ=255)                                                                                         actgcACCCCCCATGTCGTGGAACAGCAAGTTTCAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTCCACGGGAATCCCAGCACCATAAAACCGCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                                                      <  2:356859/145‑1 (MQ=255)                                                                                         aCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCAGCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                                                      <  1:509629/149‑1 (MQ=255)                                                                                          cGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaaa                                                                                                                                      >  2:590462/1‑68 (MQ=255)                                                                                          cGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaaa                                                                                                                                      <  1:590462/68‑1 (MQ=255)                                                                                           gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCCCTTATTgcggc                                                    >  1:537002/1‑149 (MQ=255)                                                                                                 ccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCGGCACCATAAAGCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg                                              >  2:413361/1‑149 (MQ=255)                                                                                                   ccACGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCCACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGcc                                            <  1:39993/149‑1 (MQ=255)                                                                                                          ggggtACTGTAAGTTTCAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACttt                                     <  2:579285/147‑1 (MQ=255)                                                                                                             gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                                  <  1:216596/149‑1 (MQ=255)                                                                                                             gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                                  <  1:438097/149‑1 (MQ=255)                                                                                                             gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAg                                  <  2:920714/149‑1 (MQ=255)                                                                                                              tACAGTTAGTGTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAgg                                 >  2:829315/1‑149 (MQ=255)                                                                                                              tACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGCgg                                 >  1:240792/1‑149 (MQ=255)                                                                                                              tACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTGGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAgg                                 >  2:5623/1‑149 (MQ=255)                                                                                                               aCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGc                                <  1:823434/149‑1 (MQ=255)                                                                                                               aCAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTACGGCGTTAGGCCACCATTTGAggg                                >  2:535129/1‑148 (MQ=255)                                                                                                                  gTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAACCCTCGATCGACCCAGGCTTTCTGAGTATTCTCCAACTGCTCATATTGCGGCCTTAAGCCACACTTTGAGGCGAc                             >  1:80605/1‑149 (MQ=255)                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATCTTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                        >  1:29438/1‑149 (MQ=255)                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                        <  1:719726/149‑1 (MQ=255)                                                                                                                         tttaTCTCTTCATCGCTGCCCGGTTC‑TGATCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACa                      <  2:625919/146‑1 (MQ=255)                                                                                                                          tAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACaa                     >  2:654871/1‑149 (MQ=255)                                                                                                                           aaTCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCACATTGATTGGCGTGAAGCCACACTTTGAGGCGACATTGACAAt                    <  1:1080888/149‑1 (MQ=255)                                                                                                                              ctctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                 <  2:46479/149‑1 (MQ=255)                                                                                                                               ttgtCATCCCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAAc                <  2:982110/146‑1 (MQ=255)                                                                                                                                ctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACa               <  1:886199/149‑1 (MQ=255)                                                                                                                                       gCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACAtt        <  1:679896/149‑1 (MQ=255)                                                                                                                                       gCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGATTTCTGAAGATTCTCCAACTGACACCATCATCTCTTTAAGCCACCCCTCGAGACGTCAAGGACAATCAACAGCACAtt        >  1:790345/1‑149 (MQ=255)                                                                                                                                         tGCCCGGTTC‑TGGTCCCGAAACTGGTTGCACGGGAAGCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTAc      <  1:399480/149‑1 (MQ=255)                                                                                                                                            ccGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCgg   <  1:543291/149‑1 (MQ=255)                                                                                                                                             cGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc  <  1:879895/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGC  >  NC_000913/1793751‑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: ‑