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.2%	248.9 / ‑2.3	58	coding (213/552 nt)	btuE	glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base . (30/26);  minor base A (0/1);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793751‑1794041                                                                                                                                                   |                                                                                                                                                 gcgACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCTTCATTTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgt                                                                                                                                                <  1:246038/149‑3 (MQ=255)   gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                              <  2:1052098/149‑5 (MQ=255)   gACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                              <  1:353774/149‑5 (MQ=255)   gACGGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCtgttc                                                                                                                                              <  2:124569/149‑5 (MQ=255)    aCTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTcc                                                                                                                                            >  2:909994/1‑149 (MQ=255)        cGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaa                                                                                                                                        >  1:507589/1‑149 (MQ=255)         ggTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGaaa                                                                                                                                       <  2:616711/149‑1 (MQ=255)           tCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGCAACt                                                                                                                                     >  2:420374/1‑149 (MQ=255)             ggCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTTCAGGAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTgg                                                                                                                                   <  2:320792/149‑1 (MQ=255)                     ggCAATCAATTTCTGATACAGCGGATGGCGGCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTT‑ATGTTCCAGAAACTGGTTGCACgg                                                                                                                           <  2:302321/149‑1 (MQ=255)                           caatTTTTGCTACAGCGGATGGCGTCCTTCTCCATCAACTTCAATCGTAATGCTCATCGGGCTCGTCTGCCCCCTCGCGCTCGTCCAGTAAGTTTTAATCTCCTCATCGCTGCTCGGCTC‑TGTTCCAGAAACTGGTTACACGGGAATcc                                                                                                                     <  2:206256/149‑1 (MQ=255)                              tttttGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAg                                                                                                                  <  2:883278/149‑1 (MQ=255)                                     tACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATa                                                                                                           >  1:739695/1‑149 (MQ=255)                                            gATTGCCTCCTTGCCCTTTAACTCCAATCTTCTTGAACATCGGGATCGTCTCCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCGTCACCGCTGCCCGTTTC‑TCTTCCAGAAACTGTTTGCACGGGAATCCCAGCACCATAAAACCtc                                                                                                    <  2:375308/149‑1 (MQ=255)                                              tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcga                                                                                                  <  1:121117/149‑1 (MQ=255)                                              tGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcga                                                                                                  <  1:1194989/149‑1 (MQ=255)                                               ggCGTCGTTCGCCATTAACTTCAATCTTACTGAACATCGGGAGCGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAGCTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgat                                                                                                 <  2:305805/149‑1 (MQ=255)                                                gCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCTGCACCATAAAACCtcgatc                                                                                                <  1:903032/149‑1 (MQ=255)                                                 cGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCCGACCCATAAAACCtgcatcg                                                                                               >  2:1165806/1‑149 (MQ=255)                                                 cGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCtcgatcg                                                                                               >  2:981675/1‑149 (MQ=255)                                                      ttCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                                          >  1:496126/1‑149 (MQ=255)                                                          ccATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCt                                                                                      <  2:739695/149‑1 (MQ=255)                                                                  ttCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAat                                                                              <  1:868470/149‑1 (MQ=255)                                                                       tCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtct                                                                         >  1:401779/1‑149 (MQ=255)                                                                       tCTTACTGAACATCGGGAACGTCACCCCCCATGTGGGGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGGatcgtct                                                                         >  2:1205863/1‑144 (MQ=255)                                                                        ctgaCTGACCATCGGGAACGTCGCCCCCCATGTGGTGGGACAGTAAGTTTTAGTCTCTTCAGCGCTGCCCGGTTC‑CGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGCATATtctc                                                                        <  2:686530/146‑1 (MQ=255)                                                                          tACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTCCTGAATATTCACCa                                                                      >  2:103896/1‑149 (MQ=255)                                                                            cTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCCCCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAAc                                                                    >  2:427627/1‑149 (MQ=255)                                                                             tGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGCTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATCTTCCCCAAcc                                                                   >  2:1186911/1‑148 (MQ=255)                                                                               aaCATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGGATATTCTCCAACTGc                                                                 >  1:1229714/1‑149 (MQ=255)                                                                                  aTCGGGAACGTCACCCTCCATGGGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGcccc                                                              <  2:1141158/149‑4 (MQ=255)                                                                                  aTCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGCATCCCAGCACCCTCAAACCTCGATCGACCCGGGCTTTCTGAATATTCTCCACCTGCTCa                                                              >  1:812070/1‑149 (MQ=255)                                                                                     gggAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCCACTGCTCatat                                                           >  1:616457/1‑149 (MQ=255)                                                                                      ggAACGTCA‑CCCCCATGTGGTGGTACAGTAAGAATTCATTCGATCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTg                                                         >  1:301708/1‑149 (MQ=255)                                                                                       gAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTg                                                         <  2:1018308/149‑1 (MQ=255)                                                                                           gTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCCGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcggc                                                     <  1:1097584/149‑1 (MQ=255)                                                                                                cccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATCTTCTCCAACTGCTCCTATTGCGGCCTTaa                                                >  2:1057414/1‑149 (MQ=255)                                                                                                 ccccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAg                                               <  2:507589/149‑1 (MQ=255)                                                                                                  cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGGATCCCAGCACCATAAAACCTCGATCGACCCAGGACTTTTACATATTCTCCAACTGCTCAGTTTGCGGCGTTAAGc                                              <  1:916097/149‑1 (MQ=255)                                                                                                  cccATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGc                                              <  1:644513/149‑1 (MQ=255)                                                                                                    cATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGccc                                            >  2:65048/1‑148 (MQ=255)                                                                                                            ggtACAGTAATTTTTAATCTCTCCACCGCTGCCCGGTTC‑TGTCCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCAGGCTTTACGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                    <  1:436363/149‑1 (MQ=255)                                                                                                            ggtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGa                                    <  1:54740/149‑1 (MQ=255)                                                                                                            ggtACAGGAAGTTTTAATCTCTTAATCGCTGCCCGGCTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGACCACATTGCGGCGTTAGGCCACACTTTGa                                    >  2:846862/1‑149 (MQ=255)                                                                                                             gtACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATCTTTCGCCGTTAAGCCACGCTTTCAg                                   >  1:862798/1‑149 (MQ=255)                                                                                                                   tAAGTTTTAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACCCTTTGAGGcacct                             >  2:621213/1‑145 (MQ=255)                                                                                                                     gattttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATACAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACAtt                           <  2:802939/147‑1 (MQ=255)                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACAttta                         >  1:50904/1‑147 (MQ=255)                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGa                         >  1:196050/1‑149 (MQ=255)                                                                                                                       ttttAATCTCTTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCACAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTACAACTGCTCATATTGCGGCGTTCAGCCACACTTTACGGCGACaatca                         >  2:342527/1‑145 (MQ=255)                                                                                                                              cttttCCTCGCTGCCCGGTTCTT‑TTCCAGAAACCCGTTGCACGGGAATCCCGCCACCATACGACCCCGATCGACCCAGGCTTTCCGAATACTCTCCAGCTGCTCATGTTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCac                  <  2:5767/146‑2 (MQ=255)                                                                                                                              ctctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCaa                  <  1:320513/149‑1 (MQ=255)                                                                                                                               tctTCATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCAGACTTTGAGGCGACATTGACAATCAAc                 <  2:536966/149‑1 (MQ=255)                                                                                                                                 ttCATCGCTGCCCGGCTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCCGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCCCGTATTGCGGCGTTAAGCCACACTTTGAGGCGACATGGAACATCAACAg               >  2:1122678/1‑149 (MQ=255)                                                                                                                                  gcATCGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGc              <  1:405510/148‑1 (MQ=255)                                                                                                                                   cATCGCTGCCCGGCTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGCATCTGCGCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGca             <  1:420374/149‑1 (MQ=255)                                                                                                                                      cGCTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACGTTGAGGCGACATTGACACTCAACAGCACAt          >  1:973973/1‑149 (MQ=255)                                                                                                                                        cTGCCCGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTa        >  1:328000/1‑149 (MQ=255)                                                                                                                                            ccGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACGCTTTGAGGCGACATTGACAATCAACAGCACATTACCgg    >  1:345031/1‑149 (MQ=255)                                                                                                                                             cGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc   <  2:1063438/149‑1 (MQ=255)                                                                                                                                             cGGTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGc   >  2:165824/1‑149 (MQ=255)                                                                                                                                              ggTTC‑TGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGAAATTGACAATCAACAGCACATTCCCGGCg  >  2:474332/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GCGACTGCGGTCGGCGCTGCGGCAATCAATTTTTGATACAGCGGATGGCGTCCTTCGCCATTAACTTCAATCTTACTGAACATCGGGAACGTCACCCCCCATGTGGTGGTACAGTAAGTTTTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATTGACAATCAACAGCACATTACCGGCG  >  NC_000913/1793751‑1794041

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