Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,323,622 G→T 9.8% 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 9.8% 528.4 / 24.9 152pseudogene (4260/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (53/84);  new base T (6/9);  total (59/93)
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.77e-01

Only 100 of 156 total aligned reads displayed.

TTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATTTTTGCC  >  CP009273/2323476‑2323769
                                                                                                                                                  |                                                                                                                                                   
ttttCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGa                                                                                                                                                    <  1:2581230/148‑1 (MQ=255)
ttttCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGa                                                                                                                                                    >  2:2788130/1‑148 (MQ=255)
ttttCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGa                                                                                                                                                    <  1:1742337/148‑1 (MQ=255)
ttttCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTGGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGa                                                                                                                                                    >  2:2788278/1‑148 (MQ=255)
 tttCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGAt                                                                                                                                                   <  2:1720109/148‑1 (MQ=255)
     cAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCa                                                                                                                                               >  2:387090/1‑148 (MQ=255)
        cAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCg                                                                                                                                            >  2:1411531/1‑148 (MQ=255)
          aTTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgag                                                                                                                                          <  2:3522095/148‑1 (MQ=255)
          aTTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgag                                                                                                                                          <  2:1010965/148‑1 (MQ=255)
           ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGTCCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                         >  1:1313811/1‑148 (MQ=255)
           ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                         <  1:2734162/148‑1 (MQ=255)
           ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                         >  2:27376/1‑148 (MQ=255)
           ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                         <  1:850132/148‑1 (MQ=255)
           ttGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCgaga                                                                                                                                         <  1:850229/148‑1 (MQ=255)
            tGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGAt                                                                                                                                        <  2:1493816/148‑1 (MQ=255)
            tGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGAt                                                                                                                                        >  2:817745/1‑148 (MQ=255)
            tGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGAt                                                                                                                                        <  2:989515/148‑1 (MQ=255)
              ccccTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCAAGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATAc                                                                                                                                      >  1:637379/1‑148 (MQ=255)
                   gctgTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAgcgc                                                                                                                                 <  1:1015584/148‑1 (MQ=255)
                   gctgTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAgcgc                                                                                                                                 <  1:2854879/148‑1 (MQ=255)
                    ctgTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCa                                                                                                                                >  2:2378436/1‑148 (MQ=255)
                       tttCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCg                                                                                                                             <  2:918207/148‑1 (MQ=255)
                       tttCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCg                                                                                                                             <  2:918015/148‑1 (MQ=255)
                             ggcgTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCa                                                                                                                       <  1:2930309/148‑1 (MQ=255)
                             ggcgTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCa                                                                                                                       <  1:2930353/148‑1 (MQ=255)
                             ggcgTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCa                                                                                                                       <  2:767225/148‑1 (MQ=255)
                                gTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGATGTCAAtc                                                                                                                    <  1:27376/148‑1 (MQ=255)
                                   gggTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCa                                                                                                                 <  2:3498759/148‑1 (MQ=255)
                                    ggTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCAt                                                                                                                <  1:3591/148‑1 (MQ=255)
                                      tttATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTg                                                                                                              <  2:812659/148‑1 (MQ=255)
                                       ttATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGc                                                                                                             >  1:1905576/1‑148 (MQ=255)
                                          ttGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTgg                                                                                                          <  1:3063155/148‑1 (MQ=255)
                                           tGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGt                                                                                                         <  2:555059/148‑1 (MQ=255)
                                           tGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGt                                                                                                         <  2:555191/148‑1 (MQ=255)
                                            gACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTc                                                                                                        <  1:3436321/148‑1 (MQ=255)
                                              ccGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCAc                                                                                                      <  2:3241208/148‑1 (MQ=255)
                                                 aaaTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCgg                                                                                                   <  2:561235/148‑1 (MQ=255)
                                                 aaaTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCgg                                                                                                   <  1:1210841/148‑1 (MQ=255)
                                                 aaaTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCgg                                                                                                   <  1:1211073/148‑1 (MQ=255)
                                                   aTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCgg                                                                                                   <  2:1903517/146‑1 (MQ=255)
                                                   aTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCgg                                                                                                   >  1:1063036/1‑146 (MQ=255)
                                                   aTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCgg                                                                                                   <  2:1063036/146‑1 (MQ=255)
                                                   aTGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGtt                                                                                                 >  2:630925/1‑148 (MQ=255)
                                                    tGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTg                                                                                                >  2:3411369/1‑148 (MQ=255)
                                                        ccATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGa                                                                                            <  1:1584751/148‑1 (MQ=255)
                                                        ccATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGa                                                                                            <  2:3212362/148‑1 (MQ=255)
                                                        ccATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGa                                                                                            <  1:1348599/148‑1 (MQ=255)
                                                        ccATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGa                                                                                            <  2:275795/148‑1 (MQ=255)
                                                         cATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGaa                                                                                           <  2:3400337/148‑1 (MQ=255)
                                                         cATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGaa                                                                                           <  1:3336295/148‑1 (MQ=255)
                                                          aTGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAAt                                                                                          >  2:245617/1‑148 (MQ=255)
                                                           tgttgtGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATa                                                                                         >  2:2933723/1‑148 (MQ=255)
                                                              tgtgCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAg                                                                                      <  2:1461262/148‑1 (MQ=255)
                                                                    cTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtc                                                                                >  1:242509/1‑148 (MQ=255)
                                                                     tCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtct                                                                               <  1:2173004/148‑1 (MQ=255)
                                                                     tCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtct                                                                               <  1:1571675/148‑1 (MQ=255)
                                                                     tCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtct                                                                               <  1:316828/148‑1 (MQ=255)
                                                                      cAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtctc                                                                              >  2:3050848/1‑148 (MQ=255)
                                                                      cAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCAtctc                                                                              >  2:1274988/1‑148 (MQ=255)
                                                                         cGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCtcttc                                                                           >  1:2591615/1‑148 (MQ=255)
                                                                         cGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCtcttc                                                                           >  2:169435/1‑148 (MQ=255)
                                                                          gTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCtcttct                                                                          <  2:946325/148‑1 (MQ=255)
                                                                          gTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCtcttct                                                                          <  2:946346/148‑1 (MQ=255)
                                                                              gCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCAc                                                                      <  2:706092/148‑1 (MQ=255)
                                                                                    ggTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGa                                                                >  2:2597492/1‑148 (MQ=255)
                                                                                    ggTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGa                                                                <  2:1851237/148‑1 (MQ=255)
                                                                                      tGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATa                                                              <  1:1143095/148‑1 (MQ=255)
                                                                                      tGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATa                                                              <  1:809305/148‑1 (MQ=255)
                                                                                       ggCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATaa                                                             <  2:1271743/148‑1 (MQ=255)
                                                                                            cGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTc                                                                                                        >  2:3058318/1‑100 (MQ=255)
                                                                                                    ccACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTg                                                <  1:1665050/148‑1 (MQ=255)
                                                                                                        cTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGt                                            <  1:1411531/148‑1 (MQ=255)
                                                                                                               tAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACgg                                     >  1:742756/1‑148 (MQ=255)
                                                                                                                aGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGt                                    <  2:3553102/148‑1 (MQ=255)
                                                                                                                  cgcAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCa                                  <  2:2329761/148‑1 (MQ=255)
                                                                                                                        aCTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCTACGGTCACCGGaa                            <  1:2555645/148‑1 (MQ=255)
                                                                                                                           gCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGt                         >  1:1008815/1‑148 (MQ=255)
                                                                                                                           gCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGt                         <  1:242874/148‑1 (MQ=255)
                                                                                                                            cATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGCCACCGGTNGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGtt                        <  1:1811646/148‑1 (MQ=255)
                                                                                                                              tCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTAc                      >  2:706516/1‑148 (MQ=255)
                                                                                                                              tCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTAc                      <  1:387090/148‑1 (MQ=255)
                                                                                                                                 tGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACttt                   <  2:383908/148‑1 (MQ=255)
                                                                                                                                  gCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTttg                  >  1:45077/1‑148 (MQ=255)
                                                                                                                                  gCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTttg                  <  2:2330795/148‑1 (MQ=255)
                                                                                                                                  gCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTttg                  <  2:242509/148‑1 (MQ=255)
                                                                                                                                     cGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTttgttg               >  1:2342187/1‑148 (MQ=255)
                                                                                                                                      gCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGa              >  1:2298871/1‑148 (MQ=255)
                                                                                                                                        tGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGc            <  1:817745/148‑1 (MQ=255)
                                                                                                                                        tGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGc            >  2:2539993/1‑148 (MQ=255)
                                                                                                                                           tAAGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAt         <  1:2072873/148‑1 (MQ=255)
                                                                                                                                           tAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAt         >  2:1482291/1‑148 (MQ=255)
                                                                                                                                           tAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAt         >  1:2451776/1‑148 (MQ=255)
                                                                                                                                           tAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAt         >  1:2451775/1‑148 (MQ=255)
                                                                                                                                            aaGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAtt        >  2:1585570/1‑148 (MQ=255)
                                                                                                                                            aaGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAtt        <  1:2933723/148‑1 (MQ=255)
                                                                                                                                             aGCGTTATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAttt       >  1:314617/1‑148 (MQ=255)
                                                                                                                                               cGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCAttttt     <  2:1065122/148‑1 (MQ=255)
                                                                                                                                                gTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATTTTTg    >  1:1640488/1‑148 (MQ=255)
                                                                                                                                                  gATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATTTTTGcc  >  1:3397707/1‑148 (MQ=255)
                                                                                                                                                  gATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATGTTTGcc  >  1:3399179/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTTTCCAGCAATTGCCCCTGCTGTTTCGCGGCGTTGGGTTTATTGACCGAAATGCCCCATGTTGTGCGCTCAACGTCGGCTCCCGGTGGCAGCGGTACTTCCACCTGACCGTAGCGCAGAACTGCATCCTGCTCGCTGGTAAGCGTGATTTCATCGAGATACAGCGCATCGCTGTCAATCTCATTGCTGGTCACCGGTTGCAGAATAAAGCTCATCTCTTCTTCACCAGGGATAAGCCGATACAACTGGCGTTCAACGGTCACCGGAATGTTACTTTGTTGAGCCATTTTTGCC  >  CP009273/2323476‑2323769

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: