Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,182,647 G→T 92.1% pseudogene (1868/2445 nt) yqiG → pseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,182,6470GT92.1% 587.4 / 42.3 216pseudogene (1868/2445 nt)yqiGpseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (7/10);  new base T (106/93);  total (113/103)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 220 total aligned reads displayed.

ATGACTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAACAAT  >  CP009273/3182503‑3182792
                                                                                                                                                |                                                                                                                                                 
aTGACTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCtaga                                                                                                                                                >  1:2946733/1‑148 (MQ=255)
  gACTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACa                                                                                                                                              <  1:669569/148‑1 (MQ=255)
   aCTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAg                                                                                                                                             >  1:200453/1‑148 (MQ=255)
   aCTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAg                                                                                                                                             <  2:1630407/148‑1 (MQ=255)
    ctctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGt                                                                                                                                            >  1:1819497/1‑148 (MQ=255)
    ctctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGt                                                                                                                                            <  2:814592/148‑1 (MQ=255)
    ctctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGt                                                                                                                                            >  2:329300/1‑148 (MQ=255)
     tctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGtt                                                                                                                                           <  2:916486/148‑1 (MQ=255)
     tctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGtt                                                                                                                                           >  1:983824/1‑148 (MQ=255)
     tctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGtt                                                                                                                                           <  2:1616978/148‑1 (MQ=255)
     tctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGtt                                                                                                                                           <  1:1643845/148‑1 (MQ=255)
      gttCCGATCGAAATAATTCGTGGAGCGTTTCTGTTTCTGGGGACAATGATGAATTCAAAGTTTTGCAAGCGTCACTACGCGCCAGTTATCAGCTTAATACCGAGAACGGTCGACTCTAGCTCTCCGGTACATCACAGCTAGACAGtta                                                                                                                                          <  1:918578/147‑1 (MQ=255)
         ccGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGttatta                                                                                                                                       <  1:2173959/148‑1 (MQ=255)
          cGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTtattat                                                                                                                                      <  1:1680739/148‑1 (MQ=255)
             tCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATtct                                                                                                                                   >  2:848750/1‑148 (MQ=255)
             tCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATtct                                                                                                                                   >  1:1130624/1‑148 (MQ=255)
              cGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATtctc                                                                                                                                  <  2:49332/148‑1 (MQ=255)
              cGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATtctc                                                                                                                                  >  2:2551298/1‑148 (MQ=255)
                 aataatTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGa                                                                                                                               <  1:85574/148‑1 (MQ=255)
                  ataatTCCTGGAGCGTTTCTGCTTCAGGCGACAATTATGAATTCAAATATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGaa                                                                                                                              <  1:3294074/148‑1 (MQ=255)
                  ataatTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGaa                                                                                                                              >  2:2398466/1‑148 (MQ=255)
                  ataatTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGaa                                                                                                                              <  1:1878346/148‑1 (MQ=255)
                      ttCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGcc                                                                                                                          <  2:2449633/148‑1 (MQ=255)
                         cTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGt                                                                                                                       <  1:1333319/148‑1 (MQ=255)
                            gAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTgg                                                                                                                    >  2:3228332/1‑148 (MQ=255)
                             aGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGa                                                                                                                   <  2:1043358/148‑1 (MQ=255)
                              gCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGaa                                                                                                                  <  2:2219681/148‑1 (MQ=255)
                               cGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGGTGGAAt                                                                                                                 >  2:3156940/1‑148 (MQ=255)
                                     tGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttca                                                                                                           >  1:2550596/1‑148 (MQ=255)
                                      gCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcat                                                                                                          >  2:2770966/1‑148 (MQ=255)
                                        ttCAGGCGACAATGCTGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcattc                                                                                                        <  2:302490/148‑1 (MQ=255)
                                         tCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcattca                                                                                                       <  1:2669300/148‑1 (MQ=255)
                                                aCAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGAc                                                                                                >  2:1904272/1‑148 (MQ=255)
                                                aCAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGAc                                                                                                >  2:612016/1‑148 (MQ=255)
                                                aCAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGAc                                                                                                >  1:2922814/1‑148 (MQ=255)
                                                 cAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACt                                                                                               >  2:2944720/1‑148 (MQ=255)
                                                  aaTGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTc                                                                                              >  1:1899849/1‑148 (MQ=255)
                                                   atgatgAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCg                                                                                             <  2:172466/148‑1 (MQ=255)
                                                      atgaATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCa                                                                                          >  2:1108280/1‑148 (MQ=255)
                                                        gaATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACg                                                                                        >  1:2233532/1‑148 (MQ=255)
                                                            tCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGc                                                                                    >  2:1325642/1‑148 (MQ=255)
                                                            tCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGc                                                                                    >  1:1008498/1‑148 (MQ=255)
                                                             cAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGcc                                                                                   >  1:1533219/1‑148 (MQ=255)
                                                                   tatGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTc                                                                             <  2:2300163/148‑1 (MQ=255)
                                                                    atGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTcc                                                                            >  1:692024/1‑148 (MQ=255)
                                                                    atGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTcc                                                                            >  1:1271519/1‑148 (MQ=255)
                                                                    atGAAAGCGCCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTcc                                                                            >  1:2764702/1‑148 (MQ=255)
                                                                      gAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCAc                                                                          >  1:517918/1‑148 (MQ=255)
                                                                       aaaGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGACACGGTGCCGCTTTCCACg                                                                         >  2:3267197/1‑148 (MQ=255)
                                                                       aaaGCGTCACTACGCGCCAGTTATCAGCAAAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACg                                                                         >  1:2813946/1‑148 (MQ=255)
                                                                          gCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAAGGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACt                                                                      >  2:865525/1‑148 (MQ=255)
                                                                           cGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACta                                                                     >  1:2413966/1‑148 (MQ=255)
                                                                             tCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACtata                                                                   >  1:3024706/1‑148 (MQ=255)
                                                                               aCTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACATCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGc                                                                 >  1:2880415/1‑148 (MQ=255)
                                                                                 tACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCgg                                                               <  1:752778/148‑1 (MQ=255)
                                                                                    gcgcCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAg                                                            <  1:2637180/148‑1 (MQ=255)
                                                                                    gcgcCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAg                                                            <  1:540264/148‑1 (MQ=255)
                                                                                      gcCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTg                                                          <  2:2957661/148‑1 (MQ=255)
                                                                                       ccAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGc                                                         >  1:2181697/1‑148 (MQ=255)
                                                                                           ttATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGAc                                                     >  1:1247010/1‑148 (MQ=255)
                                                                                            tATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACt                                                    <  2:1449952/148‑1 (MQ=255)
                                                                                              tCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTcg                                                  <  2:892320/148‑1 (MQ=255)
                                                                                                 gCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGt                                               >  1:3139533/1‑148 (MQ=255)
                                                                                                  cATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGtt                                              >  2:2113509/1‑148 (MQ=255)
                                                                                                   ataataCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGttt                                             >  2:1952832/1‑148 (MQ=255)
                                                                                                    taataCCGAGACCGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGtttt                                            <  2:2838742/148‑1 (MQ=255)
                                                                                                    taataCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGtttt                                            >  2:97243/1‑148 (MQ=255)
                                                                                                     aataCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTa                                           >  1:1907259/1‑148 (MQ=255)
                                                                                                     aataCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTa                                           <  1:2499773/148‑1 (MQ=255)
                                                                                                     aataCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTa                                           <  1:1108280/148‑1 (MQ=255)
                                                                                                       taCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATg                                         >  2:2004242/1‑148 (MQ=255)
                                                                                                       taCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATg                                         >  1:3079169/1‑148 (MQ=255)
                                                                                                           gagaACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCg                                     <  2:1819497/148‑1 (MQ=255)
                                                                                                            agaAGGGTCGTCTCTACCTATCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTAATTCACTGCGACTCGCCTCGGTGCGGCTTTCCACGAATATAGCGGTAGTGCTGACTCGCGTTTTATGATCGa                                    <  2:2813946/148‑1 (MQ=255)
                                                                                                            agaACGGCCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGa                                    <  2:1548711/148‑1 (MQ=255)
                                                                                                             gaACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGAc                                   >  1:2102020/1‑148 (MQ=255)
                                                                                                                     gACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACgg                           <  2:1130624/148‑1 (MQ=255)
                                                                                                                       ctctACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCa                         >  2:506051/1‑148 (MQ=255)
                                                                                                                       ctctACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCa                         <  1:1745348/148‑1 (MQ=255)
                                                                                                                         ctACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACt                       <  2:1215634/148‑1 (MQ=255)
                                                                                                                           aCCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGa                     <  2:983823/148‑1 (MQ=255)
                                                                                                                           aCCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGa                     <  2:1224861/148‑1 (MQ=255)
                                                                                                                           aCCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGa                     <  1:656473/148‑1 (MQ=255)
                                                                                                                            ccTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGaa                    >  1:1076857/1‑148 (MQ=255)
                                                                                                                            ccTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGaa                    >  1:957938/1‑148 (MQ=255)
                                                                                                                               ctcCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGat                 >  1:1563135/1‑148 (MQ=255)
                                                                                                                                tcCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGata                <  1:1858839/148‑1 (MQ=255)
                                                                                                                                tcCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGata                <  1:492515/148‑1 (MQ=255)
                                                                                                                                 ccGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGatat               >  2:2681747/1‑148 (MQ=255)
                                                                                                                                 ccGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGatat               <  1:1766900/148‑1 (MQ=255)
                                                                                                                                 ccGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGatat               >  1:750364/1‑148 (MQ=255)
                                                                                                                                     tACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCg           >  1:1000775/1‑148 (MQ=255)
                                                                                                                                     tACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCg           >  2:2866360/1‑148 (MQ=255)
                                                                                                                                     tACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCg           <  2:2826638/148‑1 (MQ=255)
                                                                                                                                       cATCACAGCTAGACAGTTATTATTCTCTGATTGCCAGTTGTAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGCCTATAGCGGTAGTGCTGACTCGCGTTTTTTGATCGACGCAGACGGCACTGAAGATATTCCGtt         <  2:3275778/148‑1 (MQ=255)
                                                                                                                                       cATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGtt         <  1:406077/148‑1 (MQ=255)
                                                                                                                                          cacaGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGaa      >  2:424726/1‑148 (MQ=255)
                                                                                                                                           acaGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAAc     >  2:2389476/1‑148 (MQ=255)
                                                                                                                                             aGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAACaa   >  2:1125535/1‑148 (MQ=255)
                                                                                                                                              gCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAACAAt  <  1:2202027/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ATGACTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAACAAT  >  CP009273/3182503‑3182792

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: