Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,980,852 G→T 11.8% R131S (CGT→AGT)  hemY ← putative protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,980,8520GT11.8% 1118.7 / 76.3 333R131S (CGT→AGT) hemYputative protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base G (124/169);  new base T (16/24);  total (140/193)
Fisher's exact test for biased strand distribution p-value = 8.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 358 total aligned reads displayed.

AGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCACGGCGA  >  CP009273/3980704‑3981000
                                                                                                                                                    |                                                                                                                                                    
 gcagcTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACgt                                                                                                                                                     <  1:3684259/149‑1 (MQ=255)
    gcTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACgttgt                                                                                                                                                  <  2:4626370/149‑1 (MQ=255)
    gcTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACgttgt                                                                                                                                                  <  2:3102060/149‑1 (MQ=255)
     cTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACgttgtt                                                                                                                                                 <  1:2739900/149‑1 (MQ=255)
       tATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGGTGTTgc                                                                                                                                               >  2:780956/1‑149 (MQ=255)
         tCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTgcgc                                                                                                                                             >  1:3777074/1‑149 (MQ=255)
            aCGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGc                                                                                                                                          <  2:2421058/149‑1 (MQ=255)
             cGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCt                                                                                                                                         <  2:3910823/149‑1 (MQ=255)
                  tGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGc                                                                                                                                    <  2:1924323/149‑1 (MQ=255)
                  tGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGc                                                                                                                                    <  2:3892378/149‑1 (MQ=255)
                   ggCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCa                                                                                                                                   <  2:2415260/149‑1 (MQ=255)
                    gCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAg                                                                                                                                  >  2:3822358/1‑149 (MQ=255)
                     cGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGc                                                                                                                                 <  1:613542/149‑1 (MQ=255)
                      gTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGcc                                                                                                                                >  2:2687336/1‑149 (MQ=255)
                        gcagcaTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAg                                                                                                                              >  2:286896/1‑149 (MQ=255)
                        gcagcaTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAg                                                                                                                              <  2:261904/149‑1 (MQ=255)
                         cagcaTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGt                                                                                                                             <  2:3816686/149‑1 (MQ=255)
                          agcaTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTa                                                                                                                            <  2:2793932/149‑1 (MQ=255)
                            caTGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGa                                                                                                                          >  1:26485/1‑149 (MQ=255)
                              tGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATa                                                                                                                        >  1:3852168/1‑149 (MQ=255)
                               ggTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAg                                                                                                                       >  2:4540034/1‑149 (MQ=255)
                                gTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGt                                                                                                                      >  1:35051/1‑149 (MQ=255)
                                 ttttCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGtt                                                                                                                     <  2:400226/149‑1 (MQ=255)
                                    tCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTcac                                                                                                                  <  2:3533365/149‑1 (MQ=255)
                                       ttACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTcaccac                                                                                                               >  1:4164052/1‑149 (MQ=255)
                                        tACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCaccacc                                                                                                              >  2:3400/1‑149 (MQ=255)
                                        tACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCaccacc                                                                                                              <  1:3182351/149‑1 (MQ=255)
                                        tACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCaccacc                                                                                                              <  1:4252567/149‑1 (MQ=255)
                                           gggCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgt                                                                                                           >  1:3743836/1‑149 (MQ=255)
                                            ggCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgtt                                                                                                          <  1:2973653/149‑1 (MQ=255)
                                            ggCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgtt                                                                                                          >  1:3871362/1‑149 (MQ=255)
                                             gCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgttg                                                                                                         >  2:209911/1‑149 (MQ=255)
                                              ccAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgttgt                                                                                                        <  2:1390533/149‑1 (MQ=255)
                                              ccAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGgttgt                                                                                                        <  1:2539462/149‑1 (MQ=255)
                                                 gTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTcc                                                                                                     >  2:3753309/1‑149 (MQ=255)
                                                  ttGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcg                                                                                                    >  2:4405743/1‑149 (MQ=255)
                                                  ttGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcg                                                                                                    >  2:3995575/1‑149 (MQ=255)
                                                  ttGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcg                                                                                                    >  2:2014071/1‑149 (MQ=255)
                                                   tGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcgc                                                                                                   <  1:1799410/149‑1 (MQ=255)
                                                    gCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCcgcg                                                                                                  <  2:4195516/149‑1 (MQ=255)
                                                     cAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGt                                                                                                 >  1:2877289/1‑149 (MQ=255)
                                                       gACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGa                                                                                               >  1:4492473/1‑149 (MQ=255)
                                                       gACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGa                                                                                               >  2:2820221/1‑149 (MQ=255)
                                                       gACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGa                                                                                               >  2:1302369/1‑149 (MQ=255)
                                                        aCGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGAt                                                                                              <  1:3822352/149‑1 (MQ=255)
                                                          gTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTACGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCg                                                                                            >  2:2207111/1‑149 (MQ=255)
                                                                gggTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttt                                                                                      >  2:4155748/1‑149 (MQ=255)
                                                                 ggTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAtttt                                                                                     <  1:1954126/149‑1 (MQ=255)
                                                                  gTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttttt                                                                                    >  1:3104248/1‑149 (MQ=255)
                                                                  gTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCAttttt                                                                                    >  1:711813/1‑149 (MQ=255)
                                                                      tttCTACCCGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGcc                                                                                >  2:2151685/1‑149 (MQ=255)
                                                                         cTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATc                                                                             <  2:1446609/149‑1 (MQ=255)
                                                                          tACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCa                                                                            <  2:2293262/149‑1 (MQ=255)
                                                                           aCCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAg                                                                           >  1:3609427/1‑149 (MQ=255)
                                                                              ggAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCtt                                                                        >  2:3020932/1‑149 (MQ=255)
                                                                               gAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCttt                                                                       <  2:1521237/149‑1 (MQ=255)
                                                                                 aTGGTGTCGTTGCCGGCCAGCTCCGCTGCGGGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTc                                                                     <  2:2221662/149‑1 (MQ=255)
                                                                                 aTGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTc                                                                     <  1:4139668/149‑1 (MQ=255)
                                                                                 aTGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTc                                                                     <  2:1561910/149‑1 (MQ=255)
                                                                                    gtgtCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGTTGTTGGTTGGCGCGTGCTTCATCACCACTTTGGTGCGCTGCTTGGGCAGCCAGTAGATAGTTCTCAACCGGTTGTTCCGCGTGATCTGCATTTTTAGCCATCAGCTTTTCAAc                                                                  <  1:532404/149‑1 (MQ=255)
                                                                                        cGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGc                                                              <  1:4527310/149‑1 (MQ=255)
                                                                                        cGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGc                                                              <  2:4560167/149‑1 (MQ=255)
                                                                                           tGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGa                                                           <  1:766777/149‑1 (MQ=255)
                                                                                            gCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGAt                                                          <  1:1164092/149‑1 (MQ=255)
                                                                                                gCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATc                                                      >  2:2515145/1‑149 (MQ=255)
                                                                                                  cAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGc                                                    <  2:3229387/149‑1 (MQ=255)
                                                                                                    gCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCt                                                  >  2:441462/1‑149 (MQ=255)
                                                                                                         gCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGc                                             >  2:628812/1‑149 (MQ=255)
                                                                                                           tGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCa                                           <  1:4048898/149‑1 (MQ=255)
                                                                                                           tGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCa                                           <  1:3899658/149‑1 (MQ=255)
                                                                                                             cgcgTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGt                                         <  2:3639176/149‑1 (MQ=255)
                                                                                                                gTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTc                                      >  1:949529/1‑149 (MQ=255)
                                                                                                                gTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTc                                      <  1:1295824/149‑1 (MQ=255)
                                                                                                                gTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTc                                      <  1:1501675/149‑1 (MQ=255)
                                                                                                                gTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTc                                      <  1:3057161/149‑1 (MQ=255)
                                                                                                                    cAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTcagca                                  >  2:4515742/1‑149 (MQ=255)
                                                                                                                     aGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTcagcag                                 <  1:4546232/149‑1 (MQ=255)
                                                                                                                       aTGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAgcg                               >  1:3941717/1‑149 (MQ=255)
                                                                                                                         gttggttgGCGCGTGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGcc                             <  1:3237959/149‑1 (MQ=255)
                                                                                                                         gttggttgGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGcc                             <  1:2853892/149‑1 (MQ=255)
                                                                                                                          ttggttggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCt                            >  1:2013119/1‑149 (MQ=255)
                                                                                                                            ggttggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGt                          <  1:1461474/149‑1 (MQ=255)
                                                                                                                             gttggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGtt                         <  2:4662804/149‑1 (MQ=255)
                                                                                                                             gttggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGtt                         <  1:4673578/149‑1 (MQ=255)
                                                                                                                             gttggCGCGTGCTTCATCACCACGTGGTTGCGCTGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGtt                         <  2:4061542/149‑1 (MQ=255)
                                                                                                                               tggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCg                       <  1:497601/149‑1 (MQ=255)
                                                                                                                                ggCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCgg                      >  2:2334232/1‑149 (MQ=255)
                                                                                                                                  cgcgTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTc                    <  1:1135947/149‑1 (MQ=255)
                                                                                                                                   gcgTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCt                   >  1:4433055/1‑149 (MQ=255)
                                                                                                                                   gcgTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCt                   <  2:3484404/149‑1 (MQ=255)
                                                                                                                                      tGCTTCATCACCACTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCt                >  1:3453416/1‑149 (MQ=255)
                                                                                                                                       gCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCtt               >  2:2118343/1‑149 (MQ=255)
                                                                                                                                           catcaCCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGt           >  1:2253436/1‑149 (MQ=255)
                                                                                                                                             tcaCCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGc         >  2:3123978/1‑149 (MQ=255)
                                                                                                                                             tcaCCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGc         >  1:1865706/1‑149 (MQ=255)
                                                                                                                                              caccacGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCa        <  1:336292/149‑1 (MQ=255)
                                                                                                                                               accacTTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCAc       <  1:3490050/149‑1 (MQ=255)
                                                                                                                                               accacGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCAc       >  1:2374437/1‑149 (MQ=255)
                                                                                                                                                  acGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCACGGc    <  2:4752465/149‑1 (MQ=255)
                                                                                                                                                    gttgttGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCACGGCGa  <  2:3284308/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCAGCTTATCCACGCCGTGGCGTGCAGCATGGTTTTCATTACGGGCCAGTTGCAGACGTACGCGGGTGATTTCTACCGGAATGGTGTCGTTGCCGGCCAGCTCCGCTGCGCGTTCCAGATGTTGGTTGGCGCGTGCTTCATCACCACGTTGTTGCGCGGCTTCGGCAGCCAGTAGATAGTTCACCACCGGTTGTTCCGCGTGATCGGCATTTTTCGCCATCAGCTTTTCAACTTGCTGATAATCGCCTTCCGCCAGTTTCAGCAGCGCCTGTTCGGTCTGCTTACGTGCACGGCGA  >  CP009273/3980704‑3981000

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: