Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,495,918 C→A 6.3% V683F (GTT→TTT)  fryA ← putative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,495,9180CA 6.3% 650.2 / 13.7 174V683F (GTT→TTT) fryAputative PTS enzyme, Hpr component/enzyme I component/enzyme IIA component
Reads supporting (aligned to +/- strand):  ref base C (86/77);  new base A (4/7);  total (90/84)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 179 total aligned reads displayed.

CCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCAT  >  CP009273/2495771‑2496065
                                                                                                                                                   |                                                                                                                                                   
cccGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAc                                                                                                                                                     >  1:2718945/1‑148 (MQ=255)
   gCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGtc                                                                                                                                                  <  2:1670189/148‑1 (MQ=255)
     tGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtctt                                                                                                                                                <  2:1469081/148‑1 (MQ=255)
      gCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtcttc                                                                                                                                               >  2:805555/1‑148 (MQ=255)
       ccAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGtcttct                                                                                                                                              <  2:2268213/148‑1 (MQ=255)
        cAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTcttctt                                                                                                                                             <  2:675089/148‑1 (MQ=255)
        cAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTcttctt                                                                                                                                             <  1:792187/148‑1 (MQ=255)
            cATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCgg                                                                                                                                         >  2:1229671/1‑148 (MQ=255)
              tcttctTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCg                                                                                                                                       >  1:1825702/1‑148 (MQ=255)
               cttcttCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGt                                                                                                                                      >  2:2836390/1‑148 (MQ=255)
               cttcttCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGt                                                                                                                                      <  1:1912102/148‑1 (MQ=255)
               cttcttCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGt                                                                                                                                      >  2:1184026/1‑148 (MQ=255)
                ttcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGGTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTa                                                                                                                                     <  1:596063/148‑1 (MQ=255)
                 tcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaa                                                                                                                                    <  2:2738850/148‑1 (MQ=255)
                 tcttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaa                                                                                                                                    <  1:954176/148‑1 (MQ=255)
                  cttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                                                                                   <  1:2942294/148‑1 (MQ=255)
                  cttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                                                                                   >  2:3286020/1‑148 (MQ=255)
                  cttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                                                                                   <  2:1328122/148‑1 (MQ=255)
                  cttcCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTaaa                                                                                                                                   <  2:567603/148‑1 (MQ=255)
                         cTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                            <  2:485337/148‑1 (MQ=255)
                         cTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTg                                                                                                                            >  2:2273526/1‑148 (MQ=255)
                           cAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAg                                                                                                                          >  2:383981/1‑148 (MQ=255)
                            aaaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGGACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGt                                                                                                                         >  1:1986440/1‑148 (MQ=255)
                             aaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTa                                                                                                                        <  1:1178042/148‑1 (MQ=255)
                             aaCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTa                                                                                                                        >  1:1012385/1‑148 (MQ=255)
                              aCGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTaa                                                                                                                       <  1:3062986/148‑1 (MQ=255)
                               cGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAAc                                                                                                                      >  2:1225116/1‑148 (MQ=255)
                                 gATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGc                                                                                                                    <  2:1358400/148‑1 (MQ=255)
                                  aTGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCt                                                                                                                   >  2:2129469/1‑148 (MQ=255)
                                   tGTTCGGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCtt                                                                                                                  <  1:1610876/148‑1 (MQ=255)
                                   tGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCtt                                                                                                                  <  1:1610880/148‑1 (MQ=255)
                                   tGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCtt                                                                                                                  <  2:1675145/148‑1 (MQ=255)
                                    gTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTc                                                                                                                 >  2:2672565/1‑148 (MQ=255)
                                      tCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCaa                                                                                                               <  1:3029392/148‑1 (MQ=255)
                                          tGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTc                                                                                                           <  2:64018/148‑1 (MQ=255)
                                          tGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTc                                                                                                           <  2:314633/148‑1 (MQ=255)
                                            cgcCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCt                                                                                                         >  1:841240/1‑148 (MQ=255)
                                             gcCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTg                                                                                                        <  2:1909234/148‑1 (MQ=255)
                                              cccGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTgc                                                                                                       >  1:2619267/1‑148 (MQ=255)
                                               ccGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTgcg                                                                                                      >  1:930635/1‑148 (MQ=255)
                                                cGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTgcgc                                                                                                     >  1:3204753/1‑148 (MQ=255)
                                                 gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCa                                                                                                    <  1:2848306/148‑1 (MQ=255)
                                                 gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCa                                                                                                    <  2:1612236/148‑1 (MQ=255)
                                                 gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCa                                                                                                    <  1:2848248/148‑1 (MQ=255)
                                                 gTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGATCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTCACGCTGCAATTTCCTGCGCa                                                                                                    >  2:2063062/1‑148 (MQ=255)
                                                           aGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCAtt                                                                                          <  1:937365/148‑1 (MQ=255)
                                                              ttGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcac                                                                                       <  1:589177/148‑1 (MQ=255)
                                                              ttGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcac                                                                                       <  1:3101732/148‑1 (MQ=255)
                                                               tGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcaca                                                                                      <  1:330485/148‑1 (MQ=255)
                                                               tGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTcaca                                                                                      <  1:330367/148‑1 (MQ=255)
                                                                gCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACAt                                                                                     <  1:3022252/148‑1 (MQ=255)
                                                                  cGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGc                                                                                   <  1:207007/148‑1 (MQ=255)
                                                                   gCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGcc                                                                                  >  2:381266/1‑148 (MQ=255)
                                                                     acaGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTg                                                                                <  1:1075990/148‑1 (MQ=255)
                                                                      caGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGa                                                                               >  1:459218/1‑148 (MQ=255)
                                                                        ggAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                                                                             >  2:3003037/1‑148 (MQ=255)
                                                                        ggAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcg                                                                             >  2:3002704/1‑148 (MQ=255)
                                                                        ggAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCAGTCACATGCCTGAcg                                                                             >  1:3017298/1‑148 (MQ=255)
                                                                         gAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGAcgc                                                                            >  2:2624083/1‑148 (MQ=255)
                                                                           aCTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACgcgc                                                                          <  1:1170039/148‑1 (MQ=255)
                                                                                    ccTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCg                                                                 >  2:1457642/1‑148 (MQ=255)
                                                                                     cTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGa                                                                >  2:1077821/1‑148 (MQ=255)
                                                                                       gCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCCGTTCCCGaca                                                              >  2:2346224/1‑148 (MQ=255)
                                                                                       gCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGaca                                                              <  1:3136545/148‑1 (MQ=255)
                                                                                           tttATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAAAGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGcc                                                          <  1:1225116/148‑1 (MQ=255)
                                                                                            ttATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCt                                                         >  2:790463/1‑148 (MQ=255)
                                                                                            ttATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCt                                                         <  2:1273281/148‑1 (MQ=255)
                                                                                                  cTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTAt                                                   <  2:2498778/148‑1 (MQ=255)
                                                                                                    aaaaTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATcc                                                 >  1:1547798/1‑148 (MQ=255)
                                                                                                       aTCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGt                                              >  2:574144/1‑148 (MQ=255)
                                                                                                        tCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGtt                                             >  1:1405024/1‑148 (MQ=255)
                                                                                                         ccTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACATGCCTCGCTATCCAGTTg                                            >  1:1322871/1‑148 (MQ=255)
                                                                                                         ccTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTg                                            >  1:2970321/1‑148 (MQ=255)
                                                                                                         ccTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTg                                            >  2:23924/1‑148 (MQ=255)
                                                                                                          cTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGa                                           <  2:3219483/148‑1 (MQ=255)
                                                                                                            gATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACg                                         >  1:19112/1‑148 (MQ=255)
                                                                                                              tCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGaa                                       >  2:1850023/1‑148 (MQ=255)
                                                                                                              tCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGaa                                       >  2:2017186/1‑148 (MQ=255)
                                                                                                                     aGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCt                                >  2:366083/1‑148 (MQ=255)
                                                                                                                      gATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCtt                               <  2:2507487/148‑1 (MQ=255)
                                                                                                                       atatTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCttt                              >  1:2538380/1‑148 (MQ=255)
                                                                                                                       atatTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCttt                              <  1:1247560/148‑1 (MQ=255)
                                                                                                                       atatTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCttt                              >  1:537228/1‑148 (MQ=255)
                                                                                                                           tctcCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAc                          >  2:890120/1‑148 (MQ=255)
                                                                                                                                aGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCAccgccg                     >  1:2393635/1‑148 (MQ=255)
                                                                                                                                  cgcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGa                   >  2:3123389/1‑148 (MQ=255)
                                                                                                                                  cgcCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGa                   >  2:2278577/1‑148 (MQ=255)
                                                                                                                                    ccAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTTTCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAAt                 <  1:2192496/148‑1 (MQ=255)
                                                                                                                                      agcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAc               >  2:402908/1‑148 (MQ=255)
                                                                                                                                       gcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcg              <  2:3285796/148‑1 (MQ=255)
                                                                                                                                       gcagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcg              <  2:1831627/148‑1 (MQ=255)
                                                                                                                                        cagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgc             >  2:1389680/1‑148 (MQ=255)
                                                                                                                                        cagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgc             >  2:2447280/1‑148 (MQ=255)
                                                                                                                                        cagTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATAcgc             >  2:2447224/1‑148 (MQ=255)
                                                                                                                                           tGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCggg          >  2:427025/1‑148 (MQ=255)
                                                                                                                                            gggCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGc         >  2:1721011/1‑148 (MQ=255)
                                                                                                                                                 aaCGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactc    <  1:2488836/148‑1 (MQ=255)
                                                                                                                                                   cGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcat  <  1:2273526/148‑1 (MQ=255)
                                                                                                                                                   cGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCtactcat  >  2:1521468/1‑148 (MQ=255)
                                                                                                                                                   cGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTGTTCACCGCCGGAATACGCGGGCtactcat  >  2:2953359/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCGCTGCCAGACATCTTCTTCCAGCTCAAACGGATGTTCAGTGCGCCCGTTAACGCCGAGGTTGCCGCACAGGAACTGGATCGCCTGCTCTTTATTGCTAAAATCCTGATCAACAAAGATATTCTCCAGCGCCAGCAGTGGGCGAACGTCTTCTTCCGGCGTAAAGGCGGTGAGTAACGCTTCAATTTCCTGCGCACTGCGGCATTCACATGCCTGACGCGCCAGTTCCCGACACGCCTCGCTATCCAGTTGACGAAGCTGGCTTTTCACCGCCGGAATACGCGGGCTACTCAT  >  CP009273/2495771‑2496065

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: