Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,129,437	A→T	5.7%	A40A (GCT→GCA)	gcvP‑I ←	glycine dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,129,437	0	A	T	5.7%	258.9 / 4.1	88	A40A (GCT→GCA)	gcvP‑I	glycine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (48/35);  new base T (3/2);  total (51/37)
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.95e-01

GGTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTGGCT  >  NC_002947/1129290‑1129586                                                                                                                                                    |                                                                                                                                                      ggTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCACCCGTGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGc                                                                                                                                                     >  1:738215/1‑150 (MQ=255)  gtgtTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGAGTCGGACGCGCTCTGGCCGTCTTGCGAGCCAAGCACGCTGCTACCCTTGATGCTGTCGGGGGTGACCGCAGCg                                                                                                                                                    >  1:151600/1‑150 (MQ=255)       gtagtaGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCAt                                                                                                                                               >  2:1579014/1‑150 (MQ=255)        tagtagCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGGGCGCGGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATg                                                                                                                                              <  2:172667/150‑1 (MQ=255)         agtagCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATgg                                                                                                                                             <  2:994348/150‑1 (MQ=255)          gtagCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCTGCGGTCAtccc                                                                                                                                            <  2:2354728/150‑4 (MQ=255)             gCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTc                                                                                                                                         <  2:823384/150‑1 (MQ=255)              cccTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTcc                                                                                                                                        >  1:2345058/1‑150 (MQ=255)              cccTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTcc                                                                                                                                        <  1:1414704/150‑1 (MQ=255)                   gCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGa                                                                                                                                   <  1:1367402/150‑1 (MQ=255)                         gTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGaa                                                                                                                             >  1:78349/1‑150 (MQ=255)                           aGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATg                                                                                                                                              >  1:1534659/1‑131 (MQ=255)                           aGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATg                                                                                                                                              <  2:1534659/131‑1 (MQ=255)                            gCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAgtc                                                                                                                          >  1:753717/1‑148 (MQ=255)                            gCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACCCACGTACCCGTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGcc                                                                                                                          >  2:1531848/1‑150 (MQ=255)                              ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCg                                                                                                                                                    <  2:158132/122‑1 (MQ=255)                              ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCg                                                                                                                                                    >  1:158132/1‑122 (MQ=255)                              ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCca                                                                                                                        >  2:664371/1‑150 (MQ=255)                              ttttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCca                                                                                                                        >  2:305503/1‑150 (MQ=255)                                ttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCcag                                                                                                                       <  2:2059772/149‑1 (MQ=255)                                ttGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCcag                                                                                                                       >  1:2059772/1‑149 (MQ=255)                                    aCAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCcagtgcag                                                                                                                  >  1:108794/1‑150 (MQ=255)                                              cTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGc                                                                                                        <  1:305503/150‑1 (MQ=255)                                                  ccAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAg                                                                                                              >  2:1899438/1‑140 (MQ=255)                                                  ccAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAg                                                                                                              <  1:1899438/140‑1 (MQ=255)                                                  ccAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGt                                                                                                    <  2:1448547/150‑1 (MQ=255)                                                   cAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGtt                                                                                                   >  2:680088/1‑150 (MQ=255)                                                   cAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGtt                                                                                                   >  1:1648188/1‑150 (MQ=255)                                                            ccTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCAGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAgcgg                                                                                          >  2:329274/1‑150 (MQ=255)                                                              ttGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAgcggcg                                                                                        >  2:1588417/1‑150 (MQ=255)                                                                  gcgcAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGcgcgcg                                                                                    >  1:1035748/1‑150 (MQ=255)                                                                  gcgcAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATAGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCa                                                                                              >  1:217715/1‑140 (MQ=255)                                                                  gcgcAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATAGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCa                                                                                              <  2:217715/140‑1 (MQ=255)                                                                   cgcAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCgg                                                                                   <  1:1588417/150‑1 (MQ=255)                                                                    gcAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCggg                                                                                  >  1:1193052/1‑150 (MQ=255)                                                                     cAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGc                                                                                 <  2:151600/150‑1 (MQ=255)                                                                         gagcgCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTAt                                                                             <  1:1579014/150‑1 (MQ=255)                                                                               gTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACg                                                                       >  2:73602/1‑150 (MQ=255)                                                                               gTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACg                                                                       >  2:2009831/1‑150 (MQ=255)                                                                               gTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCCGAGTTCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACg                                                                       <  1:2101883/150‑1 (MQ=255)                                                                                    ccTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGTCCTATGTGACGGGCGa                                                                  <  2:1193052/150‑1 (MQ=255)                                                                                    ccTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGa                                                                  <  1:664371/150‑1 (MQ=255)                                                                                     cTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGAt                                                                 >  1:1096214/1‑150 (MQ=255)                                                                                     cTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGAt                                                                 >  2:1217399/1‑150 (MQ=255)                                                                                      tCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATg                                                                >  1:245258/1‑150 (MQ=255)                                                                                      tCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATg                                                                >  2:244261/1‑150 (MQ=255)                                                                                         ctctGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGTCCAGTGCAGTGAGCATGGCCTGTTCGCCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAAt                                                             <  2:753717/150‑1 (MQ=255)                                                                                            tGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGATTCGAAGCCCAGTGAAGCGCGGCTGAGCGGTACGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCg                                                          >  1:1143193/1‑150 (MQ=255)                                                                                            tGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGc                                                                    <  1:1379487/140‑1 (MQ=255)                                                                                            tGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGc                                                                    >  2:1379487/1‑140 (MQ=255)                                                                                            tGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCg                                                          <  2:1035748/150‑1 (MQ=255)                                                                                              gCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGtt                                                        <  2:601439/150‑1 (MQ=255)                                                                                              gCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCGGTGGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGtt                                                        <  2:1189607/150‑1 (MQ=255)                                                                                               ccGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTg                                                       <  1:2009831/150‑1 (MQ=255)                                                                                                cgtcgtGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTgg                                                      >  2:1592290/1‑150 (MQ=255)                                                                                                cgtcgtGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGCCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTgg                                                      >  2:1023468/1‑150 (MQ=255)                                                                                                  tcgtGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCg                                                    <  2:1096214/150‑1 (MQ=255)                                                                                                        gAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGAt                                                                 >  1:2169501/1‑131 (MQ=255)                                                                                                        gAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGAt                                                                 <  2:2169501/131‑1 (MQ=255)                                                                                                          gCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAg                                            >  2:1567095/1‑150 (MQ=255)                                                                                                              aGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTg                                        >  1:1411116/1‑150 (MQ=255)                                                                                                                     tGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCa                                 >  1:1859763/1‑150 (MQ=255)                                                                                                                       gTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAgc                                                                                            <  1:2187482/89‑1 (MQ=255)                                                                                                                       gTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAgc                                                                                            >  2:2187482/1‑89 (MQ=255)                                                                                                                       gTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATg                               <  1:329274/150‑1 (MQ=255)                                                                                                                         acggttGATATTGTTGGGGATGAACGCAGCGGTCATGGCGTGCAGAGAGGCGAAGCCCAGTGCAGTAAGCATGTGCTGTTCGCCAGCGGCGCGCGGTCCTATGTAACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGc                             >  2:1128503/5‑150 (MQ=255)                                                                                                                          ccTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCg                            <  2:1859763/150‑1 (MQ=255)                                                                                                                          ccTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCg                            <  2:1648188/150‑1 (MQ=255)                                                                                                                            tttGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCggg                          >  1:2317283/1‑150 (MQ=255)                                                                                                                                aTGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGAtt                       <  2:1794635/149‑1 (MQ=255)                                                                                                                                aTGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGAtt                       >  1:1794635/1‑149 (MQ=255)                                                                                                                                aTGCTGTAGGGGATGTCCGTAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGAttt                      >  2:1133926/1‑149 (MQ=255)                                                                                                                                 tGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTcc                     >  1:740319/1‑150 (MQ=255)                                                                                                                                    tGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCGAATGCGGTAAGGAATTCCTGTTCGTCATCTGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGACATGGCGGGGTTGCTCa                  >  1:2100372/1‑150 (MQ=255)                                                                                                                                    tGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCa                  >  2:1906064/1‑150 (MQ=255)                                                                                                                                    tGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCa                  <  2:108794/150‑1 (MQ=255)                                                                                                                                       cGGGGATGACCGCAGCGGTCATGGCGTCCCTCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGc               <  1:2306663/150‑1 (MQ=255)                                                                                                                                        ggggATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGcg              >  1:784383/1‑150 (MQ=255)                                                                                                                                          ggATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGcgtc            <  2:1411116/150‑1 (MQ=255)                                                                                                                                           caTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGTATTCCTCa                  <  1:720091/142‑1 (MQ=255)                                                                                                                                           caTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGTATTCCTCa                  >  2:720091/2‑143 (MQ=255)                                                                                                                                           caTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGcgtcg           >  1:1406752/2‑150 (MQ=255)                                                                                                                                           caTGACCGCTGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGACTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGcgtcg           >  2:1565376/2‑150 (MQ=255)                                                                                                                                            aTGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGcgtcgt          >  2:889655/1‑150 (MQ=255)                                                                                                                                                ccGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATg                                     <  1:266593/119‑1 (MQ=255)                                                                                                                                                ccGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATg                                     >  2:266593/1‑119 (MQ=255)                                                                                                                                                 cGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGGGTCGTCGTTg     >  1:716438/1‑150 (MQ=255)                                                                                                                                                 cGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTg     >  2:1055047/1‑150 (MQ=255)                                                                                                                                                    aGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTGGCt  <  1:1217399/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GGTGTTGTAGTAGCCCTGGCCGATGTAGCTTTTGAACAGCTGGTTCTTGCCAGCGATGGCCTTGAGCGCAGCGAGCGCGTCGGCCTCGCTCTGGCCGTCGTGCGAGCCAAGCACGCTGGTACCTTTGATGCTGTCGGGGATGACCGCAGCGGTCATGGCGTCCAGCGAGTCGAAGCCCAGTGCAGTGAGCATGGCCTGTTCGTCAGCGGCGCGCGGGCCTATGTGACGGGCGATGAATTCGTTGGCGGTGCCGAGGTTGATGGTCATGGCGGGATTCCTCAGGCGTCGTCGTTGGCT  >  NC_002947/1129290‑1129586

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