Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,397,366 A→T 13.9% S49S (TCT→TCA)  abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,397,3660AT13.9% 615.3 / 78.0 215S49S (TCT→TCA) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base A (77/108);  new base T (12/18);  total (89/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 216 total aligned reads displayed.

GTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATC  >  CP009273/1397220‑1397513
                                                                                                                                                  |                                                                                                                                                   
gTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAg                                                                                                                                                    <  1:37718/148‑1 (MQ=255)
   ggTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGAtt                                                                                                                                                 <  1:1142608/148‑1 (MQ=255)
   ggTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGAtt                                                                                                                                                 <  2:3096918/148‑1 (MQ=255)
     ttGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTcc                                                                                                                                               <  2:3373506/148‑1 (MQ=255)
            cAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCg                                                                                                                                        >  1:2518040/1‑148 (MQ=255)
               ccTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAg                                                                                                                                     >  1:3535432/1‑148 (MQ=255)
                cTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGc                                                                                                                                    <  2:446224/148‑1 (MQ=255)
                cTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGc                                                                                                                                    <  2:446323/148‑1 (MQ=255)
                cTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGc                                                                                                                                    <  2:3271449/148‑1 (MQ=255)
                  gcttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCa                                                                                                                                  >  1:1478031/1‑148 (MQ=255)
                  gcttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCa                                                                                                                                  >  1:1628581/1‑148 (MQ=255)
                   cttgctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAg                                                                                                                                 <  2:665553/148‑1 (MQ=255)
                      gctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATg                                                                                                                              >  1:3333151/1‑148 (MQ=255)
                      gctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATg                                                                                                                              <  1:1753088/148‑1 (MQ=255)
                       ctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGc                                                                                                                             <  2:2768144/148‑1 (MQ=255)
                       ctGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGc                                                                                                                             <  1:1549566/148‑1 (MQ=255)
                           cTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCg                                                                                                                         <  2:2907075/148‑1 (MQ=255)
                             tAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCt                                                                                                                       >  2:1864479/1‑148 (MQ=255)
                               aaCCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGa                                                                                                                     >  1:2166893/1‑148 (MQ=255)
                                aCCTGCCAGGGCGTCATAGTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGAc                                                                                                                    <  1:2922873/148‑1 (MQ=255)
                                   tGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAg                                                                                                                 <  2:3385905/148‑1 (MQ=255)
                                    gCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGa                                                                                                                <  2:3215639/148‑1 (MQ=255)
                                     ccAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGaa                                                                                                               >  1:3298770/1‑148 (MQ=255)
                                       aGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAAct                                                                                                             >  1:1173160/1‑148 (MQ=255)
                                              cATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCaaa                                                                                                      <  2:807137/148‑1 (MQ=255)
                                              cATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCaaa                                                                                                      >  2:2988812/1‑148 (MQ=255)
                                                   tctcCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACgtg                                                                                                 <  1:182704/148‑1 (MQ=255)
                                                    ctcCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACgtgt                                                                                                <  1:1032577/148‑1 (MQ=255)
                                                        cagcagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCt                                                                                            <  1:1502350/148‑1 (MQ=255)
                                                           cagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                         <  1:2812508/148‑1 (MQ=255)
                                                           cagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                         >  1:1916570/1‑148 (MQ=255)
                                                           cagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                         <  1:675519/148‑1 (MQ=255)
                                                           cagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                         >  1:2685776/1‑148 (MQ=255)
                                                           cagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                         <  2:2930250/148‑1 (MQ=255)
                                                           cagGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGga                                                                                         >  1:2541487/1‑148 (MQ=255)
                                                             gggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatg                                                                                       <  1:2309304/148‑1 (MQ=255)
                                                              ggCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGgatga                                                                                      <  2:797374/148‑1 (MQ=255)
                                                                cGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATc                                                                                    <  2:1749454/148‑1 (MQ=255)
                                                                cGATAACCGGTTGGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATc                                                                                    <  2:1620969/148‑1 (MQ=255)
                                                                  aTAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATccc                                                                                  >  1:1484103/1‑148 (MQ=255)
                                                                     aCCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGa                                                                               <  2:2046406/148‑1 (MQ=255)
                                                                      ccGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGAt                                                                              <  1:2512746/148‑1 (MQ=255)
                                                                      ccGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGAt                                                                              <  2:88584/148‑1 (MQ=255)
                                                                           ttGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGAt                                                                         >  1:945453/1‑148 (MQ=255)
                                                                           ttGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTACAAAACGTGTTTCTGGATGATCCCAGATTTGAt                                                                         >  2:1088751/1‑148 (MQ=255)
                                                                             gCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCg                                                                       <  2:1268117/148‑1 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    <  1:54027/148‑1 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    >  2:1612235/1‑148 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    >  2:1611866/1‑148 (MQ=255)
                                                                                ttGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCa                                                                    <  2:1505460/148‑1 (MQ=255)
                                                                                 tGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCaa                                                                   <  2:3495320/148‑1 (MQ=255)
                                                                                     cAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATc                                                               <  2:54273/148‑1 (MQ=255)
                                                                                       aaCGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCgg                                                             >  2:3344709/1‑148 (MQ=255)
                                                                                       aaCGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCgg                                                             >  1:1005600/1‑148 (MQ=255)
                                                                                         cGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGta                                                           >  1:2927722/1‑148 (MQ=255)
                                                                                         cGAAGCAATAAAGGCATTTGGGATAGTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGta                                                           <  2:2155393/148‑1 (MQ=255)
                                                                                            aGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATaa                                                        <  2:1650263/148‑1 (MQ=255)
                                                                                              cAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACg                                                      <  1:922537/148‑1 (MQ=255)
                                                                                               aaTAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGt                                                     <  1:3059037/148‑1 (MQ=255)
                                                                                               aaTAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGt                                                     <  1:2482976/148‑1 (MQ=255)
                                                                                                 tAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTg                                                   <  2:520363/148‑1 (MQ=255)
                                                                                                  aaaGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTgg                                                  <  2:2504105/148‑1 (MQ=255)
                                                                                                  aaaGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTgg                                                  >  2:1946368/1‑148 (MQ=255)
                                                                                                    aGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCg                                                <  1:3428322/148‑1 (MQ=255)
                                                                                                     ggCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGa                                               >  1:1052328/1‑148 (MQ=255)
                                                                                                     ggCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGa                                               <  1:2892015/148‑1 (MQ=255)
                                                                                                     ggCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGa                                               <  2:2546115/148‑1 (MQ=255)
                                                                                                         tttGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCgg                                           <  1:689163/148‑1 (MQ=255)
                                                                                                         tttGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCgg                                           >  1:1151812/1‑148 (MQ=255)
                                                                                                           tGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCt                                         <  2:1787272/148‑1 (MQ=255)
                                                                                                               atatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCaa                                     <  2:1057300/148‑1 (MQ=255)
                                                                                                               atatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCaa                                     <  1:1408186/148‑1 (MQ=255)
                                                                                                               atatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCaa                                     <  2:2006074/148‑1 (MQ=255)
                                                                                                                tatTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAAt                                    <  2:928073/148‑1 (MQ=255)
                                                                                                                 atTGCCTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATc                                   >  1:113322/1‑148 (MQ=255)
                                                                                                                 atTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATc                                   >  1:2289796/1‑148 (MQ=255)
                                                                                                                  ttGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCg                                  >  1:3480863/1‑148 (MQ=255)
                                                                                                                   tGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGc                                 >  2:2184659/1‑148 (MQ=255)
                                                                                                                   tGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGc                                 >  1:2313774/1‑148 (MQ=255)
                                                                                                                   tGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGc                                 <  2:1478031/148‑1 (MQ=255)
                                                                                                                    gCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCa                                <  1:3156802/148‑1 (MQ=255)
                                                                                                                     ccTACGTTGCGGGTAACGGTGAAGCCTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAt                               >  2:2419661/1‑148 (MQ=255)
                                                                                                                       tACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcg                             <  1:1895899/148‑1 (MQ=255)
                                                                                                                        aCGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgt                            >  2:1143823/1‑148 (MQ=255)
                                                                                                                        aCGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgt                            >  2:1143869/1‑148 (MQ=255)
                                                                                                                        aCGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgt                            <  2:2685776/148‑1 (MQ=255)
                                                                                                                         cGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCAtcgtc                           >  2:137964/1‑148 (MQ=255)
                                                                                                                               gggTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaa                     <  1:1810594/148‑1 (MQ=255)
                                                                                                                                ggTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATaaaa                    <  2:258190/148‑1 (MQ=255)
                                                                                                                                 gTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAAc                   >  2:2065973/1‑148 (MQ=255)
                                                                                                                                    aCGGTGAAGCCTGCAGAGTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCTGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGAt                <  2:2395454/148‑1 (MQ=255)
                                                                                                                                     cGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATa               >  2:367569/1‑148 (MQ=255)
                                                                                                                                     cGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATa               >  2:2360834/1‑148 (MQ=255)
                                                                                                                                     cGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATa               <  2:678859/148‑1 (MQ=255)
                                                                                                                                            gCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTcc        <  1:1538709/148‑1 (MQ=255)
                                                                                                                                             ccTGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCt       >  1:135998/1‑148 (MQ=255)
                                                                                                                                               tGCTGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGc     >  1:1289299/1‑148 (MQ=255)
                                                                                                                                                 cAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCAt   >  1:3003330/1‑148 (MQ=255)
                                                                                                                                                 cAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCAt   >  1:3033217/1‑148 (MQ=255)
                                                                                                                                                  aGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATc  <  1:734287/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTAGGTTGCGCGCAACCTGCTTGCTGACTTAAACCTGCCAGGGCGTCATATTCTCCCAGCAGGGCGATAACCGGTTTGCCTTGACCAAACGAAGCAATAAAGGCATTTGGGATATTGCCTACGTTGCGGGTAACGGTGAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATC  >  CP009273/1397220‑1397513

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: ‑