Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,842,798	G→A	6.2%	L680L (CTG→TTG)	gacS ←	sensor protein GacS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,842,798	0	G	A	6.2%	188.7 / 2.4	64	L680L (CTG→TTG)	gacS	sensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/2);  total (33/31)
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.12e-01

CAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGG  >  NC_002947/1842650‑1842946                                                                                                                                                     |                                                                                                                                                     cAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAgg                                                                                                                                                     <  1:591177/150‑1 (MQ=255)     cgcgGCCGTTCATGGGGAGCAGCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGTCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTgg                                                                                                                                                 <  2:815254/150‑1 (MQ=255)          ccGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCa                                                                                                                                            >  1:653071/1‑150 (MQ=255)            gTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCAtt                                                                                                                                          >  1:1820536/1‑150 (MQ=255)             tCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTg                                                                                                                                         >  2:1879100/1‑150 (MQ=255)               cATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTgtc                                                                                                                                       <  2:170821/150‑1 (MQ=255)                    cGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGa                                                                                                                                  >  2:42633/1‑150 (MQ=255)                          tCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCaga                                                                                                                            >  1:1216157/1‑150 (MQ=255)                            tGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAAGTTGGCGGCATTGTCGTCGACGCAGAGt                                                                                                                          >  1:1424862/1‑150 (MQ=255)                            tGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAAGTTGGCGGCATTGTCGTCGACGCAGAGt                                                                                                                          >  1:1425250/1‑150 (MQ=255)                                 atccatcAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTAtttt                                                                                                                     >  1:684812/1‑150 (MQ=255)                                   ccatcAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCACCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCg                                                                                                                   <  1:1189485/150‑1 (MQ=255)                                      tcAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtc                                                                                                                >  1:70985/1‑150 (MQ=255)                                      tcAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGgtc                                                                                                                <  1:634332/150‑1 (MQ=255)                                              aGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAAGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTgg                                                                                                        <  2:1425250/150‑1 (MQ=255)                                              aGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAAGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTgg                                                                                                        <  2:1424862/150‑1 (MQ=255)                                                    aaCAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTg                                                                                                  <  2:1668288/150‑1 (MQ=255)                                                    aaCAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTg                                                                                                  <  2:1457012/150‑1 (MQ=255)                                                    aaCAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTg                                                                                                  >  1:357447/1‑150 (MQ=255)                                                      cAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCg                                                                                                                   <  2:303311/131‑1 (MQ=255)                                                      cAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCg                                                                                                                   >  1:303311/1‑131 (MQ=255)                                                       aaTTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTgggg                                                                                               <  1:1283952/150‑1 (MQ=255)                                                           cttcCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGACTTCGAGCAGGGCTTTTACCAGCAGCAGGGTGGCGGCATTGTCGTCGACGCACAATATTATAGCTCGGCGCTCGCCGTATGGTTTGc                                                                                           >  2:873990/1‑150 (MQ=255)                                                            ttcCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTgtc                                                                                                                                       <  1:617062/105‑1 (MQ=255)                                                            ttcCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTgtc                                                                                                                                       >  2:617062/1‑105 (MQ=255)                                                                    cGGCCTGTACCGTCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTACCCTCgg                                                                                  <  2:514889/150‑1 (MQ=255)                                                                          gattcgcctcATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCATTGGGGTTGCTCCCCTCGGCCAGCg                                                                            >  1:393327/5‑150 (MQ=255)                                                                            aCCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGa                                                                          <  2:70985/150‑1 (MQ=255)                                                                               gccgcATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACgg                                                                       >  1:1247554/1‑150 (MQ=255)                                                                                 cgcATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCg                                                                     >  2:1368550/1‑150 (MQ=255)                                                                                   cATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAg                                                                   <  1:997293/150‑1 (MQ=255)                                                                                      aaCCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGtt                                                                <  2:158505/150‑1 (MQ=255)                                                                                       aCCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTg                                                               <  2:1577461/150‑1 (MQ=255)                                                                                        ccATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGc                                                              <  1:1879100/150‑1 (MQ=255)                                                                                         cATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCa                                                             >  1:9348/1‑150 (MQ=255)                                                                                         cATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGACGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCa                                                             >  1:1259033/1‑150 (MQ=255)                                                                                          aTTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCAc                                                            <  2:393787/150‑1 (MQ=255)                                                                                             gctgACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAg                                                         >  1:1301520/3‑150 (MQ=255)                                                                                              ttgACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGATTGCTACTCTAGGCAAGCGGGCGGCGAGTCTGCagcagc                                                        >  2:1853196/1‑150 (MQ=255)                                                                                                 aCGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTcc                                                     <  2:684812/150‑1 (MQ=255)                                                                                                 aCGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGGTTTTCCCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTcc                                                     >  2:1262708/1‑150 (MQ=255)                                                                                                  cGGCGAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCa                                                    <  2:189507/150‑1 (MQ=255)                                                                                                  cGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGCCGATGCAGAGTATTTCCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCa                                                    <  1:124589/150‑1 (MQ=255)                                                                                                   ggCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGCGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGCCGGCGAGGTTGCACCAGCTCCAg                                                   >  2:213332/1‑150 (MQ=255)                                                                                                           ccTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGTGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGgc                                           >  1:700913/1‑150 (MQ=255)                                                                                                           ccTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTcc                                                     <  2:765856/140‑1 (MQ=255)                                                                                                           ccTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTcc                                                     >  1:765856/1‑140 (MQ=255)                                                                                                            cTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGgcg                                          >  2:957536/1‑150 (MQ=255)                                                                                                             tCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGgcgc                                         <  2:1648770/150‑1 (MQ=255)                                                                                                              cGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGcgcc                                        >  1:1786735/1‑150 (MQ=255)                                                                                                               ggCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGcgccg                                       >  2:1491144/1‑150 (MQ=255)                                                                                                                gCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGcgccgc                                      >  1:229717/1‑150 (MQ=255)                                                                                                                 cACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCa                                     >  1:1748107/1‑150 (MQ=255)                                                                                                                    ccAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCAc                                                            >  1:928882/1‑124 (MQ=255)                                                                                                                    ccAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCAc                                                            <  2:928882/124‑1 (MQ=255)                                                                                                                    ccAGGTCTTCGAGCAGGGTGTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCt                                  >  2:1762110/1‑150 (MQ=255)                                                                                                                         tCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCgg                             <  1:184771/150‑1 (MQ=255)                                                                                                                         tCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCgg                             <  2:9348/150‑1 (MQ=255)                                                                                                                           ttCGAGCAGGGGTTTTGCCAGCAGCAGGTTGGCGGCATTGTCGTGGGCGCAGAGTAGTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGCTTGCACCCGCTCCCGCAAGCGGCGCCGCAGCTTGCgggg                           >  2:1116397/1‑149 (MQ=255)                                                                                                                              gAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCa                        <  2:1301520/150‑1 (MQ=255)                                                                                                                                     gTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTc                                                                                                                    >  2:1129328/1‑51 (MQ=255)                                                                                                                                     gTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTc                                                                                                                    <  1:1129328/51‑1 (MQ=255)                                                                                                                                        tttACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTg              <  2:229717/150‑1 (MQ=255)                                                                                                                                        tttACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTg              <  1:873990/150‑1 (MQ=255)                                                                                                                                          tACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGTGTGCAGGTCGGCTTGGa            <  1:1491144/150‑1 (MQ=255)                                                                                                                                             cagcagcagGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAg         >  1:115037/1‑150 (MQ=255)                                                                                                                                                    gggTTTGCTGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCTGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTgg  <  2:1572827/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGGCCAGTACCTCGGCACCCAGGTCTTCGAGCAGGGTTTTTACCAGCAGCAGGTTGGCGGCATTGTCGTCGACGCAGAGTATTTTCGGTCGTCGCTGGCCGTTGGGGTTGCTCCCCTCGGCCAGCGGACGGCGAGGTTGCACCAGCTCCAGCAAGAGGCGCCGCAGCTTGCGGGTGCAGGTCGGCTTGGACAGCAACTGG  >  NC_002947/1842650‑1842946

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