Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,141,257	T→G	15.8%	A590A (GCA→GCC)	ppc ←	phosphoenolpyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,141,257	0	T	G	15.8%	192.5 / 2.4	99	A590A (GCA→GCC)	ppc	phosphoenolpyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base T (47/29);  new base G (5/14);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 8.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAA  >  CP009273/4141110‑4141404                                                                                                                                                      |                                                                                                                                                    cAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTgc                                                                                                                                                   <  1:932064/150‑1 (MQ=255) cAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTg                                                                                                                                                    <  1:79050/149‑1 (MQ=255)  aGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTgcg                                                                                                                                                  >  1:1018497/1‑150 (MQ=255)   gAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTgcgc                                                                                                                                                 <  1:180005/150‑1 (MQ=255)    aaTCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGcgcc                                                                                                                                                <  1:526600/150‑1 (MQ=255)     aTCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGcgccg                                                                                                                                               >  1:905048/1‑150 (MQ=255)        gCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGGGcgccgccg                                                                                                                                            >  1:514608/1‑150 (MQ=255)        gCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGGGcgccgccg                                                                                                                                            >  2:404199/1‑150 (MQ=255)          cccGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGcgccgccgc                                                                                                                                           >  1:1079243/1‑149 (MQ=255)                 taAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAAt                                                                                                                                    <  1:168630/149‑1 (MQ=255)                  aaaGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATg                                                                                                                                   >  2:725018/1‑149 (MQ=255)                   aaGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATgg                                                                                                                                  <  1:324320/149‑1 (MQ=255)                          aGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAAccgccg                                                                                                                          >  1:865504/1‑150 (MQ=255)                                  gACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTg                                                                                                                  >  1:918317/1‑150 (MQ=255)                                  gACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTg                                                                                                                  >  1:671711/1‑150 (MQ=255)                                   agggTGATTTCGGGCAGACCATATTAAAGGCGGTTCATCCCGCCTTGTTGGGTAG‑CGCGCAGGCCGCTTTTCAGCCTGCTGGGGGG‑TGGTGACAGCGGCGCCGCAGGGGGGGGGGCGCGGCGGCGACCAATGGAACCGCCGCGACCGTgg                                                                                                                 <  2:277585/148‑1 (MQ=255)                                   aCGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTg                                                                                                                  >  2:1036266/1‑149 (MQ=255)                                     ggTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGaa                                                                                                               >  2:156241/1‑150 (MQ=255)                                      ttgttttCGGGCAGCCAATATTAAAAGGGGATCACCCCCCCTTTTTGGGTTC‑CGCGC‑GGCCCCCTTTTCGGGTTTCTGGGGGTTTGGGACAGCGGCGCCGCATGGGCGGGTCCGCCGCCGCACCAAATGGAACCGCCGCGACCGTGGaac                                                                                                              <  2:514608/146‑1 (MQ=255)                                        gATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGaacaa                                                                                                            >  2:153677/1‑150 (MQ=255)                                        gATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑CTGTGACAGCAGCGCCGCATGAGCAGGCGCGCCGCCGCGACCAATGGAACCGCCGCGACCGGGGaacaa                                                                                                            >  1:707631/1‑150 (MQ=255)                                         aTTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCCCGACCGTGGaacaac                                                                                                           >  2:777722/1‑150 (MQ=255)                                          tttCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCACCGCCGCATGAGCAGGTGCGCCCCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACg                                                                                                          >  2:480917/1‑150 (MQ=255)                                          tttCTGGCAGACAATATTTAAGGCGGATCACCTCGCCCTGTTCGGTTC‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGGGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACg                                                                                                          <  1:188279/150‑1 (MQ=255)                                                    aCCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGACCAATGGAACCGCCGCGAACGTGGAACAACGTCAGCTCaa                                                                                                 >  1:456850/1‑149 (MQ=255)                                                       atatTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATAcc                                                                                             >  2:1078232/1‑150 (MQ=255)                                                       atatTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATAcc                                                                                             >  2:359256/1‑150 (MQ=255)                                                        tatTAAAGGGGGATCATCCCCCCTTGTTGGGTTA‑CGCGCGGGCCGCTTTCCGGGTTCCCGGGCGG‑TTGTGACAGCAGCGCCGCATGAGCGGGGGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATAc                                                                                              <  1:201290/148‑1 (MQ=255)                                                         atTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCCCAATACCCg                                                                                           >  2:667524/1‑150 (MQ=255)                                                          ttaaaaGCGTATCCTCCCCCCCTGTTGGGTTA‑CGCGCAGGCCGCCTTTCAGGTTTCCTGGGGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCCCCGCCGCGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCg                                                                                           <  2:1027407/146‑1 (MQ=255)                                                             aaggCGGTTCACCCCGCCCT‑TTTGGTTACCGCGCAGGCGCCCTTTCAGGCTCCCGGGGGG‑TTGTGACAGCAGCGCCCGAGGAGCGGGGGCGCCGCCGCACCCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCtt                                                                                        <  2:1082218/146‑1 (MQ=255)                                                              aaGGGGATCATCTCCCCTTTTTGGGTTC‑CGCGCGGGCCGCTTTTCAGGTTTCCGGGGGT‑TTGGGACAGCAGCGCCGCAGGGGCGGGTGCGCCCCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCtttt                                                                                      <  1:455990/150‑1 (MQ=255)                                                              aaGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCGGGCGGCCTTTCGGGCTTCCTGGCGG‑TTGTGACAGCGGCGCCGCATGAGCGGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCttt                                                                                       <  2:649687/149‑1 (MQ=255)                                                                gCGGTTCATCTCCCCTTGTTCGGTCG‑GGCGCAGCCCCCCTTTCGGGTTCCCGGGGGG‑TTGTGACACCGGCGCCGCATGACGAGGGGCGCCGCCGCGACCAGTGGACCCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTc                                                                                    <  2:186781/150‑1 (MQ=255)                                                                gCGGATCATCCCGCCGTGTTCGGTTA‑CGCGCAGGCCCCTTTCCAGGCTCCGGGGGGG‑TTGTGACACCGGCGCCGAATGAGCGGGTCCGCGGCCGCCACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTc                                                                                    <  1:753811/150‑1 (MQ=255)                                                                 cGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCg                                                                                   >  2:106159/1‑150 (MQ=255)                                                                    atcatcTCGCCTTGTTGGGTTC‑CGCGCGGCCCGCTTTCCAGGCTTCCGGGGGG‑TTGTGACAGCAGCGCCGAAGGGGCGGGTGCGCCGCCGCGCCAAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAg                                                                                <  2:403564/150‑1 (MQ=255)                                                                    atcaccCCGCCCGGTGCGGATA‑CGCGGGGGCCGCCTTTAAGGCTTCCGGGGGG‑TGGAGACAGGAGCGCGCAAAGAAAAGGTGCGCCGCCGCGACAAAGGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAg                                                                                <  2:61049/150‑1 (MQ=255)                                                                     tcatcTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAgg                                                                               >  1:451209/1‑150 (MQ=255)                                                                     tcatcTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAgg                                                                               >  1:868961/1‑150 (MQ=255)                                                                       atcTCGCCCTGTTCGGTTA‑CGCTCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGtt                                                                             >  2:552819/1‑150 (MQ=255)                                                                       acctcGCCCTTTTCGGTTC‑CGCGCGGCCCGCCTTTCGGGCTTCCTGGGGG‑TTGTAACAGCGGCGCCGCATGGGCGGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGtt                                                                             <  2:865504/148‑1 (MQ=255)                                                                        tctcGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGttt                                                                            >  1:957726/1‑150 (MQ=255)                                                                        tctcGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGttt                                                                            >  1:611143/1‑150 (MQ=255)                                                                         ctcGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGtttt                                                                           >  1:733595/1‑150 (MQ=255)                                                                         ctcGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGtttt                                                                           >  2:38348/1‑150 (MQ=255)                                                                            gCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAt                                                                        >  1:122640/1‑150 (MQ=255)                                                                            gCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAt                                                                        >  1:924787/1‑150 (MQ=255)                                                                             cccTGTTCGGTTC‑CGCGCGGGCCCCCTTGCGGCCTTCCGGGGGG‑TTGGGACGGCAGCGCCGCATGAGCGGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAtt                                                                       <  2:1079243/150‑1 (MQ=255)                                                                             cccTGTTCGGTTA‑CGCGCAGGCCGCTTTTCAGGCTTCCGGGGGG‑TTGTGACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAtt                                                                       <  1:438228/150‑1 (MQ=255)                                                                             cccTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAt                                                                        >  1:579825/1‑149 (MQ=255)                                                                              ccTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCCTGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGCGGAACAACGTCAGCTCAATACCCGCCTTTTCGGCGGTTTTGAtt                                                                       >  2:1006834/1‑149 (MQ=255)                                                                              ccTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGAtt                                                                       >  1:695336/1‑149 (MQ=255)                                                                                 gtggggTTC‑CGCGCAGGCCGCCTTCCGGGTTTCCGGGGGG‑TTGGGCCAGCGGCGCCGCGTGGGAGGGGGCGCCGCCGCGACAATTGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAAt                                                                    <  1:359256/145‑1 (MQ=255)                                                                                 gTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATg                                                                   >  1:188570/1‑150 (MQ=255)                                                                                  ttCGGTTA‑CGCGCGGGCCGCCTTTCGGGCTTCCGGGCGG‑TTGGGACAGCAGCGCCGAAGGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGc                                                                  <  2:46255/150‑1 (MQ=255)                                                                                    tggtttccgcgcgGGCCGCCTTTCAGGCTTCGGGGGGG‑TTGTAACAGCAGCCCGGCAAGGGCGGGCGCCCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGc                                                                  <  1:726306/142‑1 (MQ=255)                                                                                    gggTTA‑CGCGCAGCCGCCCTTTCAGCCTTCGTGGTGG‑TTGTGCCGGCGCCGCCGTGAAAGCAGGTGGCCCGCCGAGACCAAAGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCa                                                                 <  2:944404/148‑1 (MQ=255)                                                                                    cGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAACAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGCGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCAt                                                                >  1:106333/1‑150 (MQ=255)                                                                                      ggtcccgcgcgGGCCGCTTTTCGGGCTGCCGGGCGG‑TTGGGAAAGCAGCGCCGCATGAGCGGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATc                                                               <  2:10212/145‑1 (MQ=255)                                                                                      gTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATcc                                                              >  1:826920/1‑150 (MQ=255)                                                                                       ttA‑CGCGCAGGCCGCTTTTCAGGCTTCCGGGGGT‑TTGGGACAGCAGCGCCGCATGAGCGGGGGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCt                                                             <  1:519723/150‑1 (MQ=255)                                                                                          ccgcgcgGCCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCGGCGCCGCAGGGGCAGGTGCGCCGCCGCGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCtg                                                            <  2:447786/148‑1 (MQ=255)                                                                                           cgcgCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCt                                                             >  1:447786/1‑147 (MQ=255)                                                                                            gcgcGGGCCCCTTTTCAGGCTTCCTGGGGG‑TTGTGCCAGCGGCGCCGCATGGGAGGGGGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGc                                                         <  2:868961/150‑1 (MQ=255)                                                                                               cgggCCGCCTTTCGGGTTTCCGGGGGT‑TTGTGACAGCAGCGCCGCATGAGCGGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTg                                                      <  1:354114/148‑1 (MQ=255)                                                                                                 ggCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGa                                                     >  1:1024033/1‑149 (MQ=255)                                                                                                 gcccgctTTTCGGGCTCCGGGGGGG‑TTGTGACAGCAGCGCCGCATGGGCGGGGGCCCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGat                                                    <  2:402103/148‑1 (MQ=255)                                                                                                        tttCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTgcgc                                             >  1:742070/1‑150 (MQ=255)                                                                                                          tCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGcc                                            >  1:1065242/1‑149 (MQ=255)                                                                                                          tCAGGCTCCCGGGGGG‑TTGGGACAGAGGCGCCGCAGGGGCAGGGGCGCCGCCGCAACAAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGcc                                            <  1:667524/149‑1 (MQ=255)                                                                                                            aGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAg                                         >  1:201507/1‑150 (MQ=255)                                                                                                             ggCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAg                                         >  1:957057/1‑149 (MQ=255)                                                                                                             ggCTTCCTGGCGG‑TTGGGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAg                                         >  1:821934/1‑149 (MQ=255)                                                                                                             ggCTTCCGGGGGG‑TTGTGACAGCAGCGCAGCAGGAGCGGGGGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAgg                                        <  2:918317/150‑1 (MQ=255)                                                                                                              gCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGa                                       <  1:725018/150‑1 (MQ=255)                                                                                                               cTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTTCGCCCAgg                                        >  2:520473/1‑148 (MQ=255)                                                                                                                ttCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGGGCGCCGCCGCGCCCAAAGGGACCGCCGCGACCCTGGGACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGGTTAATGCATCCTGTGCCTGATATTTCGCCCAGGAAg                                     >  2:298145/1‑150 (MQ=255)                                                                                                                  ccTGGGGT‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCCCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                   <  2:905048/150‑1 (MQ=255)                                                                                                                  ccTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCt                                   >  1:104790/1‑150 (MQ=255)                                                                                                                       ggg‑TTGTGACAGCGGCGCGGCAGGGGCGGGTGCCCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCAt                              <  2:122640/149‑1 (MQ=255)                                                                                                                        cgtttGTGACAGCAGCGCCGCATGAGCGGGTGCCCCGCCGCGCCAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCAt                              <  1:157936/147‑1 (MQ=255)                                                                                                                           ttGGGACGGCCGCCCCGCGGGAGCGGGCGCGCCGCCGCGCCCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCAc                           <  2:604681/150‑1 (MQ=255)                                                                                                                             gtgACAGCACCGCCGCAGGAGCAGGTGCGCCCCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACt                          <  2:957726/149‑1 (MQ=255)                                                                                                                             ggggcaccagcGCCGGATGGGGGGGGGCGCCGCCGCACCAAAGGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTc                         <  1:480917/143‑1 (MQ=255)                                                                                                                                 cagcagCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGc                     <  1:466290/150‑1 (MQ=255)                                                                                                                                  cgcagcGCCGCATGAGCAGGGGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCa                    <  2:675381/149‑1 (MQ=255)                                                                                                                                   gcggcgcCGCAGGGGCGGGTGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCAt                   <  2:1051857/147‑1 (MQ=255)                                                                                                                                     agcGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCt                 >  1:557408/1‑150 (MQ=255)                                                                                                                                     agcGCCGCAGGAGCGGGTGCCCCGCCGCGACAAATGAAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCt                 <  1:777722/150‑1 (MQ=255)                                                                                                                                        gccgcATGAGCAGGGGCGCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCtttt              <  1:1036266/150‑1 (MQ=255)                                                                                                                                        gccgcAGGAGCGGGCGCCCCGCCGCGACAAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCtttt              <  2:611143/150‑1 (MQ=255)                                                                                                                                             aTGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCt           >  1:504633/1‑148 (MQ=255)                                                                                                                                             aTGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCt           >  2:480534/1‑148 (MQ=255)                                                                                                                                             aTGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCt           >  2:950044/1‑148 (MQ=255)                                                                                                                                                aGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTc      >  1:168167/1‑150 (MQ=255)                                                                                                                                                aGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTc      >  2:379083/1‑150 (MQ=255)                                                                                                                                                   aGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCgg    >  1:131506/1‑149 (MQ=255)                                                                                                                                                    ggTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGaa  >  1:692276/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CAGAATCGCCCCGGTATAAAGCGACAGGCTGCTGACGGTGATTTCTGGCAGACCATATTTAAAGCGGATCATCTCGCCCTGTTCGGTTA‑CGCGCAGGCCGCCTTTCAGGCTTCCTGGCGG‑TTGTGACAGCAGCGCCGCATGAGCAGGTGCGCCGCCGCGACCAATGGAACCGCCGCGACCGTGGAACAACGTCAGCTCAATACCCGCTTTTTCGCAGGTTTTGATTAATGCATCCTGTGCCTGATATTGCGCCCAGGAAGCTGCCATCACTCCCGCATCTTTTGCTGAGTCGGAA  >  CP009273/4141110‑4141404

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