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%	79.7 / 23.6	67	intergenic (+47/‑149)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (27/21);  new base C (2/17);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 7.88e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160920‑161213                                                                                                                                                        |                                                                                                                                                  cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CAt                                                                                                                                                  >  2:240676/1‑150 (MQ=255) cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CAt                                                                                                                                                  >  2:240691/1‑150 (MQ=255)   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCa                                                                                                                                                >  2:2160509/1‑150 (MQ=255)    gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAg                                                                                                                                               >  1:4721589/1‑150 (MQ=255)    gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAg                                                                                                                                               >  1:4721485/1‑150 (MQ=255)    gAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAg                                                                                                                                               >  1:4425266/1‑150 (MQ=255)     aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAgc                                                                                                                                              >  1:267020/1‑150 (MQ=255)     aGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAgc                                                                                                                                              >  1:267059/1‑150 (MQ=255)       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCgcgc                                                                                                                                            >  1:3902661/1‑150 (MQ=255)       gTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGGATTCGTAAGGGTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CCTGAgcgc                                                                                                                                            >  1:3902710/1‑150 (MQ=255)         gCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAgcgcgc                                                                                                                                          >  1:2714881/1‑150 (MQ=255)          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCGCGCGCa                                                                                                                                         >  1:345230/1‑150 (MQ=255)           aaaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCGCGCGCAc                                                                                                                                        >  1:6596/1‑150 (MQ=255)            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACg                                                                                                                                       >  2:1906020/1‑150 (MQ=255)            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACg                                                                                                                                       >  1:1616711/1‑150 (MQ=255)             aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCCGCGCGCACGt                                                                                                                                      >  2:751756/1‑150 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtt                                                                                                                                     >  1:5125952/1‑150 (MQ=255)                gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCTACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCACTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACggtgt                                                                                                                                   >  2:369694/1‑146 (MQ=255)                 ttaTGAATGGGGGTTATCCCAAACATGTCTGGCCGGACGACCGGGAAAAT‑ACTGCCCCGACGCGACGGACTAAAATGTTTTGGTATGGTTGGTTTTCTCCCTCTCCCTGTGGGAGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTtcac                                                                                                                                  <  2:3926215/148‑1 (MQ=25)                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtgggg                                                                                                                                  >  1:1135123/1‑146 (MQ=255)                  gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CCTCAGCGCGCACGtggggg                                                                                                                                 >  1:29913/1‑145 (MQ=255)                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACgggggcgc                                                                                                                                >  1:689679/1‑143 (MQ=255)                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtgggccc                                                                                                                                >  1:689527/1‑144 (MQ=255)                    tgactgggCTTTCTCCCAAACTTGTGTGGCGGGACTCCCCGGAAATT‑CCTGCACCGACGCGAGTGACTAAATTTTTTACTAATGGTTGGTTTACACCCACTCCCTGTGGTTTGGGGCGGGGTTGCGGG‑‑CATCAGCGCGCACGTtcaccct                                                                                                                               <  2:4367411/145‑1 (MQ=255)                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtggggggc                                                                                                                               >  1:4492161/1‑143 (MQ=255)                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGtggcccca                                                                                                                               >  2:4784897/1‑143 (MQ=255)                      ataGGGCGTTATCCCAAACTTGTCTGGCGGCACGACCCGTAAAAT‑CCGGCACCGACGCGACGGACTAAAATGTATTCGTATGGTTGTTTTTCTCCCCCACCCTGTGGGCGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTtcaccctca                                                                                                                             <  2:3862594/148‑1 (MQ=255)                        aGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTgcacccccacc                                                                                                                           >  2:3247833/1‑150 (MQ=255)                          gggtttACCCAAAACTTGTCTGGCCGCCCG‑CCCCGCAAATCAGTCCCCCGACCCGCCGCACTAAAAAGTTTTCGTATGGTTGTTTTACTCCCCCCCCGTGGGGGAGTGGGCGGGGGTGTGGC‑‑CATCAGCGCGCACGTTcaccctcaccct                                                                                                                         <  2:2043161/146‑1 (MQ=255)                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTgggcccgcacccacaccaa                                                                                                                   >  2:3718055/1‑143 (MQ=38)                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTggacccgcaccccagacca                                                                                                                   >  2:2039189/1‑143 (MQ=38)                                aTCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTGGCCCCGCACCCCAAaccc                                                                                                                   >  2:698298/1‑149 (MQ=38)                                        cgtgggtGGCGGGCCGCCCGGGCAATT‑ACCGCCCCGCCGCGACGGCCGAAAATGTATGGGACTGGTGGTTTTTCCCCGCCGCCGGGGGGCTGTGGCCGGGGGTGTGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                           <  2:3379490/144‑1 (MQ=255)                                        cATGTGGGGCCGGACGCCCCGACAAAA‑ACTCCCCCCCCGCCCCGGACCAAATTTTTTTGTTATGTTGGGTTTTCCCCCACCCCGGGTGGGTGGGGGCGGGGGTGGGGC‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCa                                                                                                           <  2:4172100/150‑1 (MQ=255)                                         ttggggggCGGGACGACCCGCAAAAT‑CCTCCCCCGCCGCCGGGGCCGAAAATGTTTTCGTTGGGTTGTTTTTCGCCCCCGCCGGGTGGTCGTGGGGGGGGGGGGGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                          <  2:2408545/144‑1 (MQ=255)                                           gtgggcccGGCCGCCCCGGAAAAT‑CCTGCCCCGCCGCGAGGTACTAAAATTTTTGCGAAAGGTGGGTTTTCCCCCGCGCCGGGTGGTTGTGGGCCGGGGTGGGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAgg                                                                                                        <  2:2830631/144‑1 (MQ=255)                                                    gacgccCCGCCAAAT‑ACTGCCCCGCCGCGCCGGCCTAAAAAGTTTTCGTACGGTTGGTTTTCCCCCTCGCCCTGGGGGAGTGGGCCGGGGTGAGGG‑‑CATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                               <  2:3307403/150‑1 (MQ=37)                                                    gacgccCCGCCAAAT‑ACGGCCCCGACGCGCCGGACTAAAATGTATTCGTACGGTTGGTTTTCGCCCCCCCCCTGGGGGAGCGGGCCGGGGTGAGGC‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGa                                                                                               <  2:4366834/150‑1 (MQ=255)                                                       gccccGGCAAAT‑CCTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAAGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                            <  1:4501169/148‑1 (MQ=255)                                                       gacCCGGCAAAT‑ACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAAGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCg                                                                                            <  1:4501227/150‑1 (MQ=255)                                                                aaT‑ACGGCCCCGCCGCGCCGGACGAAAAATGTTTCGTAGGGTTGTTTTGCCCCCTCGCCCGGGGGGTGGGGGCGGGGGTGAGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGcac                                                                                   <  2:5148172/150‑1 (MQ=255)                                                                     cggccccgccgCGCGGGACTAAAATGTTTTCGTATGGTTGGGTGTCCCCCCCCCCGTGTGGGAGGGGGCGGGGGTGCGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                               <  1:3423828/145‑1 (MQ=255)                                                                      tGCCCCCACGCGACGGACTAAAAAGTATTCGTTCGGTTGGTTGTCTCCCTCCCCCTGGGGTTGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                              <  2:3818401/150‑1 (MQ=255)                                                                              cgcgccgGACGAAAAATTTTTCTTTGGGTTGGTTTTCGCCCACGCCGGGGGGTTGGGGGCGGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagaga                                                                      <  2:2403733/150‑1 (MQ=255)                                                                               gcgaggGACTAAAATTTATTCTTTTGTTTGTTTTTCACCCTCCCCGTGTGGTTGTGGGCGGGGGGGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                     <  2:3210416/150‑1 (MQ=255)                                                                                    ggactAAAATGTTTTCGTATGGTTGTTTTTCTCCCTCTCCGTGGGGGTGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtctt                                                                <  2:3374050/150‑1 (MQ=255)                                                                                         aaattgTATTGTTATGGTGGGTTTTCTCCCTCGCCGTGGGGGAGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttcttct                                                           <  2:820846/145‑1 (MQ=255)                                                                                          aattttATTGGTATGGTTTTTTTTCTCCTTATCCTGTTGGTTGGGGGCCGGGTTTTGGC‑‑CATCAGCGATCACGTTCACACTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTATTCTg                                                          <  1:2160509/145‑1 (MQ=255)                                                                                                 ttCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTGCACCCTAACACAAAAACAAACCCACAACGGAGAGGGGTACGATCAAGCAAAAAATTTGAGATATATAATATTATTTCTTGtt                                                   >  2:1203682/1‑149 (MQ=12)                                                                                                   cGTATGTTTGTTTTACCCCCTCCCCGGGTGGGTGTGGGCCGGGGTGCGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                 <  1:3251046/150‑1 (MQ=255)                                                                                                                ttCCCCCACCCCGTGTGGCTGCGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                    <  1:2173436/150‑1 (MQ=255)                                                                                                                tactccctctccctGTGGTTGCGGGCGGGGGTGAGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                    <  2:267059/148‑1 (MQ=255)                                                                                                                taccccctccccttGTGGTTGTGGGCGGGGGTGCGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaaa                                    <  1:2173555/146‑1 (MQ=255)                                                                                                                 tcgccccctccccGGGGGGGAGGGCGGGGGGGGGGGGCCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGaa                                     <  1:3786592/147‑1 (MQ=255)                                                                                                                  ccccctcccccgGGGGGTGTGGGCGGGGGGGTGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAAt                                  <  1:3489389/148‑1 (MQ=255)                                                                                                                         tccgtGTGGTTGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtt                           <  2:2405599/150‑1 (MQ=255)                                                                                                                          ccctGTGGTTGTGGGCGGGGGTGTGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                          <  2:4656722/150‑1 (MQ=255)                                                                                                                                 ggtgtgggCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                   <  2:803118/145‑1 (MQ=255)                                                                                                                                    gtgggCGGGGGTGAGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:3621272/148‑1 (MQ=255)                                                                                                                                    gggggCGGGGGTGCGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:1616711/148‑1 (MQ=255)                                                                                                                                    gcgggCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:2001330/148‑1 (MQ=255)                                                                                                                                    gcgggCCGGGGGGTGGG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:3619545/148‑1 (MQ=255)                                                                                                                                    gagGGCCGGGGTGAGTG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:2001180/150‑1 (MQ=255)                                                                                                                                    gagGGCCGGGGTGAGCG‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                <  2:3619552/150‑1 (MQ=255)                                                                                                                                           ggggTGCGGC‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCg         <  1:46482/150‑1 (MQ=255)                                                                                                                                            gggggcggg‑‑CCTCAGCGCGCACGGTCACCCTCACCCTCACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTGCTTCTGTCTTGTAACAGAAGAACAGAAACTCGGGCTTTTGCGCCTGAATATTGCgg        <  2:3300276/144‑1 (MQ=255)                                                                                                                                                  ggg‑‑CCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaa  <  2:4492161/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                  CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAAT‑ACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGG‑‑CATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAA  >  CP009273/160920‑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: ‑