Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,107,728 C→A 23.3% E27D (GAG→GAT wcaN ← putative regulatory subunit for GalU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,107,7280CA23.3% 538.6 / 119.2 223E27D (GAG→GATwcaNputative regulatory subunit for GalU
Reads supporting (aligned to +/- strand):  ref base C (87/84);  new base A (29/23);  total (116/107)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 239 total aligned reads displayed.

GAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATT  >  CP009273/2107581‑2107875
                                                                                                                                                   |                                                                                                                                                   
gAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtc                                                                                                                                                     >  1:1022814/1‑148 (MQ=255)
gAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtc                                                                                                                                                     >  1:1023395/1‑148 (MQ=255)
gAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCAtc                                                                                                                                                     <  1:3706758/148‑1 (MQ=255)
   tGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATctct                                                                                                                                                  >  1:2496512/1‑148 (MQ=255)
   tGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATctct                                                                                                                                                  <  2:208349/148‑1 (MQ=255)
   tGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATctct                                                                                                                                                  <  2:2137597/148‑1 (MQ=255)
    gATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCtt                                                                                                                                                 <  2:2454825/148‑1 (MQ=255)
      ttCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTg                                                                                                                                               >  1:2613198/1‑148 (MQ=255)
          aaCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTa                                                                                                                                           <  1:1302682/148‑1 (MQ=255)
          aaCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTa                                                                                                                                           >  2:867734/1‑148 (MQ=255)
           aCTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTAt                                                                                                                                          <  1:3299060/148‑1 (MQ=255)
              cATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGc                                                                                                                                       >  2:2566003/1‑148 (MQ=255)
               aTAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGcc                                                                                                                                      <  2:112873/148‑1 (MQ=255)
                tAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCt                                                                                                                                     <  2:1989912/148‑1 (MQ=255)
                    agGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGt                                                                                                                                 <  2:3319200/148‑1 (MQ=255)
                    agGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGt                                                                                                                                 >  1:895861/1‑148 (MQ=255)
                    agGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGt                                                                                                                                 >  1:3316998/1‑148 (MQ=255)
                         tCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTggcag                                                                                                                            >  2:165939/1‑148 (MQ=255)
                             aGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCaa                                                                                                                        >  2:3198991/1‑148 (MQ=255)
                                ggTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACat                                                                                                                     <  1:165939/148‑1 (MQ=255)
                                 gTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACata                                                                                                                    <  2:1207060/148‑1 (MQ=255)
                                  ttttCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACatat                                                                                                                   <  2:1525799/148‑1 (MQ=255)
                                          cgcgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCg                                                                                                           >  2:443337/1‑148 (MQ=255)
                                            cgTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCgag                                                                                                         <  2:1552644/148‑1 (MQ=255)
                                             gTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCgaga                                                                                                        <  1:296768/148‑1 (MQ=255)
                                              ttCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGAc                                                                                                       >  2:53438/1‑148 (MQ=255)
                                              ttCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGAc                                                                                                       >  2:3710195/1‑148 (MQ=255)
                                               tCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGAcc                                                                                                      >  1:2145123/1‑148 (MQ=255)
                                                  tGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGc                                                                                                   <  2:741910/148‑1 (MQ=255)
                                                  tGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGc                                                                                                   >  1:2339978/1‑148 (MQ=255)
                                                       gcgTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAg                                                                                              <  1:1140142/148‑1 (MQ=255)
                                                         gTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGa                                                                                            <  2:2872711/148‑1 (MQ=255)
                                                         gTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGa                                                                                            <  1:1254017/148‑1 (MQ=255)
                                                         gTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGa                                                                                            <  1:734440/148‑1 (MQ=255)
                                                          tGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGaa                                                                                           <  1:1188239/148‑1 (MQ=255)
                                                           gAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAAt                                                                                          >  1:984627/1‑148 (MQ=255)
                                                           gAGTTACCAGGAGGATTTCTTTGATCCCTGAAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAAt                                                                                          >  2:1120914/1‑148 (MQ=255)
                                                              ttACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAAc                                                                                       <  2:1774278/148‑1 (MQ=255)
                                                                 ccAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGc                                                                                    >  2:956405/1‑148 (MQ=255)
                                                                   aggaggATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtt                                                                                  >  1:2258246/1‑148 (MQ=255)
                                                                   aggaggATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtt                                                                                  <  2:3572869/148‑1 (MQ=255)
                                                                   aggaggATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtt                                                                                  <  2:1205691/148‑1 (MQ=255)
                                                                     gaggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCtttt                                                                                >  1:2878238/1‑148 (MQ=255)
                                                                       ggaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTaa                                                                              <  2:3129513/148‑1 (MQ=255)
                                                                        gaTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTaaa                                                                             >  2:2520322/1‑148 (MQ=255)
                                                                         aTTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAAt                                                                            <  2:2262357/148‑1 (MQ=255)
                                                                          tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAAtt                                                                           >  2:229350/1‑148 (MQ=255)
                                                                          tttCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAAtt                                                                           >  2:2357840/1‑148 (MQ=255)
                                                                            tCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCg                                                                         <  2:266956/148‑1 (MQ=255)
                                                                             cTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGt                                                                        >  2:2921388/1‑148 (MQ=255)
                                                                              tttGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTc                                                                       <  1:137644/148‑1 (MQ=255)
                                                                               ttGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCa                                                                      >  1:2024741/1‑148 (MQ=255)
                                                                               ttGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCa                                                                      >  2:1479087/1‑148 (MQ=255)
                                                                                  aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCttt                                                                                 >  2:52715/1‑134 (MQ=255)
                                                                                  aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCttt                                                                                 <  1:52715/134‑1 (MQ=255)
                                                                                  aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCatta                                                                   >  2:2430953/1‑148 (MQ=255)
                                                                                  aTCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCatta                                                                   >  1:1089854/1‑148 (MQ=255)
                                                                                   tCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattat                                                                  >  2:3043279/1‑148 (MQ=255)
                                                                                    cccTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattatt                                                                 >  2:1588317/1‑148 (MQ=255)
                                                                                    cccTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCattatt                                                                 >  2:932732/1‑148 (MQ=255)
                                                                                     ccTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTAttt                                                                <  1:1749471/148‑1 (MQ=255)
                                                                                     ccTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTAttt                                                                >  2:2853691/1‑148 (MQ=255)
                                                                                       tGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCa                                                              >  2:1998748/1‑148 (MQ=255)
                                                                                        gCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCAt                                                             >  2:1723818/1‑148 (MQ=255)
                                                                                        gCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCAt                                                             <  2:1233508/148‑1 (MQ=255)
                                                                                            ccACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCAc                                                         <  2:1829956/148‑1 (MQ=255)
                                                                                            ccACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCAc                                                         >  2:1462059/1‑148 (MQ=255)
                                                                                            ccACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCAc                                                         >  1:1921770/1‑148 (MQ=255)
                                                                                              acaATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCt                                                       >  1:3497565/1‑148 (MQ=255)
                                                                                               caATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTg                                                      >  1:2821963/1‑148 (MQ=255)
                                                                                                   ctcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaaa                                                  >  1:2069973/1‑148 (MQ=255)
                                                                                                    tcGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTaaaa                                                 <  1:1006191/148‑1 (MQ=255)
                                                                                                      gTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATg                                               <  1:84824/148‑1 (MQ=255)
                                                                                                         aacaaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGtt                                            <  1:2990266/148‑1 (MQ=255)
                                                                                                          acaaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTg                                           >  2:2410732/1‑148 (MQ=255)
                                                                                                          acaaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTg                                           <  2:2496512/148‑1 (MQ=255)
                                                                                                           caaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGc                                          >  1:3069616/1‑148 (MQ=255)
                                                                                                           caaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGc                                          >  1:1164360/1‑148 (MQ=255)
                                                                                                           caaTGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGc                                          >  1:2721886/1‑148 (MQ=255)
                                                                                                               gTACTGAATCATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGaa                                      >  2:3491685/1‑148 (MQ=255)
                                                                                                                 aCTGATTGATTGGCTTGTCGACGATTGGTTTCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGTGACCCGCTACAGGAATAACTGCTTTTATATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAAtt                                    <  2:3455516/148‑1 (MQ=255)
                                                                                                                   tGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATtat                                  <  2:2759428/148‑1 (MQ=255)
                                                                                                                    gAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGTTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAACTGGTTGCTGAATtata                                 <  1:3125440/148‑1 (MQ=255)
                                                                                                                     aaTCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAg                                <  1:3717620/148‑1 (MQ=255)
                                                                                                                      aTCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGc                               >  1:1926763/1‑148 (MQ=255)
                                                                                                                       tCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACATGTAAAATGGTTGATGAATTATAGCt                              >  2:2302070/1‑148 (MQ=255)
                                                                                                                        cATTGGCTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCtt                             <  1:452152/148‑1 (MQ=255)
                                                                                                                         aTTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTg                            <  2:1023395/148‑1 (MQ=255)
                                                                                                                           tGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGtt                          >  2:3374725/1‑148 (MQ=255)
                                                                                                                           tGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGtt                          >  1:671078/1‑148 (MQ=255)
                                                                                                                             gCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCg                        <  2:1878470/148‑1 (MQ=255)
                                                                                                                             gCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCg                        <  1:3233794/148‑1 (MQ=255)
                                                                                                                              cTTGTCGACGATTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGa                       <  1:2765394/148‑1 (MQ=255)
                                                                                                                                      cgaTTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTc               >  1:1919245/1‑148 (MQ=255)
                                                                                                                                       gaTTGGTAGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCg              >  2:2648153/1‑148 (MQ=255)
                                                                                                                                          tGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCa           >  1:1791048/1‑148 (MQ=255)
                                                                                                                                             tAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCa        >  1:3584706/1‑148 (MQ=255)
                                                                                                                                              aGCATATCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCAt       >  2:2191047/1‑148 (MQ=255)
                                                                                                                                                   ctctTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAAtt  >  1:3737173/1‑148 (MQ=255)
                                                                                                                                                   ctctTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAAtt  <  2:823561/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAGTGATTCTAACTCATAAGAGGTGTCGAAGTGGTTTTCGACCGCGTTCTTGGACGCGTGAGTTACCAGGAGGATTTCTTTGATCCCTGCAGCCACAATCTCGTCAACAATGTACTGAATCATTGGCTTGTCGACGATTGGTAGCATCTCTTTGGGTATCGCCTTAGTGGCAGGCAACATATGCATCCCGAGACCCGCTACAGGAATAACTGCTTTTAAATTCGTCATTATTTCATCCACCTGTAAAATGGTTGCTGAATTATAGCTTGTTCGATTTTTTTCGCCAGCATCAATT  >  CP009273/2107581‑2107875

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: