Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,455	+G	intergenic (‑2/+1)	glpR ← / ← glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,455	1	.	G	99.0%	366.6 / ‑4.3	95	intergenic (‑2/+1)	glpR/glpR	pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (44/50);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTT  >  NC_000913/3560308‑3560598                                                                                                                                                      |                                                                                                                                                ccGATATCGATAAACAGCGTCGAGCCTTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑Gcg                                                                                                                                                 <  1:156176/149‑2 (MQ=255)  cGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑Gcgc                                                                                                                                                >  2:399642/1‑147 (MQ=255)   gATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑Gcgcagg                                                                                                                                             >  1:1063779/1‑146 (MQ=255)   gATACCGATAAACAGCGTCGAGCCTTGGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTGGCGATCGTGCCACGGCGTTTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑Gcgcag                                                                                                                                              <  2:810871/150‑5 (MQ=255)       tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATc                                                                                                                                          >  1:663288/1‑150 (MQ=255)       tCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATc                                                                                                                                          >  2:589971/1‑150 (MQ=255)        cGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATc                                                                                                                                          >  2:1157497/1‑149 (MQ=255)         gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAgg                                                                                                                                       >  1:142209/1‑151 (MQ=255)          aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGt                                                                                                                                      >  1:948427/1‑151 (MQ=255)          aTAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGt                                                                                                                                      >  1:716473/1‑151 (MQ=255)               cAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTg                                                                                                                                  <  1:704248/150‑1 (MQ=255)                  cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGGGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTcc                                                                                                                              <  1:177444/151‑1 (MQ=255)                   gTCGGGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTg                                                                                                                                  <  2:279745/146‑1 (MQ=255)                   gTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCg                                                                                                                             >  1:749271/1‑151 (MQ=255)                     cGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGcc                                                                                                                           >  1:165164/1‑151 (MQ=255)                      gggCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCg                                                                                                                             <  1:1113858/146‑1 (MQ=255)                       aGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGTCC‑CAGGCTCAGGTTTTGCTCCGcc                                                                                                                           >  2:653515/1‑149 (MQ=255)                        ggcgtgggggATTTGTGCGGCCACTTGGCGGGCGATGCGCTCTTTTTCTCCGGTCGGGGTGGCCTTGCGATCGTGCCCCGGCGTGTTAACCGATCTGGAGGGCGGCGCCGCCCCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAg                                                                                                                         <  2:419392/144‑1 (MQ=255)                          cATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCTTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                                                                                                      >  2:1066802/1‑151 (MQ=255)                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAt                                                                                                                    >  1:353550/1‑151 (MQ=255)                            ttGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAt                                                                                                                    >  1:723192/1‑151 (MQ=255)                            tgggggATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                                     <  1:588690/148‑1 (MQ=255)                             tGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                   >  1:22157/1‑151 (MQ=255)                              ggggTTTTGCCCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                   <  2:1047452/150‑1 (MQ=255)                              ggggATTTGCCCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTCCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                                                                                                  <  2:296992/151‑1 (MQ=255)                              ggggATTTCCCCCGCTTCTTTGCGGCGGTTGCGCTCTTTTTCTTCGGTCGGGGTGGCTTTGCGATCGTGCCACGGCGTGTTACCCGACCTGCAAGCCAGCGCCCCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCAtt                                                                                                                   <  2:454612/150‑1 (MQ=255)                                 ggttgggtCCGCCCCTTTGCGGGCGGTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCCGGGCGTGTTAACCAAAGGGGAAGGCAGCACCCCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGa                                                                                                                 <  2:875789/142‑1 (MQ=255)                                  aTTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                              >  1:926036/1‑151 (MQ=255)                                   tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                                                                                                              <  1:731592/150‑1 (MQ=255)                                    ttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                            >  1:994233/1‑151 (MQ=255)                                     tGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                            <  1:777123/150‑1 (MQ=255)                                      gCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc                                                                                                           >  1:938840/1‑150 (MQ=255)                                         ccGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                          <  1:660248/148‑1 (MQ=255)                                         ccGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                                         <  1:785784/149‑1 (MQ=255)                                            ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGaa                                                                                                     <  1:287803/150‑1 (MQ=255)                                                tttGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc                                                                                                >  2:602434/1‑151 (MQ=255)                                                  tGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGGGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                               <  2:199371/150‑1 (MQ=255)                                                   gCGGGGGGTGCGCTCTTTTTCTTCGGTCGGGGTGGCCTTGCGATCGTGCCGCGGCGTGTTATCTGAATTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                               <  1:1076171/149‑1 (MQ=255)                                                   gCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt                                                                                               <  2:259068/149‑1 (MQ=255)                                                     gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                           >  1:143860/1‑151 (MQ=255)                                                     gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg                                                                                            <  1:878276/150‑1 (MQ=255)                                                      ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                                          <  1:311955/151‑1 (MQ=255)                                                        cGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                        >  2:380274/1‑151 (MQ=255)                                                         gATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCt                                                                                        <  2:1028179/150‑1 (MQ=255)                                                          aTGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                                       >  2:952699/1‑150 (MQ=255)                                                             cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                   >  1:730553/1‑151 (MQ=255)                                                             cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                   >  2:852799/1‑151 (MQ=255)                                                             cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                   >  2:893584/1‑151 (MQ=255)                                                               ctctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACgg                                                                                   >  2:295777/1‑149 (MQ=255)                                                                 ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgag                                                                                 >  2:1106239/1‑149 (MQ=255)                                                                   ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaa                                                                               <  1:634307/149‑1 (MQ=255)                                                                   ttttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGAGGCGAATAATCTGCGGGCTGACGGAGAAAt                                                                             >  2:853497/1‑151 (MQ=255)                                                                     ttcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGc                                                                           >  2:1112154/1‑151 (MQ=255)                                                                          cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCAAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGctct                                                                        <  2:951165/149‑1 (MQ=255)                                                                            ggcTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGGACTGTGGGGCGGGGCCGCACCGCAACGATG‑GCGCAGGCTCAGGTTTTGCTCCGCCCGCTCATGGTGCGCGCGGCGAATAGTAGGCGAACTGACGGAGACATGTTTTAc                                                                      >  2:758852/3‑149 (MQ=255)                                                                            gTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                     >  2:16253/1‑150 (MQ=255)                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                     >  1:797769/1‑149 (MQ=255)                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                                     >  1:1052576/1‑149 (MQ=255)                                                                             tCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCa                                                                    <  1:543543/150‑1 (MQ=255)                                                                                    gccctgccGATCGGCCCGCGGCGTTTTAACCGAACTGGAAGGCATCGCCGCCCCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGATAAATGCTCTACCAGCTCTTCg                                                            <  2:171818/144‑1 (MQ=255)                                                                                          gCGATCGTGCCCCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                                      <  2:959447/151‑1 (MQ=255)                                                                                            gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTAc                                                        <  1:846642/147‑1 (MQ=255)                                                                                            gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                    >  1:323919/1‑151 (MQ=255)                                                                                            gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                                    >  1:691991/1‑151 (MQ=255)                                                                                               cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACAt                                                  <  1:10065/150‑1 (MQ=255)                                                                                                gTGCCACGGCGTGTTAACCGACCTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTTCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATaa                                                <  2:54187/151‑1 (MQ=255)                                                                                                  gCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAAc                                               <  1:523255/150‑1 (MQ=255)                                                                                                        gCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCctgctg                                        >  2:257102/1‑151 (MQ=255)                                                                                                          gtcttaacccgccTGGAAGCCACGGGCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCGAGCTGATTGAGGTCGCGGCGAATAGTCTGCGGTTTGACGGAGAAATGCTTTACCGGCTCTTCGGTACTGACATAACCctgctg                                        <  2:908928/137‑1 (MQ=255)                                                                                                              tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                  >  1:976150/1‑151 (MQ=255)                                                                                                               aaCCGAACTGGAAGGCAGCGCCGCCCCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAc                                  <  2:801516/150‑1 (MQ=255)                                                                                                               aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 >  2:996707/1‑151 (MQ=255)                                                                                                               aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 >  1:1229775/1‑151 (MQ=255)                                                                                                               aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAAcc                                 >  2:125124/1‑151 (MQ=255)                                                                                                                   gAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGtt                             >  1:865550/1‑151 (MQ=255)                                                                                                                    aaCTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                            <  2:260119/151‑1 (MQ=255)                                                                                                                    aaCTGGAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTc                            <  1:925686/151‑1 (MQ=255)                                                                                                                        ggAAGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGat                         <  1:508088/150‑1 (MQ=255)                                                                                                                          aaGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat                      <  2:487895/151‑1 (MQ=255)                                                                                                                          aaGGCAGCGCCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataa                       <  1:1117217/150‑1 (MQ=255)                                                                                                                                gcgcCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGtt                <  2:905752/151‑1 (MQ=255)                                                                                                                                gcgcCGCACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGACAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                 >  2:787461/1‑150 (MQ=255)                                                                                                                                 cgccgcgccgCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATAc                    <  2:165162/146‑1 (MQ=255)                                                                                                                                 cgccgcACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGt                 <  2:1160101/149‑1 (MQ=255)                                                                                                                                   ccgcACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt              <  1:1046784/150‑1 (MQ=255)                                                                                                                                    cgcACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGa            <  2:1073011/151‑1 (MQ=255)                                                                                                                                     gcACCGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgt              <  2:888159/148‑1 (MQ=255)                                                                                                                                      ccccGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtgtg             <  2:1237754/146‑1 (MQ=255)                                                                                                                                      cACCGGCATGATG‑GCGCAGGATTAGGTTTTGCTCCGCCAGCTCAATGAGGTCGCGGCGAATAATCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATGCCGtt                <  2:140922/145‑1 (MQ=255)                                                                                                                                       aCCGCCATGATG‑GCGCAGGACCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGTGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACAGTTGTGACGt         >  2:990736/1‑151 (MQ=255)                                                                                                                                        ccGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGtt        >  1:971110/1‑151 (MQ=255)                                                                                                                                         cGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGt         <  1:653524/149‑1 (MQ=255)                                                                                                                                         cGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg       >  2:335957/1‑151 (MQ=255)                                                                                                                                         cGCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtg       <  1:460359/151‑1 (MQ=255)                                                                                                                                          gCCATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgt      <  1:219060/151‑1 (MQ=255)                                                                                                                                           ccATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTtgtg     <  1:894304/151‑1 (MQ=255)                                                                                                                                            cATGATG‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTgtgt    <  1:531790/151‑1 (MQ=255)                                                                                                                                             atgatg‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgtt   >  1:307300/1‑151 (MQ=255)                                                                                                                                              tgatg‑GCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGtgttt  >  2:226299/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                CCGATATCGATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATG‑GC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTGACGTTGTGTTT  >  NC_000913/3560308‑3560598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑