Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,909,297	A→T	10.5%	intergenic (‑77/+92)	spuD ← / ← spuC‑II	spermidine/putrescine ABC transporter substrate‑binding protein/polyamine:pyruvate transaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,909,297	0	A	T	10.5%	178.3 / 8.4	64	intergenic (‑77/+92)	spuD/spuC‑II	spermidine/putrescine ABC transporter substrate‑binding protein/polyamine:pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (4/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGT  >  NC_002947/5909151‑5909443                                                                                                                                                   |                                                                                                                                                   cAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGaaa                                                                                                                                                    <  5:92777/147‑1 (MQ=255)  aGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAACCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACa                                                                                                                                                  <  6:328080/148‑1 (MQ=255)     cTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGTTTCAGTCg                                                                                                                                              >  8:238767/1‑149 (MQ=255)        caCAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGc                                                                                                                                             <  8:121795/146‑1 (MQ=255)         aCAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGTTTCAGTCGccgc                                                                                                                                          <  7:238767/149‑4 (MQ=255)            gCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCt                                                                                                                                       <  1:333944/149‑1 (MQ=255)                 ggCCATGGCATGCATCAGGGCGGCGGTCAGCAAGGTCTTGGCCATTTTATTCATGCGTAATGCTCCAGTATTTTTATTTAGCTTCTGTCGATCAGGGTGATCCCAATGGGATGAGCGGGTGTCGACTGATTGAGTCGCGAGTCTgtgaa                                                                                                                                  <  6:69901/149‑5 (MQ=255)                  gCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAg                                                                                                                                 >  3:226548/1‑149 (MQ=255)                   ccATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                <  4:226548/149‑1 (MQ=255)                    cATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                >  5:111893/1‑148 (MQ=255)                     aTGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCAAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTa                                                                                                                              >  1:285822/1‑149 (MQ=255)                            ccATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGcc                                                                                                                       >  2:128381/1‑149 (MQ=255)                                  gggCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCTTGTAt                                                                                                                 >  4:63466/1‑149 (MQ=255)                                      ggcggcCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAg                                                                                                             <  8:114308/149‑1 (MQ=255)                                                aaCGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAggg                                                                                                   <  6:181990/149‑1 (MQ=255)                                                   gTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCa                                                                                                >  5:41883/1‑149 (MQ=255)                                                    tCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCTTGTATCAAGGAAAGCAGGGCCAt                                                                                               <  3:63466/149‑1 (MQ=255)                                                    tCTTGCCCATTTTATGCATGCGTACTGCTGCTGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGTTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCAt                                                                                               <  6:307934/149‑1 (MQ=255)                                                      ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCaa                                                                                                                                  >  3:200091/1‑112 (MQ=255)                                                      ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCaa                                                                                                                                  <  4:200091/112‑1 (MQ=255)                                                      ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCaa                                                                                                              <  7:322309/132‑1 (MQ=255)                                                      ttGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCaa                                                                                                              >  8:322309/1‑132 (MQ=255)                                                        gCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTaa                                                                                           >  2:268618/1‑149 (MQ=255)                                                        gCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTa                                                                                            <  6:41883/148‑1 (MQ=255)                                                           cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGt                                                                                                                                                <  6:232430/93‑1 (MQ=255)                                                           cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGt                                                                                                                                                >  5:232430/1‑93 (MQ=255)                                                           cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGa                                                                                        >  4:146201/1‑149 (MQ=255)                                                            aTTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          >  7:177001/1‑146 (MQ=255)                                                            aTTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          <  8:177001/146‑1 (MQ=255)                                                                 aTTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTaa                                                                                   <  1:268618/148‑1 (MQ=255)                                                                     aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                <  7:80900/99‑1 (MQ=255)                                                                     aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                >  8:80900/1‑99 (MQ=255)                                                                     aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTaa                                                                                   >  2:138246/1‑144 (MQ=255)                                                                     aTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTaa                                                                                   <  1:138246/144‑1 (MQ=255)                                                                             tGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGc                                                                        >  3:64393/1‑147 (MQ=255)                                                                                tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          >  5:104148/1‑126 (MQ=255)                                                                                tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAAc                                                                                          <  6:104148/126‑1 (MQ=255)                                                                                tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                   >  5:129451/1‑149 (MQ=255)                                                                                tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                   >  2:204213/1‑149 (MQ=255)                                                                                  cAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                   >  4:224769/1‑147 (MQ=255)                                                                                        ttttATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAAc                                                           >  6:282352/1‑149 (MQ=255)                                                                                                aGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCGGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTAtt                                                    >  6:260033/1‑148 (MQ=255)                                                                                                   aaCTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAAt                                                <  4:64393/149‑1 (MQ=255)                                                                                                             cagcCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAgg                                      <  1:204213/149‑1 (MQ=255)                                                                                                              agcCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAggg                                     >  1:131685/1‑149 (MQ=255)                                                                                                                 cAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTAtt                                                    >  3:69642/1‑131 (MQ=255)                                                                                                                 cAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTAtt                                                    <  4:69642/131‑1 (MQ=255)                                                                                                                 cAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTc                                   <  1:128381/148‑1 (MQ=255)                                                                                                                   gAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATc                                               <  1:323068/134‑1 (MQ=255)                                                                                                                   gAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATc                                               >  2:323068/1‑134 (MQ=255)                                                                                                                     aCCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTc                              <  6:129451/149‑1 (MQ=255)                                                                                                                        cACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAg                           <  2:285822/149‑1 (MQ=255)                                                                                                                                agagCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgc                   <  3:146201/149‑1 (MQ=255)                                                                                                                                 gagCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGt                               <  8:21853/136‑1 (MQ=255)                                                                                                                                 gagCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGt                               >  7:21853/1‑136 (MQ=255)                                                                                                                                  agCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgcgc                 <  3:224769/149‑1 (MQ=255)                                                                                                                                    cgcgcgGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgcgcgc               >  4:168162/1‑149 (MQ=255)                                                                                                                                      cgcgGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCt             <  2:131685/149‑1 (MQ=255)                                                                                                                                        tagcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGa                                                  <  6:124040/108‑1 (MQ=255)                                                                                                                                        tagcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGa                                                  >  5:124040/3‑110 (MQ=255)                                                                                                                                        tagcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCttt           >  1:126716/3‑149 (MQ=255)                                                                                                                                        cggcgACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCttt           <  6:111893/149‑1 (MQ=255)                                                                                                                                          gcgACTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         >  4:174578/1‑149 (MQ=255)                                                                                                                                          gcgACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         <  5:282352/149‑1 (MQ=255)                                                                                                                                           cgACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         >  4:8921/1‑148 (MQ=255)                                                                                                                                                 aaaCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGt  >  8:161005/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGACTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGT  >  NC_002947/5909151‑5909443

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