Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,491,629 C→A 8.3% P148Q (CCG→CAG)  cysG → fused siroheme synthase 1,3‑dimethyluroporphyriongen III dehydrogenase and siroheme ferrochelatase/uroporphyrinogen methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,491,6290CA 8.3% 1029.5 / 51.5 287P148Q (CCG→CAG) cysGfused siroheme synthase 1,3‑dimethyluroporphyriongen III dehydrogenase and siroheme ferrochelatase/uroporphyrinogen methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (112/151);  new base A (8/16);  total (120/167)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 288 total aligned reads displayed.

TTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGCTGGCAAACAACGATCAGA  >  CP009273/3491481‑3491778
                                                                                                                                                    |                                                                                                                                                     
ttctGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCg                                                                                                                                                      >  2:2385550/1‑150 (MQ=255)
  ctGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGTCCGCTCACCGCTCATGGAAGCGGTCTCCAATGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTAAATCACAGATGCCGtt                                                                                                                                                    >  2:4860399/1‑150 (MQ=255)
  ctGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGtt                                                                                                                                                    <  2:4589135/150‑1 (MQ=255)
    gTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTAc                                                                                                                                                  >  1:3046099/1‑150 (MQ=255)
     tAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGATACa                                                                                                                                                 >  2:2317871/1‑150 (MQ=255)
        cGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCt                                                                                                                                              <  1:5144895/150‑1 (MQ=255)
        cGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGGTACATCt                                                                                                                                              >  2:4759212/1‑150 (MQ=255)
            gTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGc                                                                                                                                          >  2:5115379/1‑150 (MQ=255)
                aTGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAgg                                                                                                                                      <  2:3997904/150‑1 (MQ=255)
                  gcgcCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTa                                                                                                                                    >  1:1934102/1‑150 (MQ=255)
                     ccGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCg                                                                                                                                 <  1:4656519/150‑1 (MQ=255)
                     ccGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCg                                                                                                                                 >  1:2944839/1‑150 (MQ=255)
                      cGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGa                                                                                                                                >  1:1342518/1‑150 (MQ=255)
                         aaGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAAt                                                                                                                             <  1:4348655/150‑1 (MQ=255)
                           gccgccAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATAc                                                                                                                           <  2:2627665/150‑1 (MQ=255)
                            ccgccAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACg                                                                                                                          >  1:4554827/1‑150 (MQ=255)
                             cgccAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGc                                                                                                                         >  2:1740157/1‑150 (MQ=255)
                              gccAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGcc                                                                                                                        >  2:3737746/1‑150 (MQ=255)
                              gccAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGcc                                                                                                                        <  1:3611653/150‑1 (MQ=255)
                              gccAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGcc                                                                                                                        >  2:2529110/1‑150 (MQ=255)
                               ccAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCACCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCg                                                                                                                       >  1:3726743/1‑150 (MQ=255)
                                cAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCgg                                                                                                                      >  1:4188277/1‑150 (MQ=255)
                                    tttATTCTGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCaa                                                                                                                  <  1:4914803/150‑1 (MQ=255)
                                       attatGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTa                                                                                                               <  1:2759088/150‑1 (MQ=255)
                                       attatGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTa                                                                                                               <  2:1263029/150‑1 (MQ=255)
                                        ttatGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTAc                                                                                                              <  1:4864849/150‑1 (MQ=255)
                                           tGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTAcgcg                                                                                                           <  2:2324685/150‑1 (MQ=255)
                                             ccGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCggg                                                                                                         <  2:64837/150‑1 (MQ=255)
                                              cGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCTCCTCTCCGGTTCTGGCACGCCTGTTGCGCGAACAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGc                                                                                                        <  2:843856/150‑1 (MQ=255)
                                              cGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGc                                                                                                        <  2:3157109/150‑1 (MQ=255)
                                               gTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCCGGTAGCGAAATACGCCGGGCAATTACGCGGGCg                                                                                                       >  2:5074882/1‑150 (MQ=255)
                                                tCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGa                                                                                                      >  1:1019002/1‑150 (MQ=255)
                                                 cGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAg                                                                                                     >  1:2078809/1‑150 (MQ=255)
                                                  gATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGt                                                                                                    >  1:608504/1‑150 (MQ=255)
                                                   attagtGACCGCTCACCGCTCAGGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTg                                                                                                   <  2:2716124/150‑1 (MQ=255)
                                                    ttattGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGa                                                                                                  >  1:2242364/1‑150 (MQ=255)
                                                    ttattGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGa                                                                                                  >  1:3654429/1‑150 (MQ=255)
                                                        tGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAaaca                                                                                              <  1:4611124/150‑1 (MQ=255)
                                                          accgctcaccgctcATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCAGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAaacaac                                                                                            <  1:1845945/150‑1 (MQ=255)
                                                          accgctcaccgctcATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCAGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAaacaac                                                                                            <  1:4452601/150‑1 (MQ=255)
                                                          accgctcaccgctcATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCAGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAaacaac                                                                                            <  2:1203914/150‑1 (MQ=255)
                                                           ccgctcaccgctcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAacaaca                                                                                           >  2:2924292/1‑150 (MQ=255)
                                                            cgctcaccgctcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAg                                                                                          >  2:358827/1‑150 (MQ=255)
                                                            cgctcaccgctcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAg                                                                                          >  1:593471/1‑150 (MQ=255)
                                                                caccgctcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCg                                                                                      >  1:500747/1‑150 (MQ=255)
                                                                    gctcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCAc                                                                                  <  2:4125916/150‑1 (MQ=255)
                                                                     ctcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACg                                                                                 <  1:4765913/150‑1 (MQ=255)
                                                                      tcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGa                                                                                >  1:4322155/1‑150 (MQ=255)
                                                                      tcaTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCAGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGa                                                                                >  2:1166026/1‑150 (MQ=255)
                                                                       caTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGAt                                                                               >  1:1999883/1‑150 (MQ=255)
                                                                        aTGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATg                                                                              <  2:4386396/150‑1 (MQ=255)
                                                                         tGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATgg                                                                             <  1:4218696/150‑1 (MQ=255)
                                                                          ggTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATggg                                                                            >  1:2593978/1‑150 (MQ=255)
                                                                           gTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGTACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGt                                                                           >  2:2832154/1‑150 (MQ=255)
                                                                           gTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACTCCTGTTGCGCGACAAACTTGAATCACGGCTGCCGTTACGTCTGGGCCAGGTAGCGAATTACGCCGGGCATGTACGCGGGCGAGTGTACCTCCAGGGGGCCACGAAGGGt                                                                           >  2:3668433/1‑150 (MQ=255)
                                                                              gCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAg                                                                        <  2:3863433/150‑1 (MQ=255)
                                                                                  tctcCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTc                                                                    <  1:1635643/150‑1 (MQ=255)
                                                                                  tctcCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTc                                                                    <  1:5191737/150‑1 (MQ=255)
                                                                                   ctcctcTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCg                                                                   >  2:576001/1‑150 (MQ=255)
                                                                                    tcctcTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGc                                                                  >  1:711438/1‑150 (MQ=255)
                                                                                     cctcTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGcc                                                                 <  2:119677/150‑1 (MQ=255)
                                                                                     cctcTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGcc                                                                 <  1:435575/150‑1 (MQ=255)
                                                                                        ctGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGtt                                                              <  2:527184/150‑1 (MQ=255)
                                                                                        ctGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGtt                                                              >  1:5218778/1‑150 (MQ=255)
                                                                                           gcggcACCTCTCCGGTTCTGGCACGCCTGGTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCt                                                           <  2:2185607/150‑1 (MQ=255)
                                                                                            cggcACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTg                                                          <  2:3870203/150‑1 (MQ=255)
                                                                                             ggcACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTgg                                                         >  2:2032649/1‑150 (MQ=255)
                                                                                             ggcACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTgg                                                         >  1:31210/1‑150 (MQ=255)
                                                                                             ggcACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATAACTGCTGCCGTTACAACTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAAAAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTgt                                                         >  2:4091485/1‑149 (MQ=255)
                                                                                                  ctctCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGaa                                                    >  1:2613866/1‑150 (MQ=255)
                                                                                                    ctcCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAAt                                                  >  1:1829625/1‑150 (MQ=255)
                                                                                                    ctcCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCAGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAAt                                                  <  1:3798145/150‑1 (MQ=255)
                                                                                                       cGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGt                                               <  2:1088375/150‑1 (MQ=255)
                                                                                                       cGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGt                                               <  1:3334634/150‑1 (MQ=255)
                                                                                                        ggTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGtt                                              <  2:1093667/150‑1 (MQ=255)
                                                                                                             tGGCACGCCTGTTTCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTTGGCCAGGTAGCGAAATACGACGGGCAATTACTCGTGCGAGTGATACAACAGTTCGCCAAGATGGGTGAGCGTCTCCGTTTCTGGGAGAAATTGTTCGTTa                                         <  1:1762234/150‑1 (MQ=255)
                                                                                                             tGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTa                                         <  2:3102672/150‑1 (MQ=255)
                                                                                                             tGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTa                                         <  2:3139044/150‑1 (MQ=255)
                                                                                                             tGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTa                                         <  1:1058104/150‑1 (MQ=255)
                                                                                                              ggCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTaa                                        >  1:775661/1‑150 (MQ=255)
                                                                                                               gCACGCCTGTTGCGCGAAAAACTTGACTCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAAc                                       <  2:3313335/150‑1 (MQ=255)
                                                                                                               gCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCAGGTACACCTGGGCCCCCCCCCCAAACACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAAc                                       <  1:1893279/150‑1 (MQ=255)
                                                                                                                cACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAACTTGTTCGTTAACg                                      >  2:1369748/1‑150 (MQ=255)
                                                                                                                cACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACg                                      >  2:1970128/1‑150 (MQ=255)
                                                                                                                  cGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGAc                                    >  2:3067244/1‑150 (MQ=255)
                                                                                                                  cGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCAGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGAc                                    <  1:4543730/150‑1 (MQ=255)
                                                                                                                    ccTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGTAGAAATTGTTCGTTAACGACCg                                  >  2:2549700/1‑150 (MQ=255)
                                                                                                                     cTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGc                                 >  1:3563171/1‑150 (MQ=255)
                                                                                                                       gTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCt                               <  2:2113821/150‑1 (MQ=255)
                                                                                                                       gTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCt                               <  2:1457815/150‑1 (MQ=255)
                                                                                                                         tGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTgg                             <  2:1904972/150‑1 (MQ=255)
                                                                                                                           cgcgAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGgcg                           >  2:4240379/1‑150 (MQ=255)
                                                                                                                             cgAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCa                         >  1:946877/1‑150 (MQ=255)
                                                                                                                               aaaaaCTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGt                       >  2:1303481/1‑150 (MQ=255)
                                                                                                                                  aaCTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGc                    <  2:1408891/150‑1 (MQ=255)
                                                                                                                                   aCTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGCt                   >  1:3188428/1‑150 (MQ=255)
                                                                                                                                   aCTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGCt                   >  1:1582405/1‑150 (MQ=255)
                                                                                                                                            aCTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGCTGGCAaacaa          >  2:4911052/1‑150 (MQ=255)
                                                                                                                                             ctgctgCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGCTGGCAaacaac         <  1:5159307/150‑1 (MQ=255)
                                                                                                                                             ctgctgCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGCTGGCAaacaac         >  1:3966904/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTCGATTATTGACCGCTCACCGCTCATGGTAGCGGTCTCCTCTGGCGGCACCTCTCCGGTTCTGGCACGCCTGTTGCGCGAAAAACTTGAATCACTGCTGCCGTTACATCTGGGCCAGGTAGCGAAATACGCCGGGCAATTACGCGGGCGAGTGAAACAACAGTTCGCCACGATGGGTGAGCGTCGCCGTTTCTGGGAGAAATTGTTCGTTAACGACCGCCTGGCGCAGTCGCTGGCAAACAACGATCAGA  >  CP009273/3491481‑3491778

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: