Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	877,229	C→T	6.6%	R318H (CGC→CAC)	hprA ←	glycerate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	877,229	0	C	T	6.6%	173.2 / 2.7	61	R318H (CGC→CAC)	hprA	glycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base T (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

TCGGCGGGGGCGCCGGCGATCAGCAGCGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCGGGTTGCCGTTGACCGGC  >  NC_002947/877082‑877378                                                                                                                                                    |                                                                                                                                                      tcgGCGGGGGCGCCGGCGATCAGCAGCGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACgcgg                                                                                                                                                     <  1:669035/150‑1 (MQ=255)    gcggGGGCGCCGGCGATCAGCAGCGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGt                                                                                                                                                  <  2:1122863/150‑1 (MQ=255)     cggGGGCGCCGGCGATCAGCAGCGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTg                                                                                                                                                 >  2:636494/1‑150 (MQ=255)               ggCGATCAGCAGCGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGc                                                                                                                                       <  2:421484/150‑1 (MQ=255)                      agcagcGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGaagaag                                                                                                                                <  2:673392/150‑1 (MQ=255)                           cGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATCGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTg                                                                                                                           <  2:939044/150‑1 (MQ=255)                            ggCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGTCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTgg                                                                                                                          <  1:260946/150‑1 (MQ=255)                            ggCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTgg                                                                                                                          <  1:636494/150‑1 (MQ=255)                             gCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTggg                                                                                                                         <  1:823663/150‑1 (MQ=255)                               cccTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCg                                                                                                                       >  1:493686/1‑150 (MQ=255)                                 cTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGtt                                                                                                                     >  2:810105/1‑150 (MQ=255)                                      aCAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGc                                                                                                                >  2:1166305/1‑150 (MQ=255)                                      aCAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGc                                                                                                                >  2:1166491/1‑150 (MQ=255)                                          ctctGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAg                                                                                                            >  1:264217/1‑150 (MQ=255)                                               ccTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCgg                                                                                                                                        <  1:665231/117‑1 (MQ=255)                                               ccTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCgg                                                                                                                                        >  2:665231/1‑117 (MQ=255)                                                      gAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGAt                                                                                                >  1:500563/1‑150 (MQ=255)                                                      gAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGAt                                                                                                <  1:414581/150‑1 (MQ=255)                                                         cAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAgg                                                                                                                               >  2:913615/1‑116 (MQ=255)                                                         cAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAgg                                                                                                                               <  1:913615/116‑1 (MQ=255)                                                          aaCACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCACCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATgcgc                                                                                            <  1:75700/150‑1 (MQ=255)                                                            cacTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTg                                                                                          <  1:143558/150‑1 (MQ=255)                                                             acTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTgg                                                                                         >  2:237604/1‑150 (MQ=255)                                                                  aCTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCt                                                                                           <  2:994900/143‑1 (MQ=255)                                                                  aCTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCt                                                                                           >  1:994900/1‑143 (MQ=255)                                                                  aCTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGa                                                                                    >  1:349950/1‑150 (MQ=255)                                                                    tGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCgaag                                                                                                                                   <  2:166894/101‑1 (MQ=255)                                                                    tGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCgaag                                                                                                                                   >  1:166894/1‑101 (MQ=255)                                                                     gcgcGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTc                                                                                 <  2:878778/150‑1 (MQ=255)                                                                           gTCCATGGACGACGTGTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGGGGGCGCCGGGCTCAGCTGCCCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGGTGAGCTGGCCGACGATGGGCTGGGGCGAATCCACGGc                                                                           >  2:267215/1‑150 (MQ=255)                                                                           gTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGc                                                                           >  1:323632/1‑150 (MQ=255)                                                                           gTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGc                                                                           >  1:61100/1‑150 (MQ=255)                                                                           gTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGc                                                                           >  2:1206637/1‑150 (MQ=255)                                                                               aTGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCAccc                                                                       >  2:645545/1‑150 (MQ=255)                                                                                   acgacgTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAgg                                                                   >  2:251451/1‑150 (MQ=255)                                                                                   acgacgTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAgg                                                                   >  2:825436/1‑150 (MQ=255)                                                                                   acgacgTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAgg                                                                   <  2:751572/150‑1 (MQ=255)                                                                                        gTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCtg                                                              >  2:1014731/1‑150 (MQ=255)                                                                                           tcCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGc                                                           >  2:1005326/1‑150 (MQ=255)                                                                                           tcCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGc                                                           >  1:69957/1‑150 (MQ=255)                                                                                             cTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGTGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCgg                                                         >  2:1138166/1‑150 (MQ=255)                                                                                             cTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGTGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCgg                                                         >  2:1138221/1‑150 (MQ=255)                                                                                             cTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCgg                                                         <  2:61100/150‑1 (MQ=255)                                                                                              tGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCggg                                                        >  1:496959/1‑150 (MQ=255)                                                                                                gAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGt                                                      >  1:480138/1‑150 (MQ=255)                                                                                                     aaGTGGCGCAGTGTAACAGAGGCGGGCGCGTAGATCAGATGACCACCCGGCGGGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGTCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGtgatga                                                 <  2:264217/150‑1 (MQ=255)                                                                                                        tGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGTGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCa                                              <  1:1138166/150‑1 (MQ=255)                                                                                                        tGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGTGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCa                                              <  1:1138221/150‑1 (MQ=255)                                                                                                              aGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGt                                                      >  2:90426/1‑136 (MQ=255)                                                                                                              aGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGt                                                      <  1:90426/136‑1 (MQ=255)                                                                                                              aGTGCAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGGAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGCCTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGcgcg                                        <  1:155425/150‑1 (MQ=255)                                                                                                                    aCAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGc                                  <  1:1166305/150‑1 (MQ=255)                                                                                                                    aCAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGc                                  <  1:1166491/150‑1 (MQ=255)                                                                                                                      agagGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCg                                <  1:251451/150‑1 (MQ=255)                                                                                                                       gagGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGAAGATCAGGCGCGGGATGCCgg                               >  1:1102278/1‑150 (MQ=255)                                                                                                                        agGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCggg                              <  1:825436/150‑1 (MQ=255)                                                                                                                                    ggCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCgg                  >  2:877122/1‑150 (MQ=255)                                                                                                                                    ggCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGAAGATCAGGCGCGGGATGCCGGGCTCCAGCAGCgg                  <  2:1102278/150‑1 (MQ=255)                                                                                                                                     gCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCggg                 <  1:1040164/150‑1 (MQ=255)                                                                                                                                         aGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCGGGTTGc             <  2:941817/150‑1 (MQ=255)                                                                                                                                              aCCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCGGGTTGCCGTTg        >  2:1126818/1‑150 (MQ=255)                                                                                                                                                 aCGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCGGGTTGCCGTTGAcc     >  1:479458/1‑150 (MQ=255)                                                                                                                                                    cggcgTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCGGGTTGCCGTTGACcggc  <  1:645545/150‑1 (MQ=255)                                                                                                                                                    cggcgTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCGGGTTGCCGTTGACcggc  <  2:349950/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TCGGCGGGGGCGCCGGCGATCAGCAGCGGCCCCTGAAACAGCTCTGCCTGACGGAGCAACACTTCACTGCGCGGGTCCATGGACGACGTCTCCTGGAAAAAAGTGGCGCAGTGTAACAGAGGCGGGCGCCTGGCTCAGCTGACCACGCGGCGTGGCTGCCCGGCGAAGAAGGCCTGGGCGTTCTCGCTGAGCTGGCCGACGATGCGCTGGCGCGACTCCACGGCACCCCAGGCGCTGTGCGGGGTGATGATCAGGCGCGGGATGCCGGGCTCCAGCAGCGGGTTGCCGTTGACCGGC  >  NC_002947/877082‑877378

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