Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,087,484 A→T 6.1% Y88N (TAC→AAC)  pgaA ← biofilm adhesin polysaccharide PGA secretin; OM porin; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine export protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,087,4840AT 6.1% 659.0 / 24.8 211Y88N (TAC→AAC) pgaAbiofilm adhesin polysaccharide PGA secretin; OM porin; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine export protein
Reads supporting (aligned to +/- strand):  ref base A (94/104);  new base T (2/11);  total (96/115)
Fisher's exact test for biased strand distribution p-value = 4.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 219 total aligned reads displayed.

CCCGTTGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAATGCATC  >  CP009273/1087337‑1087631
                                                                                                                                                   |                                                                                                                                                   
cccGTTGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTa                                                                                                                                                     >  2:3251956/1‑148 (MQ=255)
  cGTTGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTaaa                                                                                                                                                   >  2:1740443/1‑148 (MQ=255)
  cGTTGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTaaa                                                                                                                                                   >  2:1741562/1‑148 (MQ=255)
    ttGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACg                                                                                                                                                 >  2:1367947/1‑148 (MQ=255)
    ttGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACg                                                                                                                                                 >  2:3047107/1‑148 (MQ=255)
     tGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTTAACgg                                                                                                                                                >  1:1899724/1‑148 (MQ=255)
        tAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGtaa                                                                                                                                             >  1:783485/1‑148 (MQ=255)
        tAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGtaa                                                                                                                                             >  1:1049474/1‑148 (MQ=255)
        tAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGtaa                                                                                                                                             >  1:1251448/1‑148 (MQ=255)
           tCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGtaataa                                                                                                                                          >  1:2526121/1‑148 (MQ=255)
            ccTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTTAACGGTAATAAc                                                                                                                                         <  1:581623/148‑1 (MQ=255)
            ccTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAAc                                                                                                                                         <  1:1586107/148‑1 (MQ=255)
             cTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAAcc                                                                                                                                        >  2:143723/1‑148 (MQ=255)
              ttATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCt                                                                                                                                       >  1:3024450/1‑148 (MQ=255)
               tATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTg                                                                                                                                      <  2:3763533/148‑1 (MQ=255)
                    tGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTAt                                                                                                                                 <  2:3147309/148‑1 (MQ=255)
                    tGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTAt                                                                                                                                 <  2:731244/148‑1 (MQ=255)
                     gCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATc                                                                                                                                <  2:1454530/148‑1 (MQ=255)
                     gCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATc                                                                                                                                <  2:1238366/148‑1 (MQ=255)
                      cGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCt                                                                                                                               <  1:3497771/148‑1 (MQ=255)
                        gCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTTAACGGTAATAACCTGTTTATCTTg                                                                                                                             <  1:1859911/148‑1 (MQ=255)
                          tCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGcc                                                                                                                           <  1:626275/148‑1 (MQ=255)
                           ccAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGccc                                                                                                                          >  2:915303/1‑148 (MQ=255)
                           ccAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGGTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGccc                                                                                                                          >  1:3570065/1‑148 (MQ=255)
                               tgcgagCGCCTTTTGCCACAGTTTAAGCTTTTTTTGCCTTTTTGGCTGTTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGcc                                                                                                                      <  1:1305777/145‑1 (MQ=255)
                               agagagCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGcc                                                                                                                      >  1:854632/1‑148 (MQ=255)
                               agagagCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGcc                                                                                                                      >  2:1659375/1‑148 (MQ=255)
                               agagagCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGcc                                                                                                                      >  1:3608330/1‑148 (MQ=255)
                                   agCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATa                                                                                                                  >  1:1733091/1‑148 (MQ=255)
                                    gcgcCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATaa                                                                                                                 <  2:2152359/148‑1 (MQ=255)
                                    gcgcCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATaa                                                                                                                 <  2:1480565/148‑1 (MQ=255)
                                    gcgcCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATaa                                                                                                                 <  1:2855648/148‑1 (MQ=255)
                                        cTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCa                                                                                                             >  2:3314804/1‑148 (MQ=255)
                                        cTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCa                                                                                                             <  2:834324/148‑1 (MQ=255)
                                        cTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCa                                                                                                             >  2:433476/1‑148 (MQ=255)
                                        cTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCa                                                                                                             >  2:433408/1‑148 (MQ=255)
                                         ttttGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCaa                                                                                                            <  1:2234069/148‑1 (MQ=255)
                                          tttGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAAt                                                                                                           <  1:3234252/148‑1 (MQ=255)
                                          tttGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAAt                                                                                                           <  2:2998282/148‑1 (MQ=255)
                                          tttGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAAt                                                                                                           <  1:2495787/148‑1 (MQ=255)
                                          tttGCCACAGTGTAAGCGAGTTTTGCAATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTTAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAAt                                                                                                           <  2:2786852/148‑1 (MQ=255)
                                              ccACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGt                                                                                                       >  1:1225839/1‑148 (MQ=255)
                                              ccACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGt                                                                                                       >  1:590212/1‑148 (MQ=255)
                                                    tgtAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAg                                                                                                 >  1:1505146/1‑148 (MQ=255)
                                                      tAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTc                                                                                               <  2:2439470/148‑1 (MQ=255)
                                                          cGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCa                                                                                           >  1:3212129/1‑148 (MQ=255)
                                                           gAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCaa                                                                                          <  1:1287900/148‑1 (MQ=255)
                                                           gAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCaa                                                                                          <  1:3066991/148‑1 (MQ=255)
                                                            aGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAAt                                                                                         <  1:1800799/148‑1 (MQ=255)
                                                            aGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAAt                                                                                         <  2:505124/148‑1 (MQ=255)
                                                                 tGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGt                                                                                    <  2:283299/148‑1 (MQ=255)
                                                                    cATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTAt                                                                                 <  1:3711411/148‑1 (MQ=255)
                                                                    cATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTAt                                                                                 <  1:2678902/148‑1 (MQ=255)
                                                                      ttgttgCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTg                                                                               <  1:292291/148‑1 (MQ=255)
                                                                       tgttgCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGc                                                                              >  1:1486664/1‑148 (MQ=255)
                                                                         ttgCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTg                                                                            <  2:2627111/148‑1 (MQ=255)
                                                                          tgCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGa                                                                           <  1:918627/148‑1 (MQ=255)
                                                                                 tACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGa                                                                    >  1:2606903/1‑148 (MQ=255)
                                                                                 tACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGa                                                                    >  1:3654947/1‑148 (MQ=255)
                                                                                  aCGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGAt                                                                   <  2:3060782/148‑1 (MQ=255)
                                                                                  aCGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGAt                                                                   >  2:1483225/1‑148 (MQ=255)
                                                                                   cGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGAtt                                                                  >  2:3563687/1‑148 (MQ=255)
                                                                                   cGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGAtt                                                                  <  2:1689939/148‑1 (MQ=255)
                                                                                   cGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGAtt                                                                  >  1:3720375/1‑148 (MQ=255)
                                                                                     aTAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGAtttt                                                                >  2:3557882/1‑148 (MQ=255)
                                                                                     aTAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGAtttt                                                                <  1:3406618/148‑1 (MQ=255)
                                                                                      tAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTa                                                               >  2:1618947/1‑148 (MQ=255)
                                                                                         gCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTg                                                            >  2:1408656/1‑148 (MQ=255)
                                                                                            aCGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCaa                                                         <  2:3596354/148‑1 (MQ=255)
                                                                                                  aCAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATg                                                   >  2:470277/1‑148 (MQ=255)
                                                                                                  aCAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATg                                                   >  2:41274/1‑148 (MQ=255)
                                                                                                  aCAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATg                                                   >  2:2392006/1‑148 (MQ=255)
                                                                                                  aCAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATg                                                   <  1:1668492/148‑1 (MQ=255)
                                                                                                      cTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACaa                                               <  1:3314804/148‑1 (MQ=255)
                                                                                                      cTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACaa                                               >  2:869050/1‑148 (MQ=255)
                                                                                                      cTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACaa                                               <  1:433450/148‑1 (MQ=255)
                                                                                                         cATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgc                                            <  2:2176178/148‑1 (MQ=255)
                                                                                                           tAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctg                                          >  1:3334760/1‑148 (MQ=255)
                                                                                                            aaCCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctgg                                         <  1:408671/148‑1 (MQ=255)
                                                                                                            aaCCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctgg                                         <  2:3285110/148‑1 (MQ=255)
                                                                                                             aCCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctggc                                        >  1:818691/1‑148 (MQ=255)
                                                                                                              ccACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAgctggct                                       >  1:2573200/1‑148 (MQ=255)
                                                                                                                 cgcgCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAg                                    <  1:1659309/148‑1 (MQ=255)
                                                                                                                  gcgcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGt                                   <  1:98429/148‑1 (MQ=255)
                                                                                                                   cgcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTa                                  <  1:111295/148‑1 (MQ=255)
                                                                                                                    gcTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTAt                                 <  1:530183/148‑1 (MQ=255)
                                                                                                                          aaTTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCt                           >  1:2463725/1‑148 (MQ=255)
                                                                                                                          aaTTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCt                           >  2:1850008/1‑148 (MQ=255)
                                                                                                                             tGCTGATGACGGTAGCGGTTGTTAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGc                        >  2:753633/1‑148 (MQ=255)
                                                                                                                              gCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCg                       <  2:1891425/148‑1 (MQ=255)
                                                                                                                                 gatgaCGGTAGCGGTTGTTAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGc                    <  2:347462/148‑1 (MQ=255)
                                                                                                                                         tAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCaataa            >  2:3349596/1‑148 (MQ=255)
                                                                                                                                           gCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATc          <  1:1740743/148‑1 (MQ=255)
                                                                                                                                            cGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCa         >  1:2193959/1‑148 (MQ=255)
                                                                                                                                              gTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAAt       <  2:1049474/148‑1 (MQ=255)
                                                                                                                                              gTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAAt       <  1:2107643/148‑1 (MQ=255)
                                                                                                                                                 gTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAATGca    <  2:1645263/148‑1 (MQ=255)
                                                                                                                                                  tAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAATGcat   <  2:2965621/148‑1 (MQ=255)
                                                                                                                                                   aaaCGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAATGcatc  <  1:502224/148‑1 (MQ=255)
                                                                                                                                                   aaaCGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAATGcatc  >  1:2723808/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCCGTTGATAATCCTTATTTTGCGGCTCCAGAGAGAGCGCCTTTTGCCACAGTGTAAGCGAGTTTTGCCATTGTTGCAGGTTACGATAAGCGACGGCGACAGCTGCATAACCACGCGCTGGTAATTGCTGATGACGGTAGCGGTTGTAAACGGTAATAACCTGTTTATCTTGCCCGGCCCATAAGGCAATCTGTAACCAGTCAGCAATTTGGTTATTGCTGAGTGCTGATTTTAGTGCAAACCATGACAAAGCTGGCTGAGTATTACCCTTGCGAGCTTCAATAATCAATGCATC  >  CP009273/1087337‑1087631

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: