Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	161,068	A→C	40.4%	intergenic (+47/‑149)	hrpB → / → mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	161,068	0	A	C	40.4%	12.2 / 33.7	59	intergenic (+47/‑149)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base C (5/19);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 5.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGA  >  CP009273/160920‑161208                                                                                                                                                         |                                                                                                                                              cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑t                                                                                                                                             >  1:1558902/1‑150 (MQ=255)           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCAc                                                                                                                                   >  1:2795065/1‑150 (MQ=255)            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACg                                                                                                                                  >  2:2084052/1‑150 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCC‑TCCGCGCGCACggg                                                                                                                                >  1:2285063/1‑148 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACggg                                                                                                                                >  2:4081560/1‑148 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACggg                                                                                                                                >  2:4081539/1‑148 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGtg                                                                                                                                >  1:2284862/1‑149 (MQ=255)               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGACAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCCGCGCGCACGtgg                                                                                                                               >  2:210472/1‑148 (MQ=255)               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCC‑TCAGCGCGCACGtgg                                                                                                                               >  1:2857884/1‑148 (MQ=255)               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCCGCGCGCACGtgg                                                                                                                               >  2:210474/1‑148 (MQ=255)               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACgggg                                                                                                                               >  1:1647549/1‑147 (MQ=255)               agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑GCAGCGCGCTCgggg                                                                                                                               >  2:209706/1‑147 (MQ=255)                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCCAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCC‑TCAGCGCGCACGtggcc                                                                                                                             >  2:1463616/1‑146 (MQ=255)                   aTGAATGGGGGTTATCCAAAACTGGTGGGGCGGGCCGACCGGTAAAATTCTGCCCCGACGCGCGGGACTAAATTGTATTCGTATGGTTGGT‑TTTC‑TCCCTCTCCGTGGGGGTGTGGGCG‑GGGG‑TGTGGGCA‑TCAGCGCGCACGTtcaccc                                                                                                                           <  2:2129864/150‑1 (MQ=255)                     gAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTggacccac                                                                                                                         >  2:2099958/1‑148 (MQ=255)                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTGGCCCCCCACCCCaa                                                                                                                  >  2:2459732/1‑150 (MQ=38)                             gTTATCCAAAACATGGCTGGCCGGCCGACCCGGCAAATACGGCACCGCCGCGACGGACTAAAACGTATTGTTAAGTTTGGT‑TTTC‑TCCATCTCCCTGGGGGTGTGGGCG‑GGGG‑TGTGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                 <  2:2137810/150‑1 (MQ=255)                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTGGCCCCCCCCCCCAAccca                                                                                                              >  2:2371686/1‑149 (MQ=38)                                    cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTGCGCCCCACCCAAAcccaccaa                                                                                                          >  2:351530/1‑145 (MQ=38)                                               ggCCGGACGCCCGGTAAAATCCGCCCCCGACGCGCCGGACTAAAAATTTTTCTAACGTTTGGT‑TTAC‑CCCCGCCCCCTGGGGGAGTGGGCC‑GGGG‑TGAGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgaga                                                                                               <  2:553776/150‑1 (MQ=255)                                               ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTGGCCCCCCCCCCAAACCCACACCCACACgggcgc                                                                                               >  2:3149325/1‑147 (MQ=25)                                                             gCAATTCCTCCACCGACGCGACGGCCTAAAAATTAGTCTAAAGGTTGGT‑TGTC‑TCCCTCCCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAg                                                                                 <  2:4257813/150‑1 (MQ=37)                                                              cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCC‑TCCGCGCGCACGTGCCCCCACACCCCCACCCCCACCCACCAGGGCGGGGGGGCCGCTCGAGc                                                                                >  2:3592846/1‑150 (MQ=25)                                                               aaaTCCTGCACCGCCGCGACGTCCTAAAAAGTATTCGTATGGTTGTT‑TTTC‑TCCCTCACCGGGGGGGCGTGGGCG‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                               <  2:3601537/150‑1 (MQ=255)                                                                       cACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCC‑TCAGCGCGCACGTTCACCCGCACCCAAACCCACACCCACAACGGGGTGGGGAGCGATCCAGCACCACtttt                                                                       >  2:1795138/1‑150 (MQ=25)                                                                          cgccgcgacgGACTAAAATGTTTTC‑TTAGGTTGGTGTTTC‑CCCCTCGCGTGGGGGGTGGGGGCG‑GGGT‑TGGGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                    <  2:2361041/147‑1 (MQ=255)                                                                                gacggcctAAAAATTATCCGTAAGGTTGGT‑TTTC‑GCCCTCCCCGTGTGGTCGTGGGCC‑GGGG‑TGCGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                              <  2:1971784/150‑1 (MQ=255)                                                                                  ctgactAAAAATTATTCGTAAGGTTGTT‑TTTC‑ACCCTTTCCCTGGGTTAGTGGGCG‑GGGG‑TGAGGACA‑TCAGCGCGCAAGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                            <  2:1960073/148‑1 (MQ=255)                                                                                  cggacgAAAATGATTTCCTATGGTTGGT‑TTTC‑TCCCTCGCCCGGTGGGAGCGGGCG‑GGGG‑TGAGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                            <  2:604445/150‑1 (MQ=255)                                                                                     acgAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTGACCCCCACCCCAAACCACACCCACCAAGGTGATGGGCCCCAAACAGCACAAAAATTGTGAGataaaataa                                                         >  2:690415/1‑144 (MQ=16)                                                                                        aaaaaGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTCACTCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                      >  1:770710/1‑150 (MQ=255)                                                                                                     atgttgggt‑tTTC‑CCCGTCTCCCTGGGGTTGTGGGCG‑GGGG‑TGCGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                         <  1:1068065/144‑1 (MQ=255)                                                                                                        gttggt‑tTTCGCCCCTCCCCGTGGGGTTGGGGGGG‑GGGG‑TGCGGCCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                       <  1:604589/150‑1 (MQ=255)                                                                                                                ttC‑TCCCTCTCCCTGTGGTTGAGGGCC‑GGGG‑TGCGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                               <  2:2686207/150‑1 (MQ=255)                                                                                                                ttC‑TCCCGCGCCCGGTGGGAGCGGGCC‑GGGG‑TGCGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                               <  2:2685280/150‑1 (MQ=255)                                                                                                                  ccgccctccccgtGTGGGAGTGGGCC‑GGGG‑TG‑GGGCACTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaaa                              <  2:1273128/147‑1 (MQ=255)                                                                                                                   caccctctccggGGGGGGGCGGCCG‑GGGG‑GGCGGCCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGCGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                             <  1:3149325/148‑1 (MQ=255)                                                                                                                    tccctctccgtGTGGGTGTGGGCG‑GGGG‑GGAGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                            <  2:2507236/150‑1 (MQ=255)                                                                                                                     ccctctcccgGTGGTAGAGGGCC‑GGGG‑TGCGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                           <  2:2107131/150‑1 (MQ=255)                                                                                                                          gcggggtgGGGGGGGGCG‑GGGG‑TGAGGCCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                      <  1:1714688/145‑1 (MQ=255)                                                                                                                          cccgggtgGGAGCGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGCTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                      <  1:2865111/145‑1 (MQ=255)                                                                                                                           ccctGTGGGAGCGGGCG‑GGGG‑TGAGGGCA‑TCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                     <  2:3927272/150‑1 (MQ=255)                                                                                                                           ccctGGGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAGTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                     <  2:3927292/150‑1 (MQ=255)                                                                                                                             cggtgGGAGAGGGCC‑GGGG‑TGCGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                   <  2:1038151/148‑1 (MQ=255)                                                                                                                               gtgGGCGCGGCCC‑GGGG‑TGAGGCCC‑TCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCACAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                 <  1:1775855/150‑1 (MQ=255)                                                                                                                               gggggAGGGGCCC‑GGGG‑TGCGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                 <  2:456693/148‑1 (MQ=255)                                                                                                                                gggggtggggCC‑GGGGTTGCGGCCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                 <  1:591704/143‑1 (MQ=255)                                                                                                                                 gggggggggCG‑GGGG‑GGCGGCCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGcc               <  1:1704435/144‑1 (MQ=255)                                                                                                                                    tgagGGCG‑GGGG‑TTGGGCCC‑TCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCATGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa            <  2:226171/149‑1 (MQ=255)                                                                                                                                    tgagGGCG‑GGGG‑TGCGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa            <  2:1002452/149‑1 (MQ=255)                                                                                                                                     tggggCCGGGGG‑GGCGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa            <  2:809876/148‑1 (MQ=255)                                                                                                                                      tgggCC‑GGGG‑TGAGGGCA‑TCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat          <  2:2861789/149‑1 (MQ=255)                                                                                                                                        ggCG‑GGGG‑TGTGGCCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat        <  1:3857066/150‑1 (MQ=255)                                                                                                                                           c‑GGGG‑TGCGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc     <  2:2790042/150‑1 (MQ=255)                                                                                                                                             gggg‑ggCGGGCC‑TCAGCGCGCCCGTTCACCCTCACCCTAACCATCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg    <  1:921936/150‑1 (MQ=255)                                                                                                                                             gggg‑TGCGGCCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg    <  1:922173/150‑1 (MQ=255)                                                                                                                                              ggg‑TGTGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg   <  2:3022790/150‑1 (MQ=255)                                                                                                                                               gg‑TGAGGGCC‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga  <  2:1364705/150‑1 (MQ=255)                                                                                                                                               gg‑TGAGGCCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAACGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga  <  2:2486738/150‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                              CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGT‑TTTC‑TCCCTCTCCCTGTGGGAGAGGGCC‑GGGG‑TGAGGGCA‑TCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGA  >  CP009273/160920‑161208

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