Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,113,761 C→A 5.2% E425D (GAG→GAT 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,113,7610CA 5.2% 1020.7 / 20.7 269E425D (GAG→GATwcaJcolanic biosynthesis UDP‑glucose lipid carrier transferase
Reads supporting (aligned to +/- strand):  ref base C (157/98);  new base A (3/11);  total (160/109)
Fisher's exact test for biased strand distribution p-value = 4.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 284 total aligned reads displayed.

CGCTGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCG  >  CP009273/2113612‑2113909
                                                                                                                                                     |                                                                                                                                                    
 gcTGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTtct                                                                                                                                                     >  1:1612763/1‑150 (MQ=255)
 gcTGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTtct                                                                                                                                                     >  1:1612827/1‑150 (MQ=255)
 gcTGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTtct                                                                                                                                                     >  1:1613067/1‑150 (MQ=255)
 gcTGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTtct                                                                                                                                                     >  1:40319/1‑150 (MQ=255)
   tGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTcc                                                                                                                                                   <  2:659630/150‑1 (MQ=255)
   tGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTcc                                                                                                                                                   >  2:1148727/1‑150 (MQ=255)
      tttttttttCACGTAAGCTAATATCAATATGCCTCTTTGTTAACGAAACCTTTGAACTCCGTCAGTAAATCGATTTTGATATCGAACCAGACGCTCCAGTCGCTGATGTACTAAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGc                                                                                                                                                <  2:4446316/147‑1 (MQ=255)
      tGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGc                                                                                                                                                >  1:3653047/1‑150 (MQ=255)
      tGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGc                                                                                                                                                >  1:525762/1‑150 (MQ=255)
        gTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCgt                                                                                                                                              <  1:3263541/150‑1 (MQ=255)
          tttttCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCgtgt                                                                                                                                            >  2:818913/1‑150 (MQ=255)
          tttttCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCgtgt                                                                                                                                            >  2:3350544/1‑150 (MQ=255)
             ttCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCgg                                                                                                                                         >  1:274445/1‑150 (MQ=255)
              tCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGt                                                                                                                                        >  1:4265479/1‑150 (MQ=255)
              tCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGt                                                                                                                                        >  1:4265236/1‑150 (MQ=255)
               cACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGtt                                                                                                                                       >  2:3903525/1‑150 (MQ=255)
                aCGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGttt                                                                                                                                      >  1:244372/1‑150 (MQ=255)
                 cGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTc                                                                                                                                     <  2:3834074/150‑1 (MQ=255)
                         cATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGcc                                                                                                                             <  2:2705254/150‑1 (MQ=255)
                           tatCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAg                                                                                                                           <  2:1372682/150‑1 (MQ=255)
                            atCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGc                                                                                                                          >  2:950313/1‑150 (MQ=255)
                            atCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGc                                                                                                                          <  1:1139069/150‑1 (MQ=255)
                                atatGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGtt                                                                                                                      >  2:169976/1‑150 (MQ=255)
                                atatGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGtt                                                                                                                      >  2:3193257/1‑150 (MQ=255)
                                 tatGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTa                                                                                                                     >  1:2957094/1‑150 (MQ=255)
                                   tGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGACAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAAt                                                                                                                   >  2:4507803/1‑150 (MQ=255)
                                     ccGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCt                                                                                                                 <  1:3236769/150‑1 (MQ=255)
                                     ccGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCt                                                                                                                 <  1:3236680/150‑1 (MQ=255)
                                      cGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTg                                                                                                                <  2:484859/150‑1 (MQ=255)
                                      cGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGAAGCGTTTTTCCATTTTCTCCAGAGTGTCGGTTTCGCCGCGCCAGACGTTAATCtt                                                                                                                >  2:3317522/1‑149 (MQ=255)
                                          ttGTTATCGAAACCTTTGAACACCGTCAGGAATACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGcc                                                                                                            <  2:2190754/150‑1 (MQ=255)
                                            gTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCa                                                                                                          >  2:491373/1‑150 (MQ=255)
                                              tAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGc                                                                                                        >  1:1892661/1‑150 (MQ=255)
                                                aCGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCg                                                                                                      <  1:2577169/150‑1 (MQ=255)
                                                 cGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGAGCCAGCCGTTAATCTGCGCCCAGCCgg                                                                                                     >  2:2125395/1‑150 (MQ=255)
                                                 cGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATACGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGAGCCAGCCGTTAATCTGCGCCCAGCCgg                                                                                                     >  2:4368476/1‑150 (MQ=255)
                                                   aaaCCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTa                                                                                                   <  1:168994/150‑1 (MQ=255)
                                                      ccTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATg                                                                                                >  1:4394272/1‑150 (MQ=255)
                                                          tGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGGGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCg                                                                                            >  1:1927381/1‑150 (MQ=255)
                                                              cacCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGttt                                                                                        >  2:3575043/1‑150 (MQ=255)
                                                                ccGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCa                                                                                      >  1:931666/1‑150 (MQ=255)
                                                                 cGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCAc                                                                                     >  2:2956829/1‑150 (MQ=255)
                                                                 cGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCAc                                                                                     >  2:2086676/1‑150 (MQ=255)
                                                                 cGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCAc                                                                                     >  2:2086817/1‑150 (MQ=255)
                                                                    cAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCtt                                                                                  <  2:2832922/150‑1 (MQ=255)
                                                                     aGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTa                                                                                 >  1:316972/1‑150 (MQ=255)
                                                                       gAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATg                                                                               >  1:345731/1‑150 (MQ=255)
                                                                         aaaCGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGgc                                                                             <  2:57682/150‑1 (MQ=255)
                                                                         aaaCGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGgc                                                                             >  2:923908/1‑150 (MQ=255)
                                                                         aaaCGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGgc                                                                             <  2:181385/150‑1 (MQ=255)
                                                                         aaaCGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGgc                                                                             <  2:179390/150‑1 (MQ=255)
                                                                          aaCGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGgcg                                                                            <  2:1223318/150‑1 (MQ=255)
                                                                          aaCGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGgcg                                                                            >  2:1938745/1‑150 (MQ=255)
                                                                            cGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCgca                                                                          >  2:3094068/1‑150 (MQ=255)
                                                                            cGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCgca                                                                          >  2:3692702/1‑150 (MQ=255)
                                                                            cGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCgca                                                                          >  1:1111664/1‑150 (MQ=255)
                                                                            cGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCgca                                                                          >  1:759173/1‑150 (MQ=255)
                                                                               ttttGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCgcagca                                                                       >  1:728287/1‑150 (MQ=255)
                                                                               ttttGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCgcagca                                                                       >  2:2798855/1‑150 (MQ=255)
                                                                                tttGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCAt                                                                      <  2:3159407/150‑1 (MQ=255)
                                                                                 ttGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATg                                                                     <  2:1982762/150‑1 (MQ=255)
                                                                                  tGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGt                                                                    <  1:3045309/150‑1 (MQ=255)
                                                                                  tGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGt                                                                    >  1:3605454/1‑150 (MQ=255)
                                                                                        cGAACCAGACGCTACATTCGCGGATGTACTCAATGTCGAACTCGAAGGGATTTTACATTTTCTCCAGCGTGTCGGGTTCGCCGCGCCATCCGTTAATCTGCGCCGATTCGGTAATGCCCTGTTTCACCTTATGGCGCAGCCTGttgagtt                                                              <  2:1907557/150‑7 (MQ=255)
                                                                                         gAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCCACCCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTc                                                             >  1:1797594/1‑150 (MQ=255)
                                                                                          aaCCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCa                                                            >  1:2677502/1‑150 (MQ=255)
                                                                                           aCCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCaa                                                           >  1:864052/1‑150 (MQ=255)
                                                                                            ccAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAAt                                                          <  1:549306/150‑1 (MQ=255)
                                                                                            ccAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAAt                                                          <  1:1399753/150‑1 (MQ=255)
                                                                                            ccAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAAt                                                          <  2:565234/150‑1 (MQ=255)
                                                                                                  gCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTg                                                    >  1:1959422/1‑150 (MQ=255)
                                                                                                  gCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTg                                                    >  1:1364551/1‑150 (MQ=255)
                                                                                                    tCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTc                                                  >  1:486945/1‑150 (MQ=255)
                                                                                                    tCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTATCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTc                                                  >  2:1970557/1‑150 (MQ=255)
                                                                                                     ccATTCGCGGATGTACTCACGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCg                                                 <  1:173707/150‑1 (MQ=255)
                                                                                                       aTTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGAt                                               <  2:4454488/150‑1 (MQ=255)
                                                                                                          cgcgGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACt                                            <  2:1466020/150‑1 (MQ=255)
                                                                                                           gcgGATGTACTCAAGGTCGAACTCGACGCGTTTGTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTg                                           >  2:3084622/1‑150 (MQ=255)
                                                                                                             ggATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGtt                                         <  1:373279/150‑1 (MQ=255)
                                                                                                             ggATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGtt                                         >  2:2567016/1‑150 (MQ=255)
                                                                                                             ggATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTATCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGtt                                         <  2:1504290/150‑1 (MQ=255)
                                                                                                             ggATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTATCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGtt                                         <  2:2522795/150‑1 (MQ=255)
                                                                                                              gATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTc                                        <  1:3881428/150‑1 (MQ=255)
                                                                                                               aTGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCg                                       >  2:1431590/1‑150 (MQ=255)
                                                                                                                tGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGt                                      <  1:1609661/150‑1 (MQ=255)
                                                                                                                 gTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCGCCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGtt                                     <  2:4062979/150‑1 (MQ=255)
                                                                                                                    cTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATg                                  <  2:1606716/150‑1 (MQ=255)
                                                                                                                        aGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCt                              >  2:1895671/1‑150 (MQ=255)
                                                                                                                           tCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACt                           <  1:1520129/150‑1 (MQ=255)
                                                                                                                            cGAACTCGACGCGTTTTTCCATTTTATCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACtg                          <  2:1180403/150‑1 (MQ=255)
                                                                                                                              aaCTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACtgcg                        <  1:1760716/150‑1 (MQ=255)
                                                                                                                                cTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACtgcgtg                      >  1:3313813/1‑150 (MQ=255)
                                                                                                                                     cgcgTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGAc                 <  2:61221/150‑1 (MQ=255)
                                                                                                                                     cgcgTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGAc                 >  2:1733237/1‑150 (MQ=255)
                                                                                                                                       cgTTTTTCCATTTTCTCCAGCGTGTCGGATTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGt               >  1:1586510/1‑150 (MQ=255)
                                                                                                                                         tttttCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTgg             >  2:817051/1‑150 (MQ=255)
                                                                                                                                          ttttCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGa            >  2:4025506/1‑150 (MQ=255)
                                                                                                                                              ccATTTTCTCCCGCGTGTCGGTTTCGCCGAGCCAGCCGTTTATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGa        <  1:923022/150‑1 (MQ=255)
                                                                                                                                              ccATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGa        >  1:2190190/1‑150 (MQ=255)
                                                                                                                                                  tttCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATTGCGCAGCATGTAGCCTTCAATGAGCTGTCGATAATGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAAt    >  1:412535/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCTGATGGTTTTTTCACGTAAGCTCATATCAATATGCCGCTTTGTTAACGAAACCTTTGAACACCGTCAGGAAAACGATTTTGATATCGAACCAGACGCTCCATTCGCGGATGTACTCAAGGTCGAACTCGACGCGTTTTTCCATTTTCTCCAGCGTGTCGGTTTCGCCGCGCCAGCCGTTAATCTGCGCCCAGCCGGTAATGCCCGGTTTCACCTTATGGCGCAGCATGTAGCCTTCAATGAGCTGTCGATACTGTTCGTTATGCGCTACTGCGTGCGGACGTGGACCGACAATCG  >  CP009273/2113612‑2113909

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: