Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	161,068	A→C	32.0%	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	32.0%	61.1 / 28.1	67	intergenic (+47/‑149)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base C (2/20);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 1.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  CP009273/160919‑161216                                                                                                                                                      |                                                                                                                                                     ccGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCa                                                                                                                                                      >  2:827837/1‑150 (MQ=255)  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  2:4649650/1‑150 (MQ=255)  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGACAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                     >  1:4756715/1‑150 (MQ=255)     gAGGTGCAAAAAGAGATGAATGGGCGTTATCCCAAACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                  <  2:4439696/150‑1 (MQ=255)     gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGcagcag                                                                                                                                                  >  2:4438421/1‑150 (MQ=255)     gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAg                                                                                                                                                  >  1:1794714/1‑150 (MQ=255)     gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                  >  1:3659276/1‑150 (MQ=255)     gAGGTGAAAAAATATATGAACGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACGAAAAAGTATTCTTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAg                                                                                                                                                  <  2:3067908/150‑1 (MQ=38)     gAGGTGAAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCATGTGGGAGAGGGCCGGGGTGAGGGCATCAg                                                                                                                                                  <  2:3067990/150‑1 (MQ=255)      tggTGAAAAAATTGATGAATGGGGGTTATCCCAAACTTGTGGGGCCGGACGACCGGGAAAATACTGAACCGCGGCGACGCACGAAAATTTATTCTTACGTTTGGTTTTCGCCCTCTCCGTGTGGTAGTGGGCCGGGGTGAGGGCATCAgc                                                                                                                                                 <  2:4008531/149‑1 (MQ=255)            aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                           >  2:3497963/1‑150 (MQ=255)             aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                          >  2:2199151/1‑150 (MQ=255)              aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgg                                                                                                                                         >  2:3152672/1‑149 (MQ=255)              aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                         >  2:3152682/1‑150 (MQ=255)                  ggaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCCGCGCGCACGtgggg                                                                                                                                     >  1:1381300/2‑146 (MQ=255)                  agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggggc                                                                                                                                     >  1:4468946/1‑145 (MQ=255)                   gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggacc                                                                                                                                    >  1:4832571/1‑145 (MQ=255)                      gactgggCGTTACCCCAAACTGGTCGGGCCGGCCGCCCGGGCAAATACTGCCCCGCGGCTCGGTACTAAAATTTTTTCTTATGGTTGGTTTTCACCCTCCCCCGGTGGTATTGGGCGGGGGTGTGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                 <  2:329183/146‑1 (MQ=255)                       aaaGGGCGTTATCCAAAACATGTCGGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctca                                                                                                                                <  2:3868227/150‑1 (MQ=255)                           ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGGgcacccccccccc                                                                                                                            >  2:2717710/1‑149 (MQ=38)                           ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTggcccctccccca                                                                                                                            >  2:3091994/1‑149 (MQ=38)                              gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTgggcccccaccccaca                                                                                                                         >  1:398886/1‑146 (MQ=38)                                tATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGTGCACCCGCaccccacccc                                                                                                                       >  2:4030581/1‑150 (MQ=38)                                       ccccggTCGGGCCGGACGACCCGGCACTGCCTCCCCCGCCGCGTCGGACTAAAAAGTATTCGTAAGGTGGGTGTTCCCCCTCCCCCGGTGGTAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                <  2:2898601/145‑1 (MQ=255)                                            ggcgggccggccGCCCCGCAAAATCCTCCCCCGCCGCGCCGGACTAAAATTTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGTGGTTGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                           <  2:1460636/146‑1 (MQ=255)                                               tGGCCGGCCGACCGGCCAAATACGGCCCCGACCCGCCGGCCGAAAAAGTTTTCTTACGTTTGTTTTCCTCCCCCCCCCTGTGGGTGAGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                        <  2:5126997/150‑1 (MQ=255)                                               gggCCGTACGCCCGGGCAAATCCTGCACCGCCGCGCCGGCCGAAAATTTAGTCGAATGGTGGGTGTACCCCCACCCCCTGTGGGTTCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGgag                                                                                                        <  2:3736485/149‑1 (MQ=255)                                                ggCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTCCCCCCACCCCGAACCCACACCCCCCAAGgaga                                                                                                       >  2:4933279/1‑150 (MQ=25)                                                    ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCACCCCCAAACACACACCCACAACGGGGAgggg                                                                                                   >  1:1578272/1‑150 (MQ=25)                                                     gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGTGCCCCACACCCGAAACCCCGCCCCCCCGGGAggggggg                                                                                                  >  2:3903906/1‑149 (MQ=25)                                                      ccgacCCGGCAACTCCGCCCCCGCCGCCCCGGCCGAAAACGCATCCGTACGGGTGGTTTTCTCCCCCCCCCGGGGGGAGCGGGCCGGGGGGCGGCCCTCAGCGCGCCCGTTCCCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAc                                                                                                 <  2:2300492/149‑1 (MQ=255)                                                        gccccGGAAAATACTGCACCGCCGCGTGGTACTAAAAATTTTTCGTATGGTTGGTTTTCCCCCTCTCCCTGTGGGTGTGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                               <  1:4649650/148‑1 (MQ=255)                                                        gacCCGGCAAATCCGGCACCGCCGCGCCGGACGAAAAAGTATTCGAACGGTTGGTTTGCCCCCTCTCCCTGGGGGCGCGGGCCGGGGTGCGGGCCTCCGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                               <  2:4646824/150‑1 (MQ=255)                                                           ccGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCGCACCCCAACCCACAACCGCAACGGGGAGGGGGcccaga                                                                                            >  2:2914898/1‑146 (MQ=25)                                                           ccGGAACATCCTGCCCCCACGCGACGGAGGTAAATGTATTGTTTGGGTTGTTTTTCTCCCTCTCCCGGGGGGTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATc                                                                                            <  2:3149658/150‑1 (MQ=255)                                                             ggCAAATACTGCCCCGACGCGACGGCCTAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGa                                                                                          <  2:5178153/150‑1 (MQ=38)                                                               cAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCACACCCCAAACCACACCCGCAAAGGAGTGGGGTCCGaggcagc                                                                                        >  2:827984/1‑144 (MQ=25)                                                                aaaTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCAAAACCACACCCACAACGGAGAGGGGGCCGAAAGAGca                                                                                       >  1:2932315/1‑150 (MQ=25)                                                                     cTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                  <  2:1702120/150‑1 (MQ=255)                                                                                 gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCCCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                      <  2:3264106/150‑1 (MQ=255)                                                                                          aaaaGTATTGGTCAGGTTGGTTTTCTCCCGCGCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTg                                                             <  2:3903666/150‑1 (MQ=255)                                                                                               ttttCTTATGGTGGGTTTTGCCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  1:4332291/148‑1 (MQ=255)                                                                                               ttttCTTATGGGTGTTTTTCCCCCTCTCCCGGTGGGTGTGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                        <  1:1738811/148‑1 (MQ=255)                                                                                                 ttCGGAGGGTGGGTCTCCCCCCCCGCCCGGGGGGAGCGGGCCGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                      <  2:4809130/150‑1 (MQ=255)                                                                                                   cttatggttgtttTTCTCCCTCTCCGTGTGGTAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  2:2890273/145‑1 (MQ=39)                                                                                                   cGTACGGTGGGTTTTCTCCCTCGCCCTGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                    <  2:3718229/150‑1 (MQ=255)                                                                                                      acggttggttTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                 <  2:3460407/148‑1 (MQ=255)                                                                                                           tgtttttCGCCCTCTCCGGGGGGGAGAGGGCGGGGGTGTGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                            <  1:4669833/147‑1 (MQ=255)                                                                                                             gttgcctccctcgcccgGTGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGCTCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                          <  2:2058616/145‑1 (MQ=255)                                                                                                             gttTTCTCCCTCGCCGTGGGGGTGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAg                                          <  1:4219129/150‑1 (MQ=255)                                                                                                                tacgccctcgccttGGGGGTGTGGGCTGGGGGGCTGCAATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGATCACAAATTTGGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                       <  2:503619/146‑1 (MQ=255)                                                                                                                 tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACATTCACCCTCACCCTAACCCACTCCCTCAAGGGAGTGGGGACCGATCGAGCACAAATTTTGAGAGATATATTCTTATTTCTTGTTAAAGAAGAACAGaaaa                                      >  1:766282/1‑150 (MQ=38)                                                                                                                 tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACAGTCACCCACACCCTAACCCTCGACCTCCAGGGAGAGGGGCCCGATGGGGCCCAAATTTTGTTTGATATATTCGTCTGTCTGGTAACAGAAGTACAGgaac                                      >  1:769956/1‑149 (MQ=21)                                                                                                                   tccctctccctGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATc                                    <  2:4756715/150‑1 (MQ=255)                                                                                                                   tcccacgccctGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCCCGGTCACCCTCCCCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGCAGAACAGAAAATc                                    <  2:4756630/150‑1 (MQ=255)                                                                                                                         tccctGTGGGAGTGGGCCGGGGTGAGGGCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                              <  2:2886197/150‑1 (MQ=255)                                                                                                                          ccgggggggCGAGGGCGGGGGGGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAGCGAGCACAAATTTTGAGAGATATCTTCTTCGGTCTTGTAACAGAAGAACAGACAATCGGGCttt                             <  1:3720827/144‑1 (MQ=255)                                                                                                                             ggtgGGAGGGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                          <  1:3402802/149‑1 (MQ=255)                                                                                                                                 ggAGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                      <  1:4150717/150‑1 (MQ=255)                                                                                                                                      gggCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCCCACCCCAAACCAAACCCCCAACGGGGGGGGGGCCCCAGGGGCCCAAATTTTGTGAGATAAAATAATATTTTTTTGTACAGAAGAACAGAAAATAGGGCGTTTGGGcccggaat                 >  2:4313502/1‑144 (MQ=255)                                                                                                                                       ggCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                <  2:2084503/150‑1 (MQ=255)                                                                                                                                       ggCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat                <  2:1819202/150‑1 (MQ=255)                                                                                                                                        gcgggggggCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  1:3007471/147‑1 (MQ=255)                                                                                                                                        gCCGGGGTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  2:3659276/150‑1 (MQ=255)                                                                                                                                        gCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt               <  2:1794714/150‑1 (MQ=255)                                                                                                                                          gggggTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc             <  1:4839543/149‑1 (MQ=255)                                                                                                                                            gggTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg           <  2:1044133/150‑1 (MQ=255)                                                                                                                                            gggTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg           <  2:1044124/150‑1 (MQ=255)                                                                                                                                                     cATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGc  <  2:2888465/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGC  >  CP009273/160919‑161216

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