Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	161,068	A→C	31.1%	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	31.1%	56.3 / 29.4	76	intergenic (+47/‑149)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (22/30);  new base C (4/19);  total (26/50)
Fisher's exact test for biased strand distribution p-value = 6.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160920‑161217                                                                                                                                                       |                                                                                                                                                      cGTGAGGTGCAAAAAGAGATTAATGGGGGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATAAGGCACCGACGCGACGGACGAAAAATTATACGTATGGTTGGTTTT‑ATCCCTCTCCCTGTGGGTGAGGGCCGGGGT‑GAGGGCAt                                                                                                                                                      <  2:901269/150‑1 (MQ=255) cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCAt                                                                                                                                                      >  2:332881/1‑150 (MQ=255) cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCAt                                                                                                                                                      >  2:332862/1‑150 (MQ=255) cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCAt                                                                                                                                                      >  2:3101194/1‑150 (MQ=255) cGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCAt                                                                                                                                                      >  2:4227478/1‑150 (MQ=255)   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCa                                                                                                                                                    >  2:3343716/1‑150 (MQ=255)   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCa                                                                                                                                                    >  2:214988/1‑150 (MQ=255)   tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCa                                                                                                                                                    >  2:578331/1‑150 (MQ=255)          cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCa                                                                                                                                             >  1:3190486/1‑150 (MQ=255)            aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACg                                                                                                                                           >  2:3941021/1‑150 (MQ=255)              atgtgaTTAATGGGCTTTATCCCAAACTTGTGTGGCGGGACGACCCGGCAAATACTGCACCGAGGCGCCGGACTAAAATTTTTTCGTATGGTTGGTTTT‑CTCCCTCTCCCGGTGGGTGTGGGCGGGGGT‑GTGGGCATCAGCGCGCACGtt                                                                                                                                         <  2:4562740/146‑1 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCCTCAGCGCGCACggg                                                                                                                                         >  1:3206597/1‑148 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCCTCAGCGCGCACGtg                                                                                                                                         >  1:1927950/1‑149 (MQ=255)              aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGtg                                                                                                                                         >  2:1906207/1‑149 (MQ=255)                gagaTGAACGGGCGTTATCCCAAACATGTCGGGCGGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCGTCAGCGCGCACGTtca                                                                                                                                       <  2:3304606/150‑1 (MQ=38)                gagaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGtggg                                                                                                                                       >  2:5446929/1‑147 (MQ=255)                 agaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCACATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCCCGtggcc                                                                                                                                      >  2:4722889/1‑146 (MQ=255)                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTCTTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGtgggccc                                                                                                                                    >  2:528734/1‑144 (MQ=255)                   aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCCTCAGCGCGCACgggggccc                                                                                                                                    >  2:528884/1‑143 (MQ=255)                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGtggcccca                                                                                                                                   >  2:2762481/1‑143 (MQ=255)                    tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTgtacccg                                                                                                                                   >  1:1691975/1‑149 (MQ=255)                       aaGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCCTCAGCGCGCACGTtggcccacac                                                                                                                                >  1:4534913/1‑150 (MQ=255)                          ggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTggacccaccccca                                                                                                                             >  2:3467778/1‑149 (MQ=255)                           gctttACCCAAAACTGGTCGGCCCGGCCGCCCCGCCAAATACGGCACCGACGCGAGGTACTAAAAATTATTCGTAGGTTTGTTTTT‑CCCCCTCTCCCTGGGGGCGTGGGCCGGGTT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTa                                                                                                                            <  2:2695287/147‑1 (MQ=255)                            cGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTgcccccgcccccaca                                                                                                                           >  2:826507/1‑147 (MQ=38)                                      aaCATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGGCCCCCCACCCAAAaccacaacca                                                                                                                 >  2:3312028/1‑143 (MQ=37)                                      aaCATGACTGGCCGGACGACCCGGAAAATACTGCACCGACGCGACGGACTAAAAATTATTCTAATGGTTGGTTTT‑CACCCTCTCCCTGTGGGTGCGGGCCGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                 <  2:5484097/150‑1 (MQ=255)                                          tgtgtgGCCGGACGCCCCGGAAAATCCGGCACCGACGCGCGGGACTAAAATGTTTTCGTACGGTTGTTTTT‑CCCCCTCGCCGTGTGGGAGCGGGCCGGGGT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAg                                                                                                             <  2:359601/150‑1 (MQ=255)                                             ggggCCGGACGCCCGGTCAAACACTGCACCGACGCGAGGGACTAAAATTTTTTCGTATGGTTGTTTTT‑CGCCCGCCCCGGGGGGGTGCGGGCGGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  2:370467/148‑1 (MQ=255)                                             cTGGCCGGACGCCCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGTTTTT‑CTCCCTCACCCTGTGGGAGTGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  2:4593020/150‑1 (MQ=255)                                             cTGGCCGGACGACCCGGCAAATACTGCACCGCCGCGAGGGACGAAAAATTATTCGTAAGGTTGGTTTG‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGga                                                                                                          <  2:5075689/150‑1 (MQ=255)                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCCGCGCGCACGTGGGCCCCCCCCCCAACCCCCCCCCCCCAGGGAGCgggg                                                                                                    >  2:4421334/1‑150 (MQ=25)                                                   ggACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGCACCCACACCCAAACCCACACCCACAACGGCGAgggg                                                                                                    >  2:4584613/1‑150 (MQ=37)                                                    gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTGGACCCGCACCCAAAACCACAACCGCAAGGGGGAggggt                                                                                                   >  1:2677222/1‑149 (MQ=37)                                                      cgacCCGCCAAACACTGCCCCGACGCGACGCACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCGCTCCCTGGGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGAcc                                                                                                 <  2:798989/150‑1 (MQ=38)                                                        ccccgcaaaaTACGCCCCCGACGCGACGCACCAAAATTTATCCGTATGGTGGGTTTT‑CACCCTCTCCGTGGGGGTGGGGGCCGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                               <  1:4059100/143‑1 (MQ=255)                                                        acCCGGCAAATACTGCACCGACGCGACGGACGAAAATGTATTCGTACGGTTGGTTTT‑CTCCCACGCCGTGTGGGAGTGGGCCGGGGT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGa                                                                                               <  2:4085070/150‑1 (MQ=255)                                                         cccGCCAAATACGGCCCCGCCGCGAGGCACTAAAATGTTTTCGTATGGTTGGTTTC‑CGCCCGCGCCGTGGGGGAGTGGGCGGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                              <  2:2226669/150‑1 (MQ=255)                                                              cacctACGGCCCCGCCGCCCCGGACGAAAAAGTTTGCGTACGGTTGGTTTG‑CCCCCTCGCCCGGTGGGGGCGGGCCGGGGT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                         <  2:497645/146‑1 (MQ=255)                                                              aaaaTACGGCCCCGACGCGCCGGACCAAAATGTATTCGAAGGGTTGGTTTT‑CCCCCTCACCCTGTGGTAGCGGGCCGGGGT‑GAGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGc                                                                                         <  2:630936/149‑1 (MQ=255)                                                                    ctccACCGCCGCGCCGGACCAAAATTTATTGTTAGGTTTTTTTTT‑CCCCCCCCCCGTGTGGGAGGGGGCGGGGGT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                   <  1:2679802/147‑1 (MQ=255)                                                                    cTGCACCGCCGCGACGGACGAAAAAGTATTCGTACGGTTGGTTTT‑CTCCCTCTCCGTGGGGTTGCGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                   <  2:3502983/150‑1 (MQ=255)                                                                     cgCACCTACGCGCCGTCAGAAAATGTTTTC‑TTATGTTGGTTTTCCGCCCTC‑CCCGGTGGGTATGGGCCCGGGTGGTGGCCATCAGCGAGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtt                                                                                  <  2:4488152/149‑1 (MQ=255)                                                                      gccccGACGCGACGGACTAAAACGTTTTCGTCTGGTTGGTTTA‑CCCCCACACCCTGGGGGAGCGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttc                                                                                 <  1:5055744/147‑2 (MQ=255)                                                                      gCACCGACGCGACGGACTAAAATTTATTCGTAGGGTTGTTTTA‑CCCCCACGCCCTGTGGGATAGGGCCGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                 <  1:5055723/150‑1 (MQ=255)                                                                      gCACCGACGCCCGGCAAGAAAAGGTATTCTTATGTTTTGTTTT‑CCCCCACGCCGGGTGGTTGTGGGCGGGGGT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAttt                                                                                 <  2:3144277/150‑1 (MQ=255)                                                                           gccgcgccgCACTAAAATTTTTTCTTATGTTTGTTTTT‑CGCCCCCTCCGGGTGGGTGTGGGCCGGGGT‑GTGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgaga                                                                            <  1:3598330/148‑1 (MQ=255)                                                                            acgcgacgGACTAAAATTTATTCGTAAGGTTGGTTTT‑CGCCCTCTCCGTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                           <  2:4048407/150‑1 (MQ=38)                                                                                 acggacgAAAAAGTTTTCGTAAGGTGGGTTTT‑CACCCTCTCCCTGGGGGTGTGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtc                                                                      <  2:1818506/150‑1 (MQ=39)                                                                                    gactAAAAAGTATGGTTAGGGTTGTTTTT‑CTCCCGCTCCGTGGGGGTGGGGGCGGGGGT‑GAGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATAtcttc                                                                   <  1:4963261/150‑1 (MQ=255)                                                                                        aaaaaGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCACACCCAAAACCACTCCCACAACGGCGCGGGGTCCGAACCAGCACCAAATTTTAGAGGTAtattcttct                                                               >  2:740927/1‑150 (MQ=21)                                                                                         taacgTATTCTTACGGTTGTTTTT‑CTCCCTCACCGGGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTg                                                              <  2:2247071/146‑1 (MQ=38)                                                                                             gTATTCGTATGGTTGGTTTT‑CTCCCTCTCCGTGTGGGAGTGGGCGGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                          <  2:5430706/150‑1 (MQ=255)                                                                                                ttggTATGGTTGGTTTT‑CTCCCTCGCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                       <  2:5056077/147‑1 (MQ=255)                                                                                                ttggTATGGTGGGTTTG‑CGCCCTCCCCCGGGGGGTGTGGGCCGGGGG‑GGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                       <  2:3228227/147‑1 (MQ=255)                                                                                                ttCGTATGGTGGGTTTT‑CCCCCACCCCCGGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTa                                                       <  2:5056040/150‑1 (MQ=39)                                                                                                   gTATGGTGGTTTTT‑CTCCCTCGCCGGGTGGGTGTGGGCGGGGGT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACa                                                    <  2:2886830/150‑1 (MQ=255)                                                                                                         ttgtttTG‑CCCCCTCGCGGTGGGGGTGTGGGCGGGGGT‑GCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAgaagaa                                              <  2:5159177/150‑1 (MQ=255)                                                                                                          tgttttt‑CTCCCTCTCCGGGTGGGTGTGGGCGGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  1:429135/147‑1 (MQ=255)                                                                                                          ggttttt‑CACCCCCCCCGTGTGGGTTTGGGCGGGGGT‑GAGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  1:3194799/147‑1 (MQ=255)                                                                                                          gggttTT‑CCCCCCCCCCCTGGGGGAGTGGCCGGGGTT‑CCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                             <  2:1456595/149‑1 (MQ=255)                                                                                                                    ccctcgccgtGGGGGTGGGGGGGGGGGT‑GCGGGCCTCAGCGCGCACGGTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCg                                    <  2:2008920/150‑1 (MQ=255)                                                                                                                       tctcccgGGGGTAGTGGGCGGGGGT‑GTGGGCCTCAGCGCGCCCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGc                                 <  2:2379744/150‑1 (MQ=255)                                                                                                                          ccctGTGGGCGTGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCttt                              <  2:3672325/150‑1 (MQ=255)                                                                                                                               gtggtttgggCCGGGGTGGCGGCCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgcg                          <  1:334799/143‑1 (MQ=255)                                                                                                                                 gtggggggCGGGGGG‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:3312028/145‑1 (MQ=255)                                                                                                                                 ggtgtgggCGGGGGT‑GTGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCt                       <  1:4143990/145‑1 (MQ=255)                                                                                                                                   tgggggCGGGGGT‑GCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGa                     <  1:451563/147‑1 (MQ=255)                                                                                                                                    gtgggCCGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGaa                    <  1:5243845/148‑1 (MQ=255)                                                                                                                                     cgggCGGGGGT‑GGGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                   <  2:2084920/149‑1 (MQ=255)                                                                                                                                     cgggCCGGGTT‑CAGGCCCTCAGCGCGCACGTTCACCCTCACCCTACCCCTCTCCCTCAAGGGAGAGGGGACCGATCGATCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAat                   <  2:2085135/149‑1 (MQ=255)                                                                                                                                        gCCGGGGT‑GCGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATAtt                <  2:2629342/150‑1 (MQ=255)                                                                                                                                         cgggggT‑GCGGGCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTg               <  2:1241688/148‑1 (MQ=255)                                                                                                                                             ggT‑GCGGCCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGga           <  1:289022/150‑1 (MQ=255)                                                                                                                                                 gtgggCCTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAATATCGGGCTTTTGCGCCTGAATATTGCGGAGaa        <  1:1285949/148‑1 (MQ=255)                                                                                                                                                 gtgggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGATCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa        <  1:1285722/148‑1 (MQ=255)                                                                                                                                                       aTCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCa  <  2:3323797/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      CGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTT‑CTCCCTCTCCCTGTGGGAGAGGGCCGGGGT‑GAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCA  >  CP009273/160920‑161217

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