Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,511,828 T→A 5.1% intergenic (‑118/‑118) alaW ← / → yfeC tRNA‑Ala/DUF1323 family putative DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,511,8280TA 5.1% 558.3 / 15.6 177intergenic (‑118/‑118)alaW/yfeCtRNA‑Ala/DUF1323 family putative DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (73/95);  new base A (6/3);  total (79/98)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 189 total aligned reads displayed.

AAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACC  >  CP009273/2511682‑2511975
                                                                                                                                                  |                                                                                                                                                   
  gcgcTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGa                                                                                                                                                  >  1:2347251/1‑148 (MQ=255)
  gcgcTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGa                                                                                                                                                  >  1:166600/1‑148 (MQ=255)
  gcgcTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGa                                                                                                                                                  <  2:1958984/148‑1 (MQ=255)
     ctctCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGaaaa                                                                                                                                               <  1:2567838/148‑1 (MQ=255)
     ctctCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGaaaa                                                                                                                                               >  1:64897/1‑148 (MQ=255)
      tctcCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAAc                                                                                                                                              >  2:618039/1‑148 (MQ=255)
      tctcCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAAc                                                                                                                                              <  1:1699406/148‑1 (MQ=255)
        tcCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCg                                                                                                                                            <  1:3441774/148‑1 (MQ=255)
                gAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCCTCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCAtt                                                                                                                                    >  1:3243768/1‑148 (MQ=255)
                  gCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCa                                                                                                                                  <  2:2104413/148‑1 (MQ=255)
                        gCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGtt                                                                                                                            <  1:860660/148‑1 (MQ=255)
                            cGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTc                                                                                                                        <  1:3790094/148‑1 (MQ=255)
                            cGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTc                                                                                                                        >  1:798654/1‑148 (MQ=255)
                                      aGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCAGAAAAAACAGGCAAATGAACAATTTTGGAAAACCGCTCGCAATCAGTAGTTATTCATGTCACGGt                                                                                                              >  2:1213788/1‑148 (MQ=255)
                                       gCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGtt                                                                                                             <  1:113303/148‑1 (MQ=255)
                                        cTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTTGAAAACCGCTCGCATTCAGTAGTTATTAATGTCACGGttt                                                                                                            >  1:97302/1‑148 (MQ=255)
                                        cTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGttt                                                                                                            >  2:525883/1‑148 (MQ=255)
                                         ttGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTc                                                                                                           <  2:2856477/148‑1 (MQ=255)
                                         ttGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTc                                                                                                           <  2:2594002/148‑1 (MQ=255)
                                         ttGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTc                                                                                                           <  2:3313458/148‑1 (MQ=255)
                                          tGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTcc                                                                                                          >  2:4572/1‑148 (MQ=255)
                                               aGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGTAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTaa                                                                                                     <  2:358837/148‑1 (MQ=255)
                                                  tGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACTGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAgt                                                                                                  <  1:3024843/148‑1 (MQ=255)
                                                   gACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAgtg                                                                                                 <  1:3774252/148‑1 (MQ=255)
                                                     cGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAgtggt                                                                                               <  2:1544985/148‑1 (MQ=255)
                                                       ggcggcATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGgtgt                                                                                             >  1:2535717/1‑148 (MQ=255)
                                                         cggcATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTta                                                                                           <  1:2473315/148‑1 (MQ=255)
                                                          ggcATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTtat                                                                                          <  1:675408/148‑1 (MQ=255)
                                                          ggcATCATATGAAGTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTtat                                                                                          <  2:3211164/148‑1 (MQ=255)
                                                                atatGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATg                                                                                    <  1:1948382/148‑1 (MQ=255)
                                                                atatGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATg                                                                                    <  1:15685/148‑1 (MQ=255)
                                                                 tatGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGa                                                                                   <  1:3589287/148‑1 (MQ=255)
                                                                  atGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGaa                                                                                  >  1:1071216/1‑148 (MQ=255)
                                                                  atGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGaa                                                                                  <  2:866058/148‑1 (MQ=255)
                                                                  atGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGa                                                                                   >  2:3530183/1‑147 (MQ=255)
                                                                    gAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAact                                                                                <  1:482640/148‑1 (MQ=255)
                                                                    gAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAact                                                                                <  1:963086/148‑1 (MQ=255)
                                                                     aaTTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTAGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAacta                                                                               >  1:3714359/1‑148 (MQ=255)
                                                                      aTTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAactac                                                                              <  2:185925/148‑1 (MQ=255)
                                                                       ttCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAactact                                                                             >  1:378644/1‑148 (MQ=255)
                                                                         ccGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTaa                                                                           <  2:3023648/148‑1 (MQ=255)
                                                                          cGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAAt                                                                          <  2:2106504/148‑1 (MQ=255)
                                                                             ccGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGa                                                                       <  2:1987072/148‑1 (MQ=255)
                                                                              cGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGAc                                                                      >  2:2245570/1‑148 (MQ=255)
                                                                              cGAATGTGTCAACGGCAAAATGCAACGTGTAGTTTCAATCGCTGAAAAAACACGCAAATGAACAATTTTGGAAAACCGCTCGCATTAAGTAGATAATCATGTCACGGTTTCCTGTAAAGTCGTGTTATAAAAAGAACTACTAATAGAc                                                                      >  2:2303220/1‑148 (MQ=255)
                                                                                  tgtgTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCcac                                                                  <  2:2961582/148‑1 (MQ=255)
                                                                                  tgtgTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCcac                                                                  <  2:2431172/148‑1 (MQ=255)
                                                                                   gtgtCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCcaca                                                                 <  1:912576/148‑1 (MQ=255)
                                                                                         aCGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACAtt                                                           >  1:972453/1‑148 (MQ=255)
                                                                                          cGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTc                                                          <  1:1027050/148‑1 (MQ=255)
                                                                                          cGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTc                                                          <  1:2724761/148‑1 (MQ=255)
                                                                                          cGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTc                                                          <  1:2938113/148‑1 (MQ=255)
                                                                                          cGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTAGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTc                                                          >  1:3183059/1‑148 (MQ=255)
                                                                                             cAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAgg                                                       >  1:2184290/1‑148 (MQ=255)
                                                                                               aaTTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTAGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGa                                                     >  1:28375/1‑148 (MQ=255)
                                                                                               aaTTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTAGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGa                                                     >  1:27680/1‑148 (MQ=255)
                                                                                                aTTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGaa                                                    <  2:3616458/148‑1 (MQ=255)
                                                                                                  tGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAAtt                                                  >  2:3171418/1‑148 (MQ=255)
                                                                                                   gCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTg                                                 <  1:2037659/148‑1 (MQ=255)
                                                                                                   gCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTg                                                 <  1:1948861/148‑1 (MQ=255)
                                                                                                   gCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTg                                                 <  2:798654/148‑1 (MQ=255)
                                                                                                     aaCGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGtt                                               <  2:3355097/148‑1 (MQ=255)
                                                                                                        gtgtAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTAGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATg                                            <  2:1503247/148‑1 (MQ=255)
                                                                                                         tgtAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGt                                           >  2:1780984/1‑148 (MQ=255)
                                                                                                           tAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTc                                         >  1:2811885/1‑148 (MQ=255)
                                                                                                           tAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTc                                         >  2:2470403/1‑148 (MQ=255)
                                                                                                           tAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTc                                         >  2:2546845/1‑148 (MQ=255)
                                                                                                                  aaTCGCTGAAAAATGAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGgaga                                  <  2:2759817/148‑1 (MQ=255)
                                                                                                                  aaTCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACTTTCAGGGAATTGTTATGTTCAAGgaga                                  <  2:1047455/148‑1 (MQ=255)
                                                                                                                   aTCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGagag                                 <  2:2347251/148‑1 (MQ=255)
                                                                                                                    tCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAgg                                <  2:520254/148‑1 (MQ=255)
                                                                                                                    tCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAgg                                <  2:557491/148‑1 (MQ=255)
                                                                                                                     cGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGa                               >  1:1235219/1‑148 (MQ=255)
                                                                                                                     cGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGa                               >  1:2641086/1‑148 (MQ=255)
                                                                                                                     cGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGa                               <  2:2518206/148‑1 (MQ=255)
                                                                                                                         gAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGac                           >  2:1200872/1‑148 (MQ=255)
                                                                                                                          aaaaaTCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGaca                          <  1:2139757/148‑1 (MQ=255)
                                                                                                                             aaTCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCa                       >  2:3452664/1‑148 (MQ=255)
                                                                                                                              aTCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAg                      <  1:618039/148‑1 (MQ=255)
                                                                                                                               tCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGa                     >  1:1294482/1‑148 (MQ=255)
                                                                                                                               tCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGa                     >  2:833969/1‑148 (MQ=255)
                                                                                                                                   gCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGaa                 >  1:684177/1‑148 (MQ=255)
                                                                                                                                    cAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAAc                <  2:3209240/148‑1 (MQ=255)
                                                                                                                                      aaTGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACtt              >  1:228551/1‑148 (MQ=255)
                                                                                                                                       aTGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTg             >  2:3557859/1‑148 (MQ=255)
                                                                                                                                       aTGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTg             >  2:50086/1‑148 (MQ=255)
                                                                                                                                        tGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGc            >  1:291689/1‑148 (MQ=255)
                                                                                                                                        tGAACAATTTAGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGc            <  1:1621243/148‑1 (MQ=255)
                                                                                                                                         gAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGcc           >  1:55537/1‑148 (MQ=255)
                                                                                                                                         gAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGcc           >  2:391158/1‑148 (MQ=255)
                                                                                                                                          aaCAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAATCTGCTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCa          >  1:3747558/1‑148 (MQ=255)
                                                                                                                                            cAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGa        <  1:2674848/148‑1 (MQ=255)
                                                                                                                                             aaTTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGAc       <  2:1425672/148‑1 (MQ=255)
                                                                                                                                              aTTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACt      >  2:1371631/1‑148 (MQ=255)
                                                                                                                                              aTTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACt      >  2:1301301/1‑148 (MQ=255)
                                                                                                                                               ttttGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTg     <  1:2026998/148‑1 (MQ=255)
                                                                                                                                                tttGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGa    >  1:2677129/1‑148 (MQ=255)
                                                                                                                                                tttGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGa    >  1:1773143/1‑148 (MQ=255)
                                                                                                                                                tttGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGa    >  2:195353/1‑148 (MQ=255)
                                                                                                                                                tttGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGa    >  1:1773155/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAAAATGAACTACTAATAGACCCACATACATTCAGGGAATTGTTATGTTCAAGGAGAGGATGACACCAGATGAACTTGCCAGACTGACC  >  CP009273/2511682‑2511975

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: