Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,010,966	G→A	5.8%	P9P (CCG→CCA)	bcsG →	endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,010,966	0	G	A	5.8%	205.2 / 2.0	69	P9P (CCG→CCA)	bcsG	endoglucanase
Reads supporting (aligned to +/- strand):  ref base G (31/34);  new base A (2/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

CCTGGCCTTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGG  >  NC_002947/3010818‑3011113                                                                                                                                                      |                                                                                                                                                     cctggcctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCg                                                                                                                                                      <  5:328498/149‑1 (MQ=255) cctggcctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCg                                                                                                                                                      >  2:318890/1‑149 (MQ=255) cctggcctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCg                                                                                                                                                      >  7:316677/1‑149 (MQ=255)   tggcctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCgg                                                                                                                                                     >  3:312584/1‑148 (MQ=255)    ggcctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTa                                                                                                                                                   >  7:82933/1‑149 (MQ=255)     gcctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTAc                                                                                                                                                  <  8:45452/149‑1 (MQ=255)       ctTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTg                                                                                                                                                <  7:120490/149‑1 (MQ=255)       ctTGTTGCTAGTGCGCCTGTGTCCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTg                                                                                                                                                >  5:380863/1‑149 (MQ=255)        ttgttgCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTgg                                                                                                                                               >  2:297477/1‑149 (MQ=255)          gttgCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGc                                                                                                                                              <  3:247960/148‑1 (MQ=255)          gttgCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGcc                                                                                                                                             >  5:356950/1‑149 (MQ=255)             gCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGc                                                                                                                                          >  7:96236/1‑149 (MQ=255)                 gTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGc                                                                                                                                      <  6:356950/149‑1 (MQ=255)                       cTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGAc                                                                                                                                >  6:70330/1‑149 (MQ=255)                         gtgtTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCCCTGGGCGTACGCCCCCGCGCAGCACGCCAGGAACCCCAGAATGACTAATACCACCATAACCCGCCCGGTACAGGCCCGCGGGCCGGGCCTg                                                                                                                              >  2:83659/1‑149 (MQ=255)                          tgtTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTTCCCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCAGTCCAGGCCAGCGGGccagtacc                                                                                                                               <  2:4226/148‑7 (MQ=255)                              cGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAgg                                                                                                                         >  3:3085/1‑149 (MQ=255)                              cGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAg                                                                                                                          <  7:303967/148‑1 (MQ=255)                               ggCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAgg                                                                                                                         <  4:312584/148‑1 (MQ=255)                                gCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGt                                                                                                                       <  2:214576/149‑1 (MQ=255)                                gCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGt                                                                                                                       <  1:297477/149‑1 (MQ=255)                                     aCGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGaa                                                                                                                   <  1:318890/148‑1 (MQ=255)                                      cGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAAc                                                                                                                  <  6:380863/148‑1 (MQ=255)                                        tGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCCCGCCAGGACCCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCGGCTGGCGTGGAGTGGGAGGGTGGAACgtc                                                                                                               >  5:304387/1‑149 (MQ=255)                                              gtTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCCCCATCACCCTCCCGGTACTGGCCCGCTGGCCGGGACTGGGAGGGTGGAACCTCTAtttt                                                                                                         >  4:233537/1‑149 (MQ=255)                                                  cgcCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttttt                                                                                                      <  1:83183/148‑1 (MQ=255)                                                  cgcCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTAttttttt                                                                                                      <  4:3085/148‑1 (MQ=255)                                                      gTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCaaa                                                                                                 <  6:167173/149‑1 (MQ=255)                                                         gCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGt                                                                                               <  5:347839/148‑1 (MQ=255)                                                          ccTT‑CCCCCCCGTTACCTGAAACCGCTGTGCGTACGCCGCAGCTCTGCACGCCAGGAACCGCAGAATGACTAATACCAACATCACCCTCCCGGTACTGGCGCGCTGGCCTTGACTGGGAGTGTGGAACCTCTATTTTTTTGCAAAGttt                                                                                             <  8:316677/149‑1 (MQ=255)                                                           cTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTc                                                                                            <  8:82933/149‑1 (MQ=255)                                                            tt‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTct                                                                                           >  2:102136/1‑149 (MQ=255)                                                            tt‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAAACTCTATTTTTTTGCAAAGTTTct                                                                                           >  3:20060/1‑149 (MQ=255)                                                                ccccccGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTctgc                                                                                         <  5:70330/148‑1 (MQ=255)                                                                     cGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTcg                                                                                   >  3:206217/1‑149 (MQ=255)                                                                         aCCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCt                                                                               >  2:99966/1‑149 (MQ=255)                                                                          ccTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTg                                                                              >  8:40850/1‑149 (MQ=255)                                                                           cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                             <  8:96236/149‑1 (MQ=255)                                                                           cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                             <  1:83659/149‑1 (MQ=255)                                                                           cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                             >  3:161681/1‑149 (MQ=255)                                                                           cTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                             <  3:233537/149‑1 (MQ=255)                                                                            tGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTggg                                                                            >  2:198766/1‑149 (MQ=255)                                                                             gAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                             >  8:190302/1‑147 (MQ=255)                                                                             gAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTgg                                                                             <  7:190302/147‑1 (MQ=255)                                                                                aCCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTc                                                                                    >  5:142181/1‑137 (MQ=255)                                                                                aCCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTc                                                                                    <  6:142181/137‑1 (MQ=255)                                                                                  cGCTGGGCGTACGCCGGCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCt                                                                      <  6:304387/149‑1 (MQ=255)                                                                                     tGGGCGTACGCCGCCGCGCTGCACGCCAGGCACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGGAAAGTTTCTGCTGGTCGCGCTGGGCATGCTc                                                                     >  3:296094/1‑147 (MQ=255)                                                                                      gggCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑Ac                                                                 >  1:194731/1‑149 (MQ=255)                                                                                                       cTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGaa                                                <  4:206217/149‑1 (MQ=255)                                                                                                          cACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAAc                                               <  1:238076/147‑1 (MQ=255)                                                                                                               cAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCtgttgt                                        >  8:186427/1‑149 (MQ=255)                                                                                                                aGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTgttgtt                                       <  1:99966/149‑1 (MQ=255)                                                                                                                  gAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCg                                     <  1:198766/149‑1 (MQ=255)                                                                                                                     ccccAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTg                                  <  7:205375/149‑1 (MQ=255)                                                                                                                        cAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGt                               <  4:20060/149‑1 (MQ=255)                                                                                                                               aCTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGc                         <  3:86371/148‑1 (MQ=255)                                                                                                                                 tAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTg                      >  1:350605/1‑149 (MQ=255)                                                                                                                                   aTACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACATGATGTTCGCTGCGTTTCTGATTGt                     >  8:165925/1‑148 (MQ=255)                                                                                                                                    tACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTcc                   <  4:357397/149‑1 (MQ=255)                                                                                                                                      ccaccaTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCt                 >  7:127578/1‑149 (MQ=255)                                                                                                                                       caccaTCACCCTCCCGCTACTGGCCCGCTGGCCTCGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCAGGCTAG‑ACTTCCAGGCGCTGCCGAACAGGTTGTTCGCTGCGTTTCTGCTTGTCCCtt                >  5:277469/1‑149 (MQ=255)                                                                                                                                          caTCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGc             >  1:299333/1‑149 (MQ=255)                                                                                                                                            tCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGTACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGcc            <  1:102136/149‑1 (MQ=255)                                                                                                                                              aCCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCggg         >  2:230743/1‑149 (MQ=255)                                                                                                                                              aCCCTCCCAGTCCAGGCCAGCGGGCCAGTACCGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCggg         >  6:112098/1‑149 (MQ=38)                                                                                                                                               cccTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTcc                   >  4:207477/1‑138 (MQ=255)                                                                                                                                               cccTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTcc                   <  3:207477/138‑1 (MQ=255)                                                                                                                                               cccTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCggg         >  2:215556/1‑148 (MQ=255)                                                                                                                                                 cTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCaa      <  4:161681/149‑1 (MQ=255)                                                                                                                                                 cTCCCGGTACAGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTATCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGATGCGTTTCTGCTTGTCCCTTTGCCGGGc        >  6:144238/1‑147 (MQ=255)                                                                                                                                                  tCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAg     >  5:25120/1‑149 (MQ=255)                                                                                                                                                     cGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTgg  <  6:25120/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCTGGCCTTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGG  >  NC_002947/3010818‑3011113

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