Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	901,737	C→T	6.2%	G42G (GGC→GGT)	trxB →	thioredoxin reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	901,737	0	C	T	6.2%	187.8 / 2.0	64	G42G (GGC→GGT)	trxB	thioredoxin reductase
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base T (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

CCGCTCGACAGGAGTTCGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAACGCTGTCGACCTGGC  >  NC_002947/901592‑901883                                                                                                                                                  |                                                                                                                                                     ccGCTCGACAGGAGTTCGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAg                                                                                                                                                   <  1:273511/149‑1 (MQ=255)   gCTCGACAGGAGTTCGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCt                                                                                                                                                 <  2:190750/149‑1 (MQ=255)       gACAGGAGTTCGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGacc                                                                                                                                             >  1:189763/1‑149 (MQ=255)          aGGAGTTCGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGaccacc                                                                                                                                          >  2:349307/1‑149 (MQ=255)               ttCGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGaccaccaccac                                                                                                                                     <  1:349307/149‑1 (MQ=255)                 cGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGaccaccaccacc                                                                                                                                    >  7:205025/1‑148 (MQ=255)                    cATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAg                                                                                                                                >  2:231490/1‑149 (MQ=255)                        tCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGa                                                                                                                            <  8:205025/149‑1 (MQ=255)                          tGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGTAGGCGGGTGGCCAGCTGACCACCACCACCGaa                                                                                                                                 >  6:135823/1‑142 (MQ=255)                          tGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGTAGGCGGGTGGCCAGCTGACCACCACCACCGaa                                                                                                                                 <  5:135823/142‑1 (MQ=255)                          tGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGAc                                                                                                                           >  4:63267/1‑148 (MQ=255)                          tGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCACTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACAACCACCGAAGTCGAc                                                                                                                           >  8:82842/1‑148 (MQ=255)                           gAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTc                                                                                                                              >  4:50675/1‑144 (MQ=255)                           gAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTc                                                                                                                              <  3:50675/144‑1 (MQ=255)                           gAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACaa                                                                                                                         >  8:185971/1‑149 (MQ=255)                            aaGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGAc                                                                                                                           <  5:85389/146‑1 (MQ=255)                             aGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGAATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACaa                                                                                                                         >  8:241204/1‑147 (MQ=255)                              gTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTg                                                                                                                      <  1:163863/149‑1 (MQ=255)                                     aTTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCggg                                                                                                                >  7:89359/1‑148 (MQ=255)                                     aTTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCggg                                                                                                                <  8:89359/148‑1 (MQ=255)                                         gcgcgTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGTCAGCTGgccacccgc                                                                                                                                       >  7:87860/1‑118 (MQ=255)                                         gcgcgTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGTCAGCTGgccacccgc                                                                                                                                       <  8:87860/121‑4 (MQ=255)                                         gcgcgTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGTCAGCTGgccaccagc                                                                                                                                       <  7:83321/121‑3 (MQ=255)                                         gcgcgTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGTCAGCTGgccaccagc                                                                                                                                       >  8:83321/1‑119 (MQ=255)                                         gcgcgTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACt                                                                                                                       >  6:91018/1‑137 (MQ=255)                                         gcgcgTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACt                                                                                                                       <  5:91018/137‑1 (MQ=255)                                          cgcgTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGAccc                                                                                                          <  2:92908/149‑1 (MQ=255)                                                atcatTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACgg                                                                                                    >  2:169315/1‑149 (MQ=255)                                                 tcatTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGc                                                                                                   >  2:180601/1‑149 (MQ=255)                                                 tcatTCTCGGTTCCGG‑CCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGcc                                                                                                  >  2:70117/1‑149 (MQ=255)                                                   atTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCt                                                                                                 >  3:82505/1‑149 (MQ=255)                                                    ttCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGAGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGAccc                                                                                                          <  8:132505/139‑1 (MQ=255)                                                    ttCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGAGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGAccc                                                                                                          >  7:132505/1‑139 (MQ=255)                                                       tcGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACgg                                                                                                    >  5:246452/1‑142 (MQ=255)                                                       tcGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACgg                                                                                                    <  6:246452/142‑1 (MQ=255)                                                                     gCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGAt                                                                             >  2:80309/1‑149 (MQ=255)                                                                             cAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGt                                                                      <  1:231490/148‑1 (MQ=255)                                                                                cgcTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGc                                                                  <  2:189763/149‑1 (MQ=255)                                                                                   tGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGc                                                                  >  8:110302/1‑146 (MQ=255)                                                                                      ggTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGc                                                            <  7:185971/149‑1 (MQ=255)                                                                                                  ccGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGtt                                                >  2:283457/1‑149 (MQ=255)                                                                                                    gTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCGCGCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGttt                                               <  1:70117/149‑1 (MQ=255)                                                                                                    gTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGtttt                                              <  3:63267/149‑1 (MQ=255)                                                                                                    gTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGtttt                                              <  1:180601/149‑1 (MQ=255)                                                                                                      cccAACCTCCAGCCGCTGCTGACCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGa                                            <  8:332220/148‑1 (MQ=255)                                                                                                       ccAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGa                                            >  6:58438/1‑148 (MQ=255)                                                                                                        cAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGaaa                                          >  1:108505/1‑149 (MQ=255)                                                                                                         aaCCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAAc                                         >  7:341503/1‑149 (MQ=255)                                                                                                            cTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACAACCACCACCGAAGTCGACAACTGGCCGGGCGCCCCCCACGGCCTGACCG‑GCCC‑CGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCg                                       >  2:188454/1‑148 (MQ=255)                                                                                                              cAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGaaa                                    >  8:129222/1‑149 (MQ=255)                                                                                                                aGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATc                                  <  5:163382/149‑1 (MQ=255)                                                                                                                        ctgATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GTCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTc                            <  7:241204/147‑1 (MQ=255)                                                                                                                        ctgATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTc                            <  4:82505/147‑1 (MQ=255)                                                                                                                        ctgATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGa                          <  1:80309/149‑1 (MQ=255)                                                                                                                         tgATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGAc                         >  8:195295/1‑149 (MQ=255)                                                                                                                           aTCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGAc                         <  1:169315/147‑1 (MQ=255)                                                                                                                              aCCGGCATGCATGCGGGTTGCCATGTGTACACCACCAACGAAGTTGACAACTTGCCGGGGGACCCCCAAGGCCTTTCCG‑GCCCGGGTTCAGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACcaca                     <  1:283457/149‑1 (MQ=255)                                                                                                                              aCCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACcaca                     >  8:89096/1‑148 (MQ=255)                                                                                                                               ccGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGACGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACAt                    >  1:249398/1‑148 (MQ=255)                                                                                                                                  gCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGCCCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAAc                >  4:316226/1‑149 (MQ=255)                                                                                                                                       cAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAACGCTGt           <  7:129222/149‑1 (MQ=255)                                                                                                                                         tggtggtggtCAGCTGGCCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAACGCTGTc          >  5:274498/5‑148 (MQ=255)                                                                                                                                                ggCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGGCCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAACGCTTTCGACCTGGc  >  8:245770/1‑149 (MQ=255)                                                                                                                                                 gCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAACGCTGTCGACCTGGc  <  5:58438/148‑1 (MQ=255)                                                                                                                                                 gCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAACGCTGTCGACCTGGc  >  1:322569/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCGCTCGACAGGAGTTCGACATGTCTGAAGTACGTCATTCGCGCGTCATCATTCTCGGTTCCGGCCCTGCCGGTTACAGCGCTGCGGTGTATGCCGCCCGTGCCAACCTCAAGCCGCTGCTGATCACCGGCATGCAGGCGGGTGGCCAGCTGACCACCACCACCGAAGTCGACAACTGGCCGGGCGACCCCCACGGCCTGACCG‑GCCC‑GGCTCTGATGCAGCGTATGCAAGAGCACGCCGAGCGTTTTGAAACCGAAATCGTCTTCGACCACATCAACGCTGTCGACCTGGC  >  NC_002947/901592‑901883

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