Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	165,929	G→T	7.9%	A326A (GCC→GCA)	pntAA ←	pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	165,929	0	G	T	7.9%	248.1 / 5.7	80	A326A (GCC→GCA)	pntAA	pyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (37/36);  new base T (4/3);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCT  >  NC_002947/165780‑166074                                                                                                                                                       |                                                                                                                                                  gCTCCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAg                                                                                                                                                   <  2:1562599/150‑1 (MQ=255)     cTTAGCCGTTTTTTTTTTTGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCAt                                                                                                                                               >  2:1443813/1‑150 (MQ=255)      ttAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCACACGCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATc                                                                                                                                              >  1:1477627/1‑150 (MQ=255)         gCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAg                                                                                                                                                   <  2:96543/142‑1 (MQ=255)         gCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAg                                                                                                                                                   >  1:96543/1‑142 (MQ=255)         gCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGc                                                                                                                                           >  2:200371/1‑150 (MQ=255)         gCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGc                                                                                                                                           >  1:1350100/1‑150 (MQ=255)          ccGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCa                                                                                                                                          <  2:1491887/150‑1 (MQ=255)             tttttGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTc                                                                                                                                       <  1:200321/150‑1 (MQ=255)                   cGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTgg                                                                                                                                 <  2:769473/150‑1 (MQ=255)                    ggACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTggg                                                                                                                                <  2:1502314/150‑1 (MQ=255)                        gacCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                            >  1:957539/1‑150 (MQ=255)                        gacCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                            <  1:200371/150‑1 (MQ=255)                          ccTGACCATCACGGCCCATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                          <  2:1350100/150‑1 (MQ=255)                          ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTgg                                                                                                                                 >  2:590783/1‑143 (MQ=255)                          ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTgg                                                                                                                                 <  1:590783/143‑1 (MQ=255)                          ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                          >  1:309146/1‑150 (MQ=255)                           cTGACCATCACGGCACATCAGGCATGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTAGAACA‑GCAGCTTCATGAAGTTCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATGGGCACTCACCTGGGCCGGCAg                                                                                                                         >  2:364435/1‑150 (MQ=255)                           cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAg                                                                                                                         >  1:1241354/1‑150 (MQ=255)                            tGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAgg                                                                                                                        >  2:202103/1‑150 (MQ=255)                              aCCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACtgta                                                                                                                                     <  1:199413/135‑4 (MQ=255)                              aCCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACtgta                                                                                                                                     >  2:199413/1‑132 (MQ=255)                                  tCACGTCACACCAGCCAGGCCGCGACTATGTCGTCTTCGAGGTGGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCTTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTc                                                                                                                  <  2:661892/150‑1 (MQ=255)                                  tCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATTCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTAAAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTc                                                                                                                  <  2:1524826/150‑1 (MQ=255)                                        cacaTCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCcg                                                                                                            >  1:738815/1‑150 (MQ=255)                                          caTCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCcgac                                                                                                          >  1:216006/1‑150 (MQ=255)                                              aggcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGAGTACAGTGCGGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGTGCCGACggtg                                                                                                      >  2:303207/1‑150 (MQ=255)                                               ggcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgatga                                                                                                     >  1:1798493/1‑150 (MQ=255)                                                gcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGtt                                                                                                                      <  1:1158813/132‑1 (MQ=255)                                                gcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGtt                                                                                                                      >  2:1158813/1‑132 (MQ=255)                                                gcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgatgat                                                                                                    <  2:365385/150‑1 (MQ=255)                                                gcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgatgat                                                                                                    <  2:365655/150‑1 (MQ=255)                                                      cgcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGcc                                                                                              <  1:1883199/150‑1 (MQ=255)                                                      cgcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGcc                                                                                              <  1:1603131/150‑1 (MQ=255)                                                       gcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCTAGGCATCGGCACTCACCTGGGCCGGCAGGTTg                                                                                                                     <  1:1223306/126‑1 (MQ=255)                                                       gcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCTAGGCATCGGCACTCACCTGGGCCGGCAGGTTg                                                                                                                     >  2:1223306/1‑126 (MQ=255)                                                       gcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCg                                                                                                                                             >  2:1513470/1‑102 (MQ=255)                                                       gcgACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCg                                                                                                                                             <  1:1513470/102‑1 (MQ=255)                                                          acgaTGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCTCGCATTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTtc                                                                                          >  1:1392948/1‑150 (MQ=255)                                                              tGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCt                                                                                      >  2:519256/1‑150 (MQ=255)                                                                tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                          <  2:1216779/112‑1 (MQ=255)                                                                tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                          >  1:1216779/1‑112 (MQ=255)                                                                tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGc                                                                                    >  1:305846/1‑150 (MQ=255)                                                                tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGc                                                                                    <  2:1798493/150‑1 (MQ=255)                                                                tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGc                                                                                    <  2:957539/150‑1 (MQ=255)                                                                tcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGc                                                                                    >  1:1065996/1‑150 (MQ=255)                                                                 cgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGGGAGTTTCA‑GCCGGCTAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGca                                                                                   >  2:477247/1‑150 (MQ=255)                                                                        gAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCCTAGGACCTCACCTGGGCCGGCAGTTTAGTCGTTCTGACGATTGTCTTGCTGTTCGCTTTCACCACCTg                                                                            <  2:1477627/150‑1 (MQ=255)                                                                        gAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTg                                                                            <  1:2700/150‑1 (MQ=255)                                                                              gatgatCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGc                                                                      <  1:519256/150‑1 (MQ=255)                                                                                 gatCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACCCACCTGGGCCGGCAGGTTGGTc                                                                                                                  >  1:375088/1‑103 (MQ=255)                                                                                 gatCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACCCACCTGGGCCGGCAGGTTGGTc                                                                                                                  <  2:375088/103‑1 (MQ=255)                                                                                            cATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGc                                                        <  2:1065996/150‑1 (MQ=255)                                                                                              tCCTTGTCGAACA‑GCAGCTTCATGAAGGCCGGCGGGTTGCGTGCGTACAGTGCCGAGGCATCGGCGCTCTCCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAg                                                      <  1:562569/150‑1 (MQ=255)                                                                                                cTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGtt                                                    >  1:1076142/1‑150 (MQ=255)                                                                                                   gTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTgcc                                                 >  1:134366/1‑150 (MQ=255)                                                                                                      gAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTTCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTgccgcc                                              >  2:1160938/1‑150 (MQ=255)                                                                                                      gAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCa                                                                                                                                          >  2:1410660/1‑58 (MQ=255)                                                                                                      gAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCa                                                                                                                                          <  1:1410660/58‑1 (MQ=255)                                                                                                        aCA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCtt                                            <  2:1518630/150‑1 (MQ=255)                                                                                                         ca‑gcagCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTg                                           <  2:1038178/150‑1 (MQ=255)                                                                                                          a‑gcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACATGGCCGATGCAGCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTgg                                          <  1:364435/150‑1 (MQ=255)                                                                                                          a‑gcagcTTCATGAAGTCCAGCACGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTgg                                          <  2:309146/150‑1 (MQ=255)                                                                                                            gcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTggg                                         <  2:1241354/150‑1 (MQ=255)                                                                                                            gcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTggg                                         <  2:1076142/150‑1 (MQ=255)                                                                                                              agcTTCATGACGTACAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCt                                       >  2:1249734/1‑150 (MQ=255)                                                                                                              agcTTCATGAAGTCCAGCAGGTTGCGTGAGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCt                                       <  1:303207/150‑1 (MQ=255)                                                                                                                 ttCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCcga                                                                                                           >  1:1681159/1‑79 (MQ=255)                                                                                                                 ttCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCcga                                                                                                           <  2:1681159/79‑1 (MQ=255)                                                                                                                    aTGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTAGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGcc                                 >  1:1612927/1‑150 (MQ=255)                                                                                                                      gAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAg                               >  2:1805681/1‑150 (MQ=255)                                                                                                                       aaGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAgg                              >  2:813985/1‑150 (MQ=255)                                                                                                                         gTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCATCGGCACTCACCTGGGCCGGCAGGTTGGACGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTc                            >  1:315559/1‑150 (MQ=255)                                                                                                                           ccAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGa                          <  2:216006/150‑1 (MQ=255)                                                                                                                            cagcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGAt                         >  1:1764477/1‑150 (MQ=255)                                                                                                                                agGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGacc                     <  2:738815/150‑1 (MQ=255)                                                                                                                                           gtACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCggg          >  2:175757/1‑150 (MQ=255)                                                                                                                                            gtgaGTGCCGATGCCTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGc         >  1:91261/4‑150 (MQ=255)                                                                                                                                              cAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCGGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCtt       >  2:1491414/1‑150 (MQ=255)                                                                                                                                               aGTGCCGATGCCTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCa                                                                                                  >  1:1111399/1‑58 (MQ=38)                                                                                                                                               aGTGCCGATGCCTCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCa                                                                                                  <  2:1111399/58‑1 (MQ=38)                                                                                                                                                  gCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATc   <  1:1805681/150‑1 (MQ=255)                                                                                                                                                   ccGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGACTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCt  >  1:17677/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                  GCTCCTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACA‑GCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTCATCT  >  NC_002947/165780‑166074

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