Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	161,068	A→C	24.2%	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	24.2%	64.0 / 8.7	65	intergenic (+47/‑149)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (20/28);  new base C (5/12);  total (25/40)
Fisher's exact test for biased strand distribution p-value = 5.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  CP009273/160922‑161215                                                                                                                                                    |                                                                                                                                                    tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                    >  2:3133564/1‑148 (MQ=255) tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCAAAACA‑TGTCTGGCCGGACGACCCGGCAATTACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGTTTTTCTCCCTCTCCCTGGGGGAGAGGGCCGGGGTGAGGGCAt                                                                                                                                                    <  2:472359/148‑1 (MQ=255)   aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCa                                                                                                                                                  >  1:503397/1‑148 (MQ=255)      tGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcg                                                                                                                                               >  2:3446110/1‑148 (MQ=255)          aaaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCa                                                                                                                                           >  1:100947/1‑148 (MQ=255)          aaaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                           >  2:2828895/1‑148 (MQ=255)          aaaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                           >  2:1396164/1‑148 (MQ=255)          aaaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCa                                                                                                                                           >  2:1395412/1‑148 (MQ=255)           aaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCCGCGCGCAc                                                                                                                                          >  1:1802009/1‑148 (MQ=255)           aaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCAc                                                                                                                                          >  1:535090/1‑148 (MQ=255)           aaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                          >  1:535140/1‑148 (MQ=255)           aaaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCAc                                                                                                                                          >  1:809031/1‑148 (MQ=255)            aaGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCCTCAGCGCGCACg                                                                                                                                         >  1:448349/1‑148 (MQ=255)             agagATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGt                                                                                                                                        >  1:2691180/1‑148 (MQ=255)                gaTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCCGCGCGCACGtgga                                                                                                                                     >  2:3450914/1‑145 (MQ=255)                 aTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAGCAGCGCGCACgggggc                                                                                                                                    >  1:3782868/1‑143 (MQ=255)                 aTGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggg                                                                                                                                    >  1:3005020/1‑144 (MQ=255)                   gAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCAACTCCGTGTGGGTGAGGGCCGGGGTGAGGGCCGCGGCGCGCCCgggggggg                                                                                                                                  >  1:1621870/1‑141 (MQ=255)                    aaaGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtggggccc                                                                                                                                 >  2:981558/1‑141 (MQ=255)                    aaaGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtgggcccc                                                                                                                                 >  1:2460614/1‑141 (MQ=255)                       gggCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTggaccccccc                                                                                                                              >  2:919884/1‑144 (MQ=255)                             atagccAAACATTGTGTGGCGGCACGACCCGCAAAAAACGCCCCCGCCGGGCGGGACTAAAATTTTTTCTTATGGTTGTTTTACGCCCTCGCCGTGGGGGAGGGGCCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                         <  2:2741609/144‑1 (MQ=255)                              aTCCCAACCA‑GGTCTGGCCGCCCGCCCCGGCAAATCCTGCACCGACGCGACGCACGAAAATGTATTCGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGGGAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAAcc                                                                                                                       <  2:2011786/148‑1 (MQ=37)                               tCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGGGCGCCCCCACCCCAAacc                                                                                                                      >  2:3576769/1‑148 (MQ=38)                                   aaaCA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGGCCCGCACCCCAAaccacac                                                                                                                  >  2:1984326/1‑144 (MQ=37)                                   aaaCA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGGCCCCGCACCCCAAcccacac                                                                                                                  >  2:209161/1‑144 (MQ=37)                                                         ccgcaaaaTCCTGCACCGCCGCGACTGACTAAAATTTATTGTTATGTTTGTTTTACACCCTCTCCGTGTGGGTGTGGGCGGGGGTGGGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                             <  2:1879540/143‑1 (MQ=255)                                                         ccGGCAAATACTCCACCGCCGCGACGTACTAAAAAGTTTTCGTAAGGTTGGTTTTCTCCCTCTCCCTTTGTGTGAGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                             <  2:1817385/148‑1 (MQ=255)                                                                acacgccACCGCCGCGCCGGACGAAAATGTATGCGTACGGTGGGTTTGCGCCCGCCCCCGGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGca                                                                                      <  2:1679751/142‑1 (MQ=255)                                                                     ccccctccgcgacgTAATAAAATTTTTTCTTATGGTTGTTTTCCCCCCTCCCGCTGGGGTTGTGGGCGGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                 <  1:2543048/141‑1 (MQ=255)                                                                         cgccgcgacgGACGAAAAATTATTCGTAAGGTTGGTTGTCTCCCTCTCCCTGTGGGATAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                             <  2:3173227/145‑1 (MQ=255)                                                                         cgacgcgacgGACTAAAAAGTATTCGTCCGGTTGGTTTTCGCCCTGTCCCTGTGGGAGCGGGCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                             <  2:2502471/148‑1 (MQ=255)                                                                           acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTGCACCCACACCCGAAACCAATCCCACAAGGGGGAGGGGAGCCCTAGTCCCCAAATTTTgag                                                                           >  2:2855812/1‑148 (MQ=21)                                                                               gacggactAAAAATTATTCTTAAGGTGGGTTTTCCCCCTCCCCCGGTGGTAGTGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                       <  2:1911364/148‑1 (MQ=38)                                                                                 cggacgAAAAGTTTTTCGTATGGTTGGTTTTCTCCCTCTCCCTGGGGGAGTGGGCCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGatat                                                                     <  2:2460614/148‑1 (MQ=255)                                                                                  ggactAAAAATTATTGGTATGGTTGGTTTTCACCCTCTCCGTGTGGTAGAGGGCGGGGGTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                    <  2:3692219/148‑1 (MQ=255)                                                                                     cgAAAAAGTATCGTTAAGTTTGGTTTACCCCCACGCCGTGGGGTCGTGGGCGGGGGGGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                 <  1:2551437/148‑1 (MQ=255)                                                                                       aaaatGTATTCGTACGGTGGGTTTTCGCCCTCCCCCTGGGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                               <  2:2410524/148‑1 (MQ=38)                                                                                       aaaaaGTATTCTTACGGTTGTTTTTCTCCCTCTCCCTGTGGTAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                               <  2:2411364/148‑1 (MQ=255)                                                                                       aaaaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttctt                                                               <  2:490564/148‑1 (MQ=255)                                                                                         aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                             <  2:3507042/148‑1 (MQ=255)                                                                                          aagaagtCGTATGGTTGTTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATAGAGCACAAATTTTGAGAGATATCTTCTTCTg                                                            <  2:135853/142‑1 (MQ=38)                                                                                               ttggTAGGGTTGGTTTCCGCCCTCGCGGGGGGGGTGTGGGCGGGGGTGGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTg                                                       <  2:2678222/145‑1 (MQ=255)                                                                                                  ttatggttggttTTCCCCCCCTCCCTGGGGGAGCGGGCGGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTaa                                                    <  2:1125788/144‑1 (MQ=255)                                                                                                   tttggttggttTTCTCCCTCTCCGGGGGGGAGGGGGCGGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                   <  2:1802009/145‑1 (MQ=255)                                                                                                   tctgtgtggttTCCCCCCTCCCCCGGGGGGGGCGGGCGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                   <  2:809031/142‑1 (MQ=255)                                                                                                      ggttggttTTCTCCCTCTCCCTGTGGGAGAGGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                <  2:3470303/148‑1 (MQ=255)                                                                                                      ggttggtgTTCCCCCTCCCCTTGTGGGTGTGGGGCGGGGTGTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                <  2:722642/148‑1 (MQ=255)                                                                                                       gttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                               <  2:1874299/148‑1 (MQ=255)                                                                                                       gtgggttTTCTCCCTCGCCCGGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagta                                               <  2:1874407/145‑3 (MQ=255)                                                                                                       gtgggttTGCGCCCTCTCCGGGTGGTTGAGGGCGGGTGTGAGGGCCTCCGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCCCAAATTTTGAGAGATTTCTTCTTCTGTCTTGTAACagaa                                               <  2:2701810/145‑1 (MQ=255)                                                                                                          ggttTTCGCCCTCGCCCTGGGGGAGAGGGCCGGGGTGCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  2:2510243/148‑1 (MQ=255)                                                                                                          gggtttCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCAACAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                            <  2:2509093/145‑1 (MQ=255)                                                                                                            ttttCCCCCCCGCCGTGGGGGCGCGGGGGGGGGTGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                          <  1:2991138/148‑1 (MQ=255)                                                                                                               tccccccctccgtGGGGGAGAGGGCCGGGGTGCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaa                                       <  2:1245671/145‑1 (MQ=255)                                                                                                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCTCCCACACCCAAACCCGATCCCGCCAGGGGGAGGGGAGCGTTCGTGGACCAATTTTGTGTGTTTTATTCTTTTGTGGTGTTAAAGAATGACAGaaaaaag                                  >  2:2541404/1‑145 (MQ=12)                                                                                                                    ctctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                  <  2:2082372/148‑1 (MQ=255)                                                                                                                      cgccgtgtgGTCGAGGGCGGGGGTGTGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCggg                                <  2:448349/143‑1 (MQ=255)                                                                                                                        ccgggggggAGAGGGCCGGGGTGTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                              <  2:2447755/142‑1 (MQ=255)                                                                                                                            gtgggcgtgggcGGGGGGGCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                          <  1:1092063/148‑1 (MQ=255)                                                                                                                                   agGGCGGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                   <  2:1621870/148‑1 (MQ=255)                                                                                                                                   agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                   <  1:612093/148‑1 (MQ=255)                                                                                                                                          gggggtgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGc            <  2:454901/142‑1 (MQ=255)                                                                                                                                            gTGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCgg          <  2:680354/148‑1 (MQ=255)                                                                                                                                                  gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa    <  1:2828895/148‑1 (MQ=255)                                                                                                                                                    aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg  <  2:406552/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    TGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACA‑TGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAG  >  CP009273/160922‑161215

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