Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,114,911 A→T 100% L42Q (CTG→CAG)  wcaJ ← colanic biosynthesis UDP‑glucose lipid carrier transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,114,9110AT98.7% 1003.9 / ‑0.1 296L42Q (CTG→CAG) wcaJcolanic biosynthesis UDP‑glucose lipid carrier transferase
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base T (125/165);  total (128/168)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 296 total aligned reads displayed.

TTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCC  >  CP009273/2114765‑2115056
                                                                                                                                                  |                                                                                                                                                 
ttAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGt                                                                                                                                                 <  1:950430/149‑1 (MQ=255)
ttAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGt                                                                                                                                                 <  2:4237424/149‑1 (MQ=255)
  aaGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTcc                                                                                                                                               >  2:3693680/1‑149 (MQ=255)
    ggTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGACTGACTGTCCGc                                                                                                                                             >  2:3756118/1‑149 (MQ=255)
    ggTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGc                                                                                                                                             <  1:4517929/149‑1 (MQ=255)
     gTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCt                                                                                                                                            <  2:3142743/149‑1 (MQ=255)
     gTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCt                                                                                                                                            <  1:1979355/149‑1 (MQ=255)
      tCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTCGAGGAATGACAGTCAGCTg                                                                                                                                           >  2:2228960/1‑149 (MQ=255)
         aGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACt                                                                                                                                        >  2:545258/1‑149 (MQ=255)
                    aaCAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAg                                                                                                                             >  2:68884/1‑149 (MQ=255)
                     aCAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGc                                                                                                                            >  1:4716337/1‑149 (MQ=255)
                     aCAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGc                                                                                                                            >  1:573459/1‑149 (MQ=255)
                       aGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCa                                                                                                                          <  1:177754/149‑1 (MQ=255)
                        gggCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCAt                                                                                                                         >  2:4526219/1‑149 (MQ=255)
                        gggCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCAt                                                                                                                         >  1:2548869/1‑149 (MQ=255)
                        gggCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCAt                                                                                                                         >  1:67014/1‑149 (MQ=255)
                         ggCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATa                                                                                                                        <  2:3400516/149‑1 (MQ=255)
                         ggCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATa                                                                                                                        >  2:2055348/1‑149 (MQ=255)
                            aaaTTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTc                                                                                                                     <  1:2393582/149‑1 (MQ=255)
                             aaTTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTcc                                                                                                                    <  2:1232638/149‑1 (MQ=255)
                                tCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGc                                                                                                                 <  1:4391071/149‑1 (MQ=255)
                                  tGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCaa                                                                                                               >  1:3454152/1‑149 (MQ=255)
                                  tGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCaa                                                                                                               >  1:3188480/1‑149 (MQ=255)
                                   gTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCaaa                                                                                                              >  1:1820812/1‑149 (MQ=255)
                                    tCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAAc                                                                                                             >  1:3278860/1‑149 (MQ=255)
                                      gCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACat                                                                                                           >  1:4211840/1‑149 (MQ=255)
                                          cccGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACatgatg                                                                                                       <  2:464288/149‑1 (MQ=255)
                                            cGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGt                                                                                                     <  1:3720311/149‑1 (MQ=255)
                                             gAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTg                                                                                                    >  1:2159341/1‑149 (MQ=255)
                                               acacCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGat                                                                                                  >  1:4442013/1‑149 (MQ=255)
                                                cacCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGata                                                                                                 <  2:1183424/149‑1 (MQ=255)
                                                  ccGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATc                                                                                               <  1:3516189/149‑1 (MQ=255)
                                                  ccGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATc                                                                                               <  2:42108/149‑1 (MQ=255)
                                                   cgcgCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCt                                                                                              <  2:3832207/149‑1 (MQ=255)
                                                    gcgcCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTg                                                                                             <  2:302333/149‑1 (MQ=255)
                                                       ccATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGaaa                                                                                          >  2:4375557/1‑149 (MQ=255)
                                                        cATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGaaaa                                                                                         >  1:3818424/1‑149 (MQ=255)
                                                         aTGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAg                                                                                        >  1:3710088/1‑149 (MQ=255)
                                                            aGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGt                                                                                     >  2:2364948/1‑149 (MQ=255)
                                                            aGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGt                                                                                     <  1:4747661/149‑1 (MQ=255)
                                                               gATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGc                                                                                  <  1:183037/149‑1 (MQ=255)
                                                                aTAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCa                                                                                 <  1:3798993/149‑1 (MQ=255)
                                                                   aaaaTCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCa                                                                              >  1:2875512/1‑149 (MQ=255)
                                                                     aaTCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATa                                                                            <  2:2647911/149‑1 (MQ=255)
                                                                       tCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATaga                                                                          >  1:1136368/1‑149 (MQ=255)
                                                                         ggTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAgaga                                                                        <  2:2591538/149‑1 (MQ=255)
                                                                         ggTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAgaga                                                                        <  2:2570497/149‑1 (MQ=255)
                                                                         ggTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAgaga                                                                        <  1:1925242/149‑1 (MQ=255)
                                                                         ggTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAgaga                                                                        >  2:1005058/1‑149 (MQ=255)
                                                                          gTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGAt                                                                       <  2:408655/149‑1 (MQ=255)
                                                                           tGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGAtt                                                                      >  2:2163906/1‑149 (MQ=255)
                                                                             aTGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTaa                                                                    <  1:3238271/149‑1 (MQ=255)
                                                                              tGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCACCAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAAc                                                                   <  1:2293812/149‑1 (MQ=255)
                                                                                 cgccCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGAt                                                                <  2:1555424/149‑1 (MQ=255)
                                                                                  gccCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATg                                                               >  1:2884558/1‑149 (MQ=255)
                                                                                  gccCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATg                                                               >  2:554293/1‑149 (MQ=255)
                                                                                  gccCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGATGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGACAAGCGTTGCACAATAGAGATTAACGATg                                                               >  1:406238/1‑149 (MQ=255)
                                                                                      aGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCAtt                                                           <  1:3748109/149‑1 (MQ=255)
                                                                                      aGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCAtt                                                           >  1:2326263/1‑149 (MQ=255)
                                                                                       gCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTg                                                          <  2:3322509/149‑1 (MQ=255)
                                                                                        cATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTgg                                                         <  1:3627590/149‑1 (MQ=255)
                                                                                         aTCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGt                                                        <  2:1586834/149‑1 (MQ=255)
                                                                                         aTCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGt                                                        <  2:1671905/149‑1 (MQ=255)
                                                                                          tCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGtt                                                       <  2:1878592/149‑1 (MQ=255)
                                                                                           cTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGttt                                                      <  2:2159341/149‑1 (MQ=255)
                                                                                            tGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGtttt                                                     <  2:3446755/149‑1 (MQ=255)
                                                                                              gAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTtcg                                                   >  2:762882/1‑149 (MQ=255)
                                                                                               aaCACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTtcgc                                                  <  1:3102970/149‑1 (MQ=255)
                                                                                                  accaccAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTtcgctcg                                               >  1:4359555/1‑149 (MQ=255)
                                                                                                   ccaccaGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTtcgctcgc                                              >  2:1551300/1‑149 (MQ=255)
                                                                                                       caGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTtcgctcgctcgc                                          <  2:4422697/149‑1 (MQ=255)
                                                                                                           gTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCtt                                      >  1:4241353/1‑149 (MQ=255)
                                                                                                             aaTCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCtttt                                    <  2:1427375/149‑1 (MQ=255)
                                                                                                               tCAGCGCCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCtttttt                                  <  1:951402/149‑1 (MQ=255)
                                                                                                                  gcgcCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGa                               >  1:1685401/1‑149 (MQ=255)
                                                                                                                  gcgcCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGa                               >  1:1230358/1‑149 (MQ=255)
                                                                                                                   cgcCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGAt                              <  2:2811258/149‑1 (MQ=255)
                                                                                                                    gcCACCAACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGAtt                             <  1:1041070/149‑1 (MQ=255)
                                                                                                                     ccaccaACAGGTGCATGTAGAGGCATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACCTGATGGTGATATCTGAACAGCGTTGCACCATAGCGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGAttt                            <  1:869804/149‑1 (MQ=255)
                                                                                                                     ccaccaACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGAttt                            >  2:318928/1‑149 (MQ=255)
                                                                                                                      caccaACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTg                           <  1:2179706/149‑1 (MQ=255)
                                                                                                                      caccaACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTg                           >  2:3091922/1‑149 (MQ=255)
                                                                                                                        ccaACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTc                         <  2:384675/149‑1 (MQ=255)
                                                                                                                         caACAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCa                        <  1:2810833/149‑1 (MQ=255)
                                                                                                                          aaCAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCAt                       >  2:2003890/1‑149 (MQ=255)
                                                                                                                          aaCAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCAt                       <  2:1798591/149‑1 (MQ=255)
                                                                                                                          aaCAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCAt                       <  2:4779557/149‑1 (MQ=255)
                                                                                                                           aCAGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATc                      >  1:1412531/1‑149 (MQ=255)
                                                                                                                             aGGTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCgt                    <  2:1377317/149‑1 (MQ=255)
                                                                                                                              ggTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCgtt                   >  1:2041174/1‑149 (MQ=255)
                                                                                                                               gTGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCgttg                  >  2:2743080/1‑149 (MQ=255)
                                                                                                                                tGCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCgttgt                 <  2:4726309/149‑1 (MQ=255)
                                                                                                                                 gCATGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCgttgtt                <  1:2944845/149‑1 (MQ=255)
                                                                                                                                    tGTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCt             <  1:3139351/149‑1 (MQ=255)
                                                                                                                                     gTAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTg            <  1:737925/149‑1 (MQ=255)
                                                                                                                                      tAGAGGAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGt           <  2:4649676/149‑1 (MQ=255)
                                                                                                                                          ggAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGttatt       <  2:1136368/149‑1 (MQ=255)
                                                                                                                                           gAATGACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCGCGCGCTTTTTTAGATTTGTGATCGTTGTTCCTGttatta      >  1:2128344/1‑149 (MQ=255)
                                                                                                                                           gAATGACTGTCCGCGGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGttatta      <  1:1374251/149‑1 (MQ=255)
                                                                                                                                               gACTGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGccc  >  1:3829802/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTAAGGTCCAGTTTTGTAGCAACAGGGCAAATTCTGTCGCTGCCCGAACACCGCGCCATGAGCGATAAAAATCGGTGATGCCGCCCAGCATCTGGAACACCACCAGCGTAATCAGCGCCACCAACAGGTGCATGTAGAGGAATGACAGTCCGCTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCC  >  CP009273/2114765‑2115056

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: