Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,542,878 A→C 5.2% intergenic (+113/+127) yfeW → / ← yfeX penicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,542,8780AC 5.2% 558.2 / 18.6 192intergenic (+113/+127)yfeW/yfeXpenicillin binding protein PBP4B; weak DD‑carboxypeptidase activity/porphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base A (91/90);  new base C (9/2);  total (100/92)
Fisher's exact test for biased strand distribution p-value = 6.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

TTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGT  >  CP009273/2542731‑2543024
                                                                                                                                                   |                                                                                                                                                  
ttGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTa                                                                                                                                                    <  1:1976194/148‑1 (MQ=255)
ttGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTa                                                                                                                                                    >  1:172707/1‑148 (MQ=255)
 tGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAg                                                                                                                                                   <  2:727192/148‑1 (MQ=255)
 tGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAg                                                                                                                                                   >  1:3018390/1‑148 (MQ=255)
      aGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGa                                                                                                                                              <  1:2557976/148‑1 (MQ=255)
      aGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGa                                                                                                                                              <  1:3123886/148‑1 (MQ=255)
       ggTGTATGCGGCGTTAAAGCAGAAGTAATACAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGAt                                                                                                                                             <  2:1064058/148‑1 (MQ=255)
         tgtATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGAttt                                                                                                                                           <  2:3036168/148‑1 (MQ=255)
         tgtATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGGTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGAttt                                                                                                                                           <  2:2343294/148‑1 (MQ=255)
                   gTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCTCGATTTCGTAGGCCgg                                                                                                                                 >  2:786700/1‑148 (MQ=255)
                   gTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCgg                                                                                                                                 <  2:2611538/148‑1 (MQ=255)
                     tAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGAt                                                                                                                               <  1:724929/148‑1 (MQ=255)
                       aaGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTCGCACGATTTCGTAGGCCGGATaa                                                                                                                             >  1:2767671/1‑148 (MQ=255)
                       aaGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATaa                                                                                                                             >  1:2824092/1‑148 (MQ=255)
                        aGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAg                                                                                                                            <  1:1093208/148‑1 (MQ=255)
                        aGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAg                                                                                                                            >  1:903565/1‑148 (MQ=255)
                        aGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAg                                                                                                                            >  2:1329794/1‑148 (MQ=255)
                          cAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGc                                                                                                                          >  2:697549/1‑148 (MQ=255)
                           aGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCg                                                                                                                         <  2:3195149/148‑1 (MQ=255)
                           aGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCg                                                                                                                         >  1:1895184/1‑148 (MQ=255)
                            gAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGt                                                                                                                        <  2:763528/148‑1 (MQ=255)
                             aaGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                                                                       >  1:1623260/1‑148 (MQ=255)
                             aaGTAATAAAAAGGCCAGTCAGTAATGTTACTGGATGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                                                                       >  2:1403851/1‑148 (MQ=255)
                                taataaAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATAGTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAc                                                                                                                    >  2:1830461/1‑148 (MQ=255)
                                   taaAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                 >  2:1682261/1‑148 (MQ=255)
                                      aaaGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                                              >  1:1914731/1‑148 (MQ=255)
                                       aaGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                             >  2:2927396/1‑148 (MQ=255)
                                       aaGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                             <  1:2036651/148‑1 (MQ=255)
                                        aGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                            >  1:1203643/1‑148 (MQ=255)
                                         ggCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                           >  1:2739911/1‑148 (MQ=255)
                                           ccAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                         >  2:1003251/1‑148 (MQ=255)
                                           ccAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATAGTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                         <  2:2210123/148‑1 (MQ=255)
                                               tcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtg                                                                                                     <  2:2955937/148‑1 (MQ=255)
                                               tcagtAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtg                                                                                                     >  2:1930893/1‑148 (MQ=255)
                                                    aaTGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACACTATATTGACTTAGCCCGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACaa                                                                                                <  2:2783834/148‑1 (MQ=255)
                                                     aTGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAAc                                                                                               >  2:548139/1‑148 (MQ=255)
                                                      tGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAAcg                                                                                              >  1:156254/1‑148 (MQ=255)
                                                         tACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACgcgc                                                                                           <  2:2265057/148‑1 (MQ=255)
                                                           ctggctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCAc                                                                                         <  2:824174/148‑1 (MQ=255)
                                                             ggctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACtt                                                                                       >  2:2160053/1‑148 (MQ=255)
                                                              gctggTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTg                                                                                      <  1:2095066/148‑1 (MQ=255)
                                                                  gTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCACCCGGCGTGTACAACGCGCACTTGGTca                                                                                  <  2:2015122/148‑1 (MQ=255)
                                                                    cTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTcaac                                                                                <  2:2855720/148‑1 (MQ=255)
                                                                     tGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTcaaca                                                                               <  2:2877084/148‑1 (MQ=255)
                                                                      gagaGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTcaacaa                                                                              <  2:2369054/148‑1 (MQ=255)
                                                                       agagTGCTGTCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAAt                                                                             >  2:1662613/1‑148 (MQ=255)
                                                                          gTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTg                                                                          <  2:976517/148‑1 (MQ=255)
                                                                           tGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTCGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                         >  2:126445/1‑148 (MQ=255)
                                                                           tGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                         >  1:1078985/1‑148 (MQ=255)
                                                                           tGCTGACAAACGCAAAACTGCCTGATGCGCTACGATTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGa                                                                         <  1:1506326/148‑1 (MQ=255)
                                                                            gCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAg                                                                        >  1:2078033/1‑148 (MQ=255)
                                                                               gACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGcc                                                                     >  1:2929840/1‑148 (MQ=255)
                                                                                aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTCGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCca                                                                    >  2:1489443/1‑148 (MQ=255)
                                                                                aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCca                                                                    >  2:957360/1‑148 (MQ=255)
                                                                                aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCca                                                                    >  2:3070365/1‑148 (MQ=255)
                                                                                aCAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCca                                                                    >  1:145892/1‑148 (MQ=255)
                                                                                    aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagtca                                                                >  2:2687111/1‑148 (MQ=255)
                                                                                    aCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagtca                                                                >  2:2525989/1‑148 (MQ=255)
                                                                                          aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATg                                                          >  1:697734/1‑148 (MQ=255)
                                                                                          aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATg                                                          >  2:2756867/1‑148 (MQ=255)
                                                                                          aaCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACGTGGTCAACAATCGGAGGCCAGTCCGTAATg                                                          >  1:2359080/1‑148 (MQ=255)
                                                                                            cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGtt                                                        <  1:369788/148‑1 (MQ=255)
                                                                                             tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTa                                                       <  2:1030454/148‑1 (MQ=255)
                                                                                             tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTa                                                       >  1:3132375/1‑148 (MQ=255)
                                                                                             tGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTa                                                       <  1:2647555/148‑1 (MQ=255)
                                                                                                 tGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGc                                                   <  2:2319108/148‑1 (MQ=255)
                                                                                                 tGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGc                                                   >  2:2860585/1‑148 (MQ=255)
                                                                                                   aTGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCt                                                 >  1:1658662/1‑148 (MQ=255)
                                                                                                      cgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                              <  2:2739911/148‑1 (MQ=255)
                                                                                                      cgcTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGc                                              <  2:1994410/148‑1 (MQ=255)
                                                                                                       gcTACGCTTATCAGGCCAACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGcc                                             <  2:252738/148‑1 (MQ=255)
                                                                                                         tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCtt                                           <  2:1174715/148‑1 (MQ=255)
                                                                                                         tACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCtt                                           <  1:1102042/148‑1 (MQ=255)
                                                                                                             cTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAAt                                       <  1:1552303/148‑1 (MQ=255)
                                                                                                               tATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCa                                     >  2:1823220/1‑148 (MQ=255)
                                                                                                                  cAGGCCTACGTGGAGCCCACAATATATTGAAGTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATg                                  <  2:2977348/148‑1 (MQ=255)
                                                                                                                    ggCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGaa                                <  2:341307/148‑1 (MQ=255)
                                                                                                                         aCGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGa                           <  2:2641848/148‑1 (MQ=255)
                                                                                                                          cGTGGAGCCCACAATATATTGAATTCGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGaa                          <  1:240557/148‑1 (MQ=255)
                                                                                                                           gTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGaaa                         <  1:1134400/148‑1 (MQ=255)
                                                                                                                           gTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGaaa                         <  2:3068639/148‑1 (MQ=255)
                                                                                                                             ggAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACg                       <  2:118597/148‑1 (MQ=255)
                                                                                                                                gCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCtt                    <  2:903849/148‑1 (MQ=255)
                                                                                                                                 cccACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATTAATCAGAAACGCTTa                   <  1:1432299/148‑1 (MQ=255)
                                                                                                                                  ccACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTAc                  >  2:3122791/1‑148 (MQ=255)
                                                                                                                                   cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACa                 >  2:312164/1‑148 (MQ=255)
                                                                                                                                   cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACa                 >  2:2508511/1‑148 (MQ=255)
                                                                                                                                   cacaATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACa                 >  1:3098244/1‑148 (MQ=255)
                                                                                                                                     caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgc               >  2:2345656/1‑148 (MQ=255)
                                                                                                                                       atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAgcgc             >  1:3126911/1‑148 (MQ=255)
                                                                                                                                           atTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATc         >  2:2552002/1‑148 (MQ=255)
                                                                                                                                             tGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGACTCAGAAACGCTTACAGCGCCATCaa       >  1:139004/1‑148 (MQ=255)
                                                                                                                                             tGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCaa       >  2:2641231/1‑148 (MQ=255)
                                                                                                                                             tGAATGAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCaa       <  2:2084549/148‑1 (MQ=255)
                                                                                                                                                 ttAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTg   >  2:606150/1‑148 (MQ=255)
                                                                                                                                                 ttAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTg   >  2:806652/1‑148 (MQ=255)
                                                                                                                                                 ttAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAAAAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTg   <  2:2824092/148‑1 (MQ=255)
                                                                                                                                                  tAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGt  >  1:1545183/1‑148 (MQ=255)
                                                                                                                                                  tAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGt  >  1:2820735/1‑148 (MQ=255)
                                                                                                                                                  tAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGt  >  2:292087/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGATCAGGTGTATGCGGCGTTAAAGCAGAAGTAATAAAAAGGCCAGTCAGTAATGTTACTGGCTGGTCTGAGAGTGCTGACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGCCTGGCCTTTAATCAATGAATCAGAAACGCTTACAGCGCCATCAACTTGT  >  CP009273/2542731‑2543024

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: