Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,010,969	A→C	5.9%	V10V (GTA→GTC)	bcsG →	endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,010,969	0	A	C	5.9%	215.0 / 3.9	68	V10V (GTA→GTC)	bcsG	endoglucanase
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base C (2/2);  total (32/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 = 8.89e-01

GGCCTTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGG  >  NC_002947/3010821‑3011113                                                                                                                                                      |                                                                                                                                                  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                   <  3:207477/138‑1 (MQ=255)                                                                                                                                            cccTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTcc                   >  4:207477/1‑138 (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)                                                                                                                                                      |                                                                                                                                                  GGCCTTGTTGCTAGTGCGCCTGTGTTCGGCATTACGTGGCGTGTTGCGCCGTCGCCTT‑CCCCCCCGTTACCTGAAACCGCTGGGCGTACGCCGCCGCGCTGCACGCCAGGAACCCCAGAATGACTAATACCACCATCACCCTCCCGGTACTGGCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCG‑ACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTGG  >  NC_002947/3010821‑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: ‑