Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,012,112	A→T	7.9%	D377V (GAC→GTC)	bcsG →	endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,012,112	0	A	T	7.9%	233.6 / 9.3	76	D377V (GAC→GTC)	bcsG	endoglucanase
Reads supporting (aligned to +/- strand):  ref base A (39/31);  new base T (3/3);  total (42/34)
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.66e-01

GTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCATGGTGATC  >  NC_002947/3011967‑3012258                                                                                                                                                  |                                                                                                                                                   gtggGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGAcggc                                                                                                                                                >  1:1214526/1‑150 (MQ=255)  tggGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGAcggcg                                                                                                                                               <  1:406530/150‑1 (MQ=255)      cTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGTCGGCGGtga                                                                                                                                           <  1:119474/150‑3 (MQ=255)        tCGACGGCTCTCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTAccc                                                                                                                                         <  2:330566/150‑1 (MQ=255)           aCGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCTGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCcgcg                                                                                                                                      >  2:822326/1‑150 (MQ=255)                tCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGTCGGCGGtgac                                                                                                                                          <  2:693630/141‑4 (MQ=255)                tCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGTCGGCGGtgac                                                                                                                                          >  1:693630/1‑138 (MQ=255)                   ccGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGAc                                                                                                                              <  1:1395349/150‑1 (MQ=255)                          ccAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCATATTGTACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCa                                                                                                                       >  1:1327503/1‑150 (MQ=255)                                  cTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGcc                                                                                                               >  1:297017/1‑150 (MQ=255)                                       tGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACTCTGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCAGGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATcc                                                                                                          >  2:597385/1‑149 (MQ=255)                                              cgcCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGttt                                                                                                   <  1:1343366/150‑1 (MQ=255)                                              cgcCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGttt                                                                                                   >  1:296929/1‑150 (MQ=255)                                               gcCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTg                                                                                                  <  2:1327503/150‑1 (MQ=255)                                               gcCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTAGCGTCACCGCCGACGGCGGTACTCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTg                                                                                                  >  1:1314347/1‑150 (MQ=255)                                                ccGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGGGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTgc                                                                                                 <  2:973314/150‑2 (MQ=255)                                                 cGTTGGTGGGGGCAACGCAAAGGCACCCCGGAAGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGc                                                                                                <  2:336771/150‑1 (MQ=255)                                                  gTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCg                                                                                               >  2:1320771/1‑150 (MQ=255)                                                     ggtggGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAAc                                                                                            <  2:1214526/150‑1 (MQ=255)                                                       tggGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGt                                                                                                                                     >  1:1403165/1‑107 (MQ=255)                                                       tggGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGt                                                                                                                                     <  2:1403165/107‑1 (MQ=255)                                                       tggGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAAt                                                                                                            <  2:508352/132‑1 (MQ=255)                                                       tggGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAAt                                                                                                            >  1:508352/1‑132 (MQ=255)                                                       tggGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACtt                                                                                          >  2:29499/1‑150 (MQ=255)                                                         ggggCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGGTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTgg                                                                                        >  1:58435/1‑150 (MQ=255)                                                           ggCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTg                                                                                                      <  2:419234/134‑1 (MQ=255)                                                           ggCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTg                                                                                                      >  1:419234/1‑134 (MQ=255)                                                              aaCGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCAcct                                                                                   >  1:500266/1‑150 (MQ=255)                                                                     aGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGa                                                                            >  1:45602/1‑150 (MQ=255)                                                                           cccGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGAcggc                                                                                                                                                <  2:546212/76‑1 (MQ=255)                                                                           cccGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGAcggc                                                                                                                                                >  1:546212/1‑76 (MQ=255)                                                                                 cGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAggg                                                                <  1:645690/150‑1 (MQ=255)                                                                                   aaCATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCa                                                              <  1:29499/150‑1 (MQ=255)                                                                                   aaCATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAAATTGGCACCTTCCTTGACGAACTGGAGGGCa                                                              <  1:1320771/150‑1 (MQ=255)                                                                                     cATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGTCGGCGGtgac                                                                                                                                          <  1:1249165/72‑4 (MQ=255)                                                                                     cATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGTCGGCGGtgac                                                                                                                                          >  2:1249165/1‑69 (MQ=255)                                                                                      aTGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTg                                                           <  2:45602/150‑1 (MQ=255)                                                                                           tCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGgccgc                                                      >  1:1309158/1‑150 (MQ=255)                                                                                           tCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGACTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGgccgc                                                      >  1:312437/1‑150 (MQ=255)                                                                                            cACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGgccgcc                                                     <  2:296929/150‑1 (MQ=255)                                                                                             aCTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTg                                                                                         <  2:203732/113‑1 (MQ=255)                                                                                             aCTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTg                                                                                         >  1:203732/1‑113 (MQ=255)                                                                                             aCTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGccgccg                                                    >  2:959222/1‑150 (MQ=255)                                                                                             aCTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGccgccg                                                    >  2:1316908/1‑150 (MQ=255)                                                                                              cTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCcgccgc                                                   <  1:416027/150‑1 (MQ=255)                                                                                               tGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATc                                                                                                           <  1:1076477/93‑1 (MQ=255)                                                                                               tGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATc                                                                                                           >  2:1076477/1‑93 (MQ=255)                                                                                               tGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGa                                                                   <  2:651563/133‑1 (MQ=255)                                                                                               tGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGa                                                                   >  1:651563/1‑133 (MQ=255)                                                                                               tGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCcgcg                                                  >  1:124362/1‑150 (MQ=255)                                                                                               tGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCcgcg                                                  >  2:1255526/1‑150 (MQ=255)                                                                                                   cTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAgt                                              <  1:1316908/150‑1 (MQ=255)                                                                                                      caacaGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAgtggt                                           >  1:609991/1‑150 (MQ=255)                                                                                                       aacaGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAgtggtg                                          >  1:767822/1‑150 (MQ=255)                                                                                                         caGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGgcc                                                        >  1:35811/1‑134 (MQ=255)                                                                                                         caGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGgcc                                                        <  2:35811/134‑1 (MQ=255)                                                                                                              tCAGCCTGCATGACGGCAACCGTATCGTCACCGCTGACGGCGGTACCCGCGTCGCCGACTATGCCACCCGCGCCAATCGCCTGCTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTg                                   >  2:738753/1‑150 (MQ=255)                                                                                                              tCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTg                                   >  2:452320/1‑150 (MQ=255)                                                                                                                aGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTc                                 >  1:69729/1‑150 (MQ=255)                                                                                                                  ccTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTccc                               >  1:1383731/1‑150 (MQ=255)                                                                                                                        tGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTcc                                <  2:419706/143‑1 (MQ=255)                                                                                                                        tGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTcc                                >  1:419706/1‑143 (MQ=255)                                                                                                                        tGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCa                         >  2:1216615/1‑150 (MQ=255)                                                                                                                             gCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATgggg                    >  1:489882/1‑150 (MQ=255)                                                                                                                                       cGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCa          <  2:297017/150‑1 (MQ=255)                                                                                                                                       cGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCa          <  2:500266/150‑1 (MQ=255)                                                                                                                                       cGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCa          >  1:743131/1‑150 (MQ=255)                                                                                                                                        gTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGtt             <  1:1346688/146‑1 (MQ=255)                                                                                                                                        gTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGtt             >  2:1346688/1‑146 (MQ=255)                                                                                                                                        gTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCAt         <  1:452320/150‑1 (MQ=255)                                                                                                                                         tCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGt                                  >  2:429979/1‑124 (MQ=255)                                                                                                                                         tCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGt                                  <  1:429979/124‑1 (MQ=255)                                                                                                                                          taCCGCCGTCGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCATgg       >  2:507816/2‑150 (MQ=255)                                                                                                                                           aCCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCATGGt      <  1:848788/150‑1 (MQ=255)                                                                                                                                             cgccgACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTCACGAACTGGAGGGCAGTGGCCGCCGCGTCGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCATGGTGa    <  2:1314347/150‑1 (MQ=255)                                                                                                                                               ccgACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCATGGTGATc  >  1:530222/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GTGGGCTTCGACGGCTCGCCGATTGCCAGTGACCTGGATGTCTTGCGCCGTTGGTGGGGGCAACGCAAAGGCACCCCGGACGAACATGTGTCACTGTTCTACAACAGCATCAGCCTGCATGACGGCAACCGTATCGTCACCGCCGACGGCGGTACCCGCGTGGCCGACTATGCCACCCGCGCCAATCGCCTGTTTGGCGAACTTGGCACCTTCCTTGACGAACTGGAGGGCAGTGGCCGCCGCGTAGTGGTGGCCATTGTCCCTGAGCATGGGGCCGCGTTGCATGGTGATC  >  NC_002947/3011967‑3012258

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