Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,323,622 G→T 7.0% pseudogene (4260/4605 nt) yfaS ← pseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,323,6220GT 7.0% 744.0 / 20.5 199pseudogene (4260/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (80/105);  new base T (5/9);  total (85/114)
Fisher's exact test for biased strand distribution p-value = 7.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 199 total aligned reads displayed.

TTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATTTTTGC  >  CP009273/2323477‑2323768
                                                                                                                                                 |                                                                                                                                                  
tttCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTAt                                                                                                                                                  >  1:413367/1‑148 (MQ=255)
    cAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCa                                                                                                                                              <  1:3278775/148‑1 (MQ=255)
         aTTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgag                                                                                                                                         >  2:850720/1‑148 (MQ=255)
          ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCCTCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                        >  2:2238022/1‑148 (MQ=255)
          ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                        >  2:2536357/1‑148 (MQ=255)
          ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                        >  2:936138/1‑148 (MQ=255)
          ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                        <  2:562331/148‑1 (MQ=255)
           tGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGAt                                                                                                                                       <  2:388135/148‑1 (MQ=255)
            gCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATa                                                                                                                                      <  1:1493572/148‑1 (MQ=255)
            gCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATa                                                                                                                                      <  1:8718/148‑1 (MQ=255)
            gCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATa                                                                                                                                      <  1:1192023/148‑1 (MQ=255)
              cccTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACa                                                                                                                                    >  2:3367445/1‑148 (MQ=255)
                   ctgTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAgcgcc                                                                                                                               >  1:341594/1‑147 (MQ=255)
                   ctgTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCa                                                                                                                               <  2:3473700/148‑1 (MQ=255)
                      tttCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCg                                                                                                                            <  2:383620/148‑1 (MQ=255)
                      tttCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCg                                                                                                                            >  2:3551329/1‑148 (MQ=255)
                         cgcgGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTg                                                                                                                         <  2:1161058/148‑1 (MQ=255)
                             gctttGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCaa                                                                                                                     <  1:793167/145‑1 (MQ=255)
                             gcgTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCaa                                                                                                                     <  2:3326661/148‑1 (MQ=255)
                             gcgTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCaa                                                                                                                     >  2:48481/1‑148 (MQ=255)
                              cgTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAAt                                                                                                                    <  1:3091876/148‑1 (MQ=255)
                               gTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAAtc                                                                                                                   <  2:2473198/148‑1 (MQ=255)
                                ttGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAAtct                                                                                                                  >  1:28714/1‑148 (MQ=255)
                                ttGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCAACGAGATACAGCGCATCGCTGTCAAtct                                                                                                                  >  1:2689953/1‑148 (MQ=255)
                                  gggTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCa                                                                                                                <  1:782637/148‑1 (MQ=255)
                                   ggTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCAt                                                                                                               <  1:2339710/148‑1 (MQ=255)
                                   ggTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCAt                                                                                                               <  2:2685343/148‑1 (MQ=255)
                                   ggTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCAt                                                                                                               <  1:3068787/148‑1 (MQ=255)
                                       tATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCt                                                                                                           <  1:901792/148‑1 (MQ=255)
                                        aTTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTg                                                                                                          <  2:586527/148‑1 (MQ=255)
                                        aTTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTg                                                                                                          >  1:2976168/1‑148 (MQ=255)
                                            aCCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCa                                                                                                      >  1:2454325/1‑148 (MQ=255)
                                             ccGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCAc                                                                                                     >  2:3421090/1‑148 (MQ=255)
                                             ccGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCAc                                                                                                     >  2:401293/1‑148 (MQ=255)
                                              cGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCAcc                                                                                                    <  2:707233/148‑1 (MQ=255)
                                                 aaTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGt                                                                                                 >  1:1050158/1‑148 (MQ=255)
                                                 aaTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGt                                                                                                 >  1:2825523/1‑148 (MQ=255)
                                                 aaTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGt                                                                                                 >  1:123896/1‑148 (MQ=255)
                                                  aTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGtt                                                                                                >  2:1726291/1‑148 (MQ=255)
                                                  aTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGtt                                                                                                >  1:1106011/1‑148 (MQ=255)
                                                    gCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGc                                                                                              <  2:1408550/148‑1 (MQ=255)
                                                      cccATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAg                                                                                            >  2:2193212/1‑148 (MQ=255)
                                                       ccATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAACGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGa                                                                                           <  1:2225903/148‑1 (MQ=255)
                                                        cATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGaa                                                                                          <  2:2379177/148‑1 (MQ=255)
                                                        cATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGaa                                                                                          >  2:2912209/1‑148 (MQ=255)
                                                         aTGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGact                                                                                         >  2:2367364/1‑146 (MQ=255)
                                                           gttgtGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATaa                                                                                       <  2:1200124/148‑1 (MQ=255)
                                                            ttgtGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATaaa                                                                                      <  1:582375/148‑1 (MQ=255)
                                                             tgtgCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAg                                                                                     <  1:3125816/148‑1 (MQ=255)
                                                             tgtgCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAg                                                                                     <  2:1906931/148‑1 (MQ=255)
                                                              gtgCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGc                                                                                    <  1:2682578/148‑1 (MQ=255)
                                                                gcgcTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTc                                                                                  <  2:986660/148‑1 (MQ=255)
                                                                 cgcTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCa                                                                                 <  1:2499846/148‑1 (MQ=255)
                                                                  gcTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAt                                                                                <  1:2365497/148‑1 (MQ=255)
                                                                  gcTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAt                                                                                <  1:3241521/148‑1 (MQ=255)
                                                                    tCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtct                                                                              >  2:2854316/1‑148 (MQ=255)
                                                                    tCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtct                                                                              >  2:266315/1‑148 (MQ=255)
                                                                    tCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtct                                                                              <  2:2341489/148‑1 (MQ=255)
                                                                      aaCGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATctct                                                                            >  1:2327670/1‑148 (MQ=255)
                                                                           cGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCttcttc                                                                       <  1:2292804/148‑1 (MQ=255)
                                                                              cTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCAcc                                                                    <  2:3507502/148‑1 (MQ=255)
                                                                                 ccGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAgg                                                                 >  2:1685649/1‑148 (MQ=255)
                                                                                  cGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAggg                                                                <  2:342684/148‑1 (MQ=255)
                                                                                   tgTGGCAGCGTTACTTCCTCCTTACCGTAGCGCAGAACTGGATCCTGCTCGCTGGTAAGCGTGATTTAATAGATATACAGCGCATCGCTGTAAATCTCATTGCTGGTCACCGGTTTCATAATTAAGCTGATCTCTTATTCACCAGGGa                                                               <  1:3297426/147‑1 (MQ=255)
                                                                                       gCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAg                                                           >  2:1224366/1‑148 (MQ=255)
                                                                                       gCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAg                                                           >  2:1841193/1‑148 (MQ=255)
                                                                                       gCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAg                                                           >  2:244147/1‑148 (MQ=255)
                                                                                        cAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGc                                                          <  1:1835832/148‑1 (MQ=255)
                                                                                              tACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATAc                                                    >  1:1600039/1‑148 (MQ=255)
                                                                                              tACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATAc                                                    >  2:1285931/1‑148 (MQ=255)
                                                                                                 ttCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAAc                                                 >  1:2198871/1‑148 (MQ=255)
                                                                                                  tCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACt                                                >  2:1976159/1‑148 (MQ=255)
                                                                                                    cACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTgg                                              <  1:1298363/148‑1 (MQ=255)
                                                                                                        tGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGtt                                          >  2:1930557/1‑148 (MQ=255)
                                                                                                        tGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGtt                                          <  1:1564403/148‑1 (MQ=255)
                                                                                                          aCCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAAGCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCa                                        >  1:1342890/1‑148 (MQ=255)
                                                                                                            cGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAAc                                      >  1:1996487/1‑148 (MQ=255)
                                                                                                                 cgcAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCa                                 <  1:401285/148‑1 (MQ=255)
                                                                                                                     gAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCgg                             >  2:3084682/1‑148 (MQ=255)
                                                                                                                      aaCTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGa                            >  2:148224/1‑148 (MQ=255)
                                                                                                                      aaCTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGa                            <  1:1036030/148‑1 (MQ=255)
                                                                                                                       aCTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGaa                           <  2:806660/148‑1 (MQ=255)
                                                                                                                       aCTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGaa                           >  1:1300582/1‑148 (MQ=255)
                                                                                                                        cTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAAt                          >  1:3350379/1‑148 (MQ=255)
                                                                                                                             tCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTNCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTAc                     >  1:1805328/1‑148 (MQ=255)
                                                                                                                             tCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTAc                     >  2:3544865/1‑148 (MQ=255)
                                                                                                                             tCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTAc                     <  1:3367638/148‑1 (MQ=255)
                                                                                                                              ccTGCTCGCTGTTAAGCGTGATTTCATCGAGATACAGCGAATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACt                    <  1:850720/148‑1 (MQ=255)
                                                                                                                               cTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACtt                   <  2:3127326/148‑1 (MQ=255)
                                                                                                                               cTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACtt                   <  2:2962858/148‑1 (MQ=255)
                                                                                                                               cTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACtt                   >  2:2355343/1‑148 (MQ=255)
                                                                                                                                   tCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTttgtt               <  2:1124441/148‑1 (MQ=255)
                                                                                                                                        ggTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTATTTTGTTGAGcc          <  1:540916/148‑1 (MQ=255)
                                                                                                                                        ggTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGcc          >  1:2527161/1‑148 (MQ=255)
                                                                                                                                        ggTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGcc          <  1:266315/148‑1 (MQ=255)
                                                                                                                                        ggTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGcc          <  1:2854316/148‑1 (MQ=255)
                                                                                                                                        ggTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGcc          >  2:1160838/1‑148 (MQ=255)
                                                                                                                                         gTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCa         >  1:1855394/1‑148 (MQ=255)
                                                                                                                                             gCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAtttt     <  2:28714/148‑1 (MQ=255)
                                                                                                                                                tGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATTTTTGc  >  1:1639107/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATTTTTGC  >  CP009273/2323477‑2323768

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: