Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	161,068	A→C	27.5%	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	27.5%	56.4 / 13.2	62	intergenic (+47/‑149)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (18/26);  new base C (3/15);  total (21/41)
Fisher's exact test for biased strand distribution p-value = 8.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160923‑161213                                                                                                                                                  |                                                                                                                                                  gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTc                                                                                                                                                 >  1:3499574/1‑148 (MQ=255)    gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCgc                                                                                                                                              >  1:1510035/1‑148 (MQ=255)           aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACg                                                                                                                                       >  1:532145/1‑148 (MQ=255)           aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACg                                                                                                                                       >  1:1870113/1‑148 (MQ=255)           aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                       >  2:893303/1‑148 (MQ=255)           aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                       >  1:3170917/1‑148 (MQ=255)           aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACg                                                                                                                                       >  1:1870330/1‑148 (MQ=255)            agagATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                      >  2:2053818/1‑148 (MQ=255)             gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                                     >  2:1825834/1‑148 (MQ=255)             gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                                     >  2:178047/1‑148 (MQ=255)             gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                     >  2:1825530/1‑147 (MQ=255)             gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                     >  2:1854890/1‑147 (MQ=255)             gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtg                                                                                                                                     >  2:3375188/1‑147 (MQ=255)              agaagaAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTCAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgc                                                                                                                                    >  2:2872693/1‑146 (MQ=255)              agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACgggg                                                                                                                                    >  1:1334694/1‑145 (MQ=255)              agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACgggg                                                                                                                                    >  1:1334457/1‑145 (MQ=255)                 tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACCGAGGAAAAAGTATGAGAGACGTAGGTTTTCACCCTCTCCCTGTGTTATAGGGGCGGGTGGTGAGCATGAGCGCGGACgagcgac                                                                                                                                 >  2:1165966/1‑142 (MQ=255)                        gCGTTATCCCAACCAGGCCGGGCCGCCCGCCCCGGAAAATACGGCCCCGACGCGACGGACGAAAATGTTTTCGTATGGTTGTTTTTCTCCCTCACCCTGTGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTtcaccctcaccc                                                                                                                          <  2:2686250/148‑1 (MQ=255)                          tttATCCAAAACATGTCTGGCCGGACCACCCGCAAAATACTGCCCCGACGCGACGGACTAAAAATTATTCGTAAGTTTGTTTTTCTCCCTCTCCGTGTGGGTGTGGGCCGGGTTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                        <  2:3544795/147‑1 (MQ=255)                              tCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGGACCCCCACCCAAAccc                                                                                                                    >  2:1402215/1‑148 (MQ=38)                                 cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCccacccccaccccag                                                                                                                 >  2:179424/1‑145 (MQ=37)                                 aaaaCTTGTGTGGCGGGACGACCCGCAAAATAGTCCCCCGCGGCGACGTACTAAAATTTTTTGTTTTGTTTGTTTTACTCCCCCTCCGTGTGGTTGTGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                 <  2:3541450/147‑1 (MQ=255)                                 aaaaCGTGTGGGGCCGGACGACCCGGAAAATACGCCCCCGACGCGACGGACGAAAAAGTTTTCTTAGGGTTGGTTTTCTCCCTCTCCGTGTGGTCGGGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCctct                                                                                                                 <  2:3541451/147‑1 (MQ=255)                                    aCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACGTGGGCCCCCCCCCAAAcccccaacc                                                                                                              >  2:1953131/1‑142 (MQ=37)                                            ggCCGCACCCCCCGCAAAATACGGCCCCGCCGCGCCGCACTAAAAATTTTTCTTAGGGTTGTTTTACCCCCACTCCCGGGGGTTGAGGGCGGGGGGGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                      <  2:829959/148‑1 (MQ=255)                                            gcccGGCCGACCGGCAAAATACTCCACCGCGGCGCCGTACGAAAATGTTTTGGTAGGGTTGGTTTGCTCCCTCGCCCTGGGGGAGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                      <  2:829986/146‑1 (MQ=255)                                                                aCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGACCCCCACCCAAACCCACTCCCACAAGGGGGAGGGGGCCCAACGAGcccc                                                                                  >  2:3503238/1‑145 (MQ=25)                                                                 cggCACCGCCGGGACGGACTAAAATGTTTTGGTAGGGTGGGTTTTCTCCCTCGCCGGGGGGGAGTGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaa                                                                                 <  2:3358316/146‑1 (MQ=255)                                                                  ggCACCGACGCGACGGACAAAAATTTATTCTTACGGTTGTTTTACACCCTCTCCGTGTGGGTGTGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACaaa                                                                                <  2:1300411/147‑1 (MQ=255)                                                                                ggaataaaaaTTTTTCGTATGTTTGGTTTTCGCCCTCGCCGTGGGGGAGTGGGCGGGGGTGGGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                  <  1:1179109/143‑1 (MQ=255)                                                                                 gactAAAAAGTTTTCGTACGGTGGGTTTTCTCCCACTCCGTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                 <  2:3210498/148‑1 (MQ=38)                                                                                 gacgAAAATGTATTCGAACGGTGGGTTTTCTCCCTCTCCCGGTGGGAGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtct                                                                 <  2:3205695/148‑1 (MQ=39)                                                                                  acgAAAAAGTTTTCTTATGGTTGGTTTTCGCCCTCGCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                <  2:1837410/148‑1 (MQ=38)                                                                                   ctaaaatGTATTCTTTTGGTGGGTTGCCCCCCTCGCCCTGTGGGAGTGGCCGGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                               <  2:3032484/146‑1 (MQ=255)                                                                                   cgAAAATTTATTGGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGTGAGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                               <  2:1956529/148‑1 (MQ=255)                                                                                         cgTATTCGTATGGTGGTTTTTCCCCCTCGCCGTGGGGGAGTGGGCGGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGt                                                         <  2:2182934/147‑1 (MQ=255)                                                                                           ttttCGTAGGGTGGGTTTACCCCCACCCCCGGGGGGAGCGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                       <  2:3105163/146‑1 (MQ=255)                                                                                           ttttCGTACGGTTGGTTTTCCCCCGCTCCCTGGGGTAGCGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                       <  2:2888167/146‑1 (MQ=255)                                                                                           tATTCGTATGTTTGGTTTCACCCCACCCCGTTGGGGTGAGGGCGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCt                                                       <  1:1484600/148‑1 (MQ=255)                                                                                                 tatgtttggttTTCGCCCTCGCCGGGTGGGTGGGGGCGGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                 <  2:1390086/143‑1 (MQ=255)                                                                                                      ttggttTACCCCGGCTCCCGGGGGTAGGGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaag                                            <  2:1361619/148‑1 (MQ=255)                                                                                                       tgtttttCCCCCTCTCCGGGTGGGAGTGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                           <  1:2919555/145‑1 (MQ=255)                                                                                                       gggtgTTCCCCCTCCCCGGGGGGGCGTGGGCCGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                           <  1:2919574/147‑1 (MQ=255)                                                                                                        ggttTTCTCCCTCTCCCTGTGGGAAAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                          <  1:3373924/148‑1 (MQ=255)                                                                                                        ggttTTCTCCCTCTCCCTGTGGGAAAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                          <  1:3374023/148‑1 (MQ=255)                                                                                                        ggttTTCTCCCTCTCCCTGTGGGAAAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                          <  1:1607791/148‑1 (MQ=255)                                                                                                          ttttCTCCCTCGCCCTGTGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                        <  2:3544431/148‑1 (MQ=255)                                                                                                          tttgccccctccccctGTGGGTGTGGGCCGGGGGGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                        <  2:1846865/142‑1 (MQ=255)                                                                                                          tttgCTCCCTCCCCGTGGGGGAGGGCGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                        <  2:1846905/148‑1 (MQ=255)                                                                                                              cccccactccgtGGGGGTGAGGGCGGGGGTGCGGCCATCAGCGCGCCCGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAGCTTCGTCTGTCTTGTAACAGAAGAACAGaaa                                    <  2:1615393/146‑1 (MQ=255)                                                                                                                       cgtgtgggtgtgggCGGGGGTGTGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                           <  1:1500186/146‑1 (MQ=255)                                                                                                                       cctGGGGGAGTGGGCCGGGTTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                           <  2:464012/148‑1 (MQ=255)                                                                                                                       ccggggggAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCGCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                           <  2:2177007/143‑1 (MQ=255)                                                                                                                          gtgGGAGCGGGCCGGGGTGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                        <  2:2637044/148‑1 (MQ=255)                                                                                                                                 ggggggggggTGCGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                 <  1:1086719/142‑1 (MQ=255)                                                                                                                                 agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                 <  2:1443089/148‑1 (MQ=255)                                                                                                                                    gcgggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAata              <  2:3078104/145‑1 (MQ=255)                                                                                                                                     ggggggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat             <  1:2080866/146‑1 (MQ=255)                                                                                                                                     cgggggTGCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAatat             <  1:893303/146‑1 (MQ=255)                                                                                                                                         ggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg         <  1:542421/148‑1 (MQ=255)                                                                                                                                         ggTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg         <  1:542467/148‑1 (MQ=255)                                                                                                                                                ccATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa  <  2:3301206/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160923‑161213

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