Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,786,991	T→G	5.3%	L135R (CTG→CGG)	trxC →	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,786,991	0	T	G	5.3%	237.5 / 3.4	75	L135R (CTG→CGG)	trxC	thioredoxin
Reads supporting (aligned to +/- strand):  ref base T (40/31);  new base G (2/2);  total (42/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 = 9.29e-01

TGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGC  >  NC_002947/5786852‑5787137                                                                                                                                            |                                                                                                                                                   tGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTg                                                                                                                                           >  5:223900/1‑149 (MQ=255)   ccgccTGGACAAGCTCGACATCTAAGCCAACCGTAACCTCGCAGGTCAACTGGTCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTgg                                                                                                                                          <  8:136939/148‑1 (MQ=255)      ccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGc                                                                                                                                      >  4:274126/1‑149 (MQ=255)          gCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTa                                                                                                                                   <  6:223900/148‑1 (MQ=255)           ccAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGt                                                                                                                                            >  5:65269/1‑138 (MQ=255)           ccAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGt                                                                                                                                            <  6:65269/138‑1 (MQ=255)            cAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGc                                                                                                                                 >  3:334964/1‑148 (MQ=255)            cAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGc                                                                                                                                 >  6:278844/1‑148 (MQ=255)               gCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAg                                                                                                                             <  5:278844/149‑1 (MQ=255)                 tCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAggg                                                                                                                           >  3:211707/1‑149 (MQ=255)                         gAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGc                                                                                                                   >  8:253305/1‑149 (MQ=255)                            gCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGtt                                                                                                                <  7:253305/149‑1 (MQ=255)                             ccAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGttt                                                                                                                       >  3:199965/1‑141 (MQ=255)                             ccAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGttt                                                                                                                       <  4:199965/141‑1 (MQ=255)                               aaCCGTAACCTCGCAGGGCAACTGGGCATAGGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttc                                                                                                             >  1:32477/1‑149 (MQ=255)                                aCCGTAACCTCGCAGGTCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttct                                                                                                            <  2:32477/149‑1 (MQ=255)                                 ccGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTaa                                                                                                                     >  3:16501/1‑139 (MQ=255)                                 ccGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTaa                                                                                                                     <  4:16501/139‑1 (MQ=255)                                        cTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAagcag                                                                                                    <  4:211707/149‑1 (MQ=255)                                            cAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGc                                                                                                                   >  8:199903/1‑130 (MQ=255)                                            cAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGc                                                                                                                   <  7:199903/130‑1 (MQ=255)                                            cAGGGACACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTGCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGaa                                                                                                <  8:140968/149‑1 (MQ=255)                                               ggCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAagcagc                                                                                                   >  3:112198/1‑143 (MQ=255)                                               ggCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAagcagc                                                                                                   <  4:112198/143‑1 (MQ=255)                                               ggCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGc                                                                                             <  4:334964/149‑1 (MQ=255)                                                  aaCTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGc                                                                                          >  5:38847/1‑149 (MQ=255)                                                  aaCTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGc                                                                                          >  2:183152/1‑149 (MQ=255)                                                   aCTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCt                                                                                         >  3:218049/1‑149 (MQ=255)                                                          aTACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACaaa                                                                                  >  4:247593/1‑149 (MQ=255)                                                              gCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgc                                                                              >  3:340220/1‑149 (MQ=255)                                                                   aTCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGt                                                                         >  8:70818/1‑149 (MQ=255)                                                                      cccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGt                                                                         >  6:345589/1‑146 (MQ=255)                                                                      cccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGt                                                                         <  5:345589/146‑1 (MQ=255)                                                                       ccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGc                                                                     >  2:114177/1‑149 (MQ=255)                                                                         aGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCtt                                                                   <  7:70818/149‑1 (MQ=255)                                                                             tgctgctgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGATCGACAAACTGCTGAGTCAGCTTGTGc                                                               >  4:54611/1‑149 (MQ=255)                                                                                  ctgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGc                                                                                             <  3:197669/114‑1 (MQ=255)                                                                                  ctgTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGc                                                                                             >  4:197669/1‑114 (MQ=255)                                                                                    gTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCt                                                        >  5:348557/1‑149 (MQ=255)                                                                                     ttCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTc                                                       <  1:114177/149‑1 (MQ=255)                                                                                      tCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTc                                                       <  3:54611/148‑1 (MQ=255)                                                                                          gAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGt                                                  >  2:168345/1‑149 (MQ=255)                                                                                          gAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGt                                                  <  3:247593/149‑1 (MQ=255)                                                                                          ccaCGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGt                                                  <  1:95142/147‑1 (MQ=255)                                                                                               gCCGCGAAATCAGCCGCCAGGCCGGGGCGCTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATc                                             <  5:316738/149‑1 (MQ=255)                                                                                                ccGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAAt                                              >  5:311325/1‑147 (MQ=255)                                                                                                   cgAAATCCGCCGCCAGGCCGGGGCGTTTCCACTCAAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCg                                         >  2:74797/1‑149 (MQ=255)                                                                                                    gAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                        >  1:6320/1‑149 (MQ=255)                                                                                                    gAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCgg                                        >  1:717/1‑149 (MQ=255)                                                                                                          aGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTTCCGCGGCTCCAGGt                                                  >  6:38484/1‑133 (MQ=255)                                                                                                           gccgccAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTTGACAAAct                                                                                >  1:224133/1‑102 (MQ=255)                                                                                                           gccgccAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTTGACAAAct                                                                                <  2:224133/102‑1 (MQ=255)                                                                                                           gccgccAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTg                                 <  3:274126/149‑1 (MQ=255)                                                                                                           gccgccAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTg                                 >  7:251730/1‑149 (MQ=255)                                                                                                               ccAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgt                             >  8:58066/1‑149 (MQ=255)                                                                                                                 aGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtgg                           <  7:58066/149‑1 (MQ=255)                                                                                                                 aGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtgg                           >  4:118229/1‑149 (MQ=255)                                                                                                                  ggCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCgtggt                          <  4:340220/149‑1 (MQ=255)                                                                                                                                cACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtt            <  4:218049/149‑1 (MQ=255)                                                                                                                                cACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAtt            >  3:89639/1‑149 (MQ=255)                                                                                                                                cACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGAAGGCTGGCGTGGTGCGATTCGCGCAtt            <  6:38847/149‑1 (MQ=255)                                                                                                                                 aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGc                               >  3:24311/1‑129 (MQ=255)                                                                                                                                 aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGc                               <  4:24311/129‑1 (MQ=255)                                                                                                                                 aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAttt           >  6:8058/1‑149 (MQ=255)                                                                                                                                 aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAttt           >  4:32427/1‑149 (MQ=255)                                                                                                                                 aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAACGGCAGGCTGGCGTGGTGCGATTCGCGCAttt           >  7:40404/1‑149 (MQ=255)                                                                                                                                 aCTCCAGCAGCGACTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgc                 >  7:269601/1‑143 (MQ=255)                                                                                                                                 aCTCCAGCAGCGACTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgc                 <  8:269601/143‑1 (MQ=255)                                                                                                                                     cAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAAt                                              >  6:176852/1‑110 (MQ=255)                                                                                                                                     cAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAAt                                              <  5:176852/110‑1 (MQ=255)                                                                                                                                      aGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCAcc      <  6:348557/149‑1 (MQ=255)                                                                                                                                       gTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGa                      >  5:235429/1‑132 (MQ=255)                                                                                                                                       gTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGa                      >  5:230172/1‑132 (MQ=255)                                                                                                                                       gTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGa                      <  6:235429/132‑1 (MQ=255)                                                                                                                                       gTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGa                      <  6:230172/132‑1 (MQ=255)                                                                                                                                          gctgctGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGc  <  4:89639/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                   TGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTTTCACCGAGC  >  NC_002947/5786852‑5787137

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: ‑