Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,087,715 C→T 77.2% M44I (ATG→ATA pgaC ← biofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,087,7150CT77.2% 281.0 / 82.8 153M44I (ATG→ATApgaCbiofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase
Reads supporting (aligned to +/- strand):  ref base C (18/17);  new base T (66/52);  total (84/69)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

TAGACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAG  >  NC_000913/1087618‑1087810
                                                                                                 |                                                                                                
  gACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATaa                                                                                              >  2:975002/1‑100 (MQ=255)
  gACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATaa                                                                                              >  2:1332730/1‑100 (MQ=255)
  gACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATaa                                                                                              >  1:621117/1‑100 (MQ=255)
    cGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATg                                                                                            <  1:2532086/100‑1 (MQ=255)
     ggATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATgg                                                                                           <  1:4067359/100‑1 (MQ=255)
       aTTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGAc                                                                                         <  1:3351480/100‑1 (MQ=255)
         tATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACAt                                                                                       >  2:623995/1‑100 (MQ=255)
            cTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATaaa                                                                                    <  2:1117389/100‑1 (MQ=255)
             tttCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATaaaa                                                                                   >  1:958968/1‑100 (MQ=255)
             tttCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATaaaa                                                                                   >  2:1644311/1‑100 (MQ=255)
             tttCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATaaaa                                                                                   >  1:3350166/1‑100 (MQ=255)
                 aaCTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACg                                                                               <  2:595804/100‑1 (MQ=255)
                   cttccgatcTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGc                                                                             <  2:3708883/92‑1 (MQ=255)
                   cTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGACGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGc                                                                             >  2:2018969/1‑100 (MQ=255)
                    tGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGcc                                                                            <  2:1122089/100‑1 (MQ=255)
                     ggggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCa                                                                           >  2:1787323/1‑100 (MQ=255)
                     ggggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCa                                                                           <  1:2920334/100‑1 (MQ=255)
                        gAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCagaa                                                                        >  1:2450601/1‑100 (MQ=255)
                           cTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAgaagaa                                                                     >  2:2749229/1‑100 (MQ=255)
                             ggTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGcagatcggaag                                                                   >  1:3708910/1‑90 (MQ=255)
                             ggTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGaaaa                                                                   <  1:1322705/100‑1 (MQ=255)
                             ggTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGaaaa                                                                   >  2:2772139/1‑100 (MQ=255)
                              gTGCGTTTTCTCCCCACGACCAGTGGCGTTCACGATAGACACAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAAc                                                                  <  2:1144074/100‑1 (MQ=255)
                              gTGCGGTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCAAGAAATAGGCGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGGAAAc                                                                  <  2:1202317/100‑1 (MQ=255)
                               tGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACg                                                                 >  2:249145/1‑100 (MQ=255)
                               tGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAAc                                                                  <  2:1329732/99‑1 (MQ=255)
                                  gTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAAc                                                              >  2:2391059/1‑100 (MQ=255)
                                      tctcCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtca                                                          <  1:130778/100‑1 (MQ=255)
                                      tctcCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtca                                                          >  2:2589896/1‑100 (MQ=255)
                                      tctcCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtca                                                          <  1:3785781/100‑1 (MQ=255)
                                      tctcCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtca                                                          <  2:90986/100‑1 (MQ=255)
                                        tcCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatc                                                        <  1:2510765/100‑1 (MQ=255)
                                        tcCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatc                                                        >  1:3310344/1‑100 (MQ=255)
                                         ccccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatca                                                       >  1:2473808/1‑100 (MQ=255)
                                         ccccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatca                                                       >  2:4009310/1‑100 (MQ=255)
                                         ccccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatca                                                       >  1:2878055/1‑100 (MQ=255)
                                          cccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatcat                                                      >  2:1007575/1‑100 (MQ=255)
                                           ccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatcatc                                                     <  1:2391074/100‑1 (MQ=255)
                                             aCGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAAcc                                                             <  2:2079217/90‑1 (MQ=255)
                                             aCGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAg                                                   <  2:3876794/100‑1 (MQ=255)
                                              cGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGt                                                  <  2:1701961/100‑1 (MQ=255)
                                                    gTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCa                                            >  1:773106/1‑100 (MQ=255)
                                                       gCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCagag                                         <  1:3178255/100‑1 (MQ=255)
                                                        ccctCACGATAGACCCAGAAATAGACGCCGCCAACACTCCATATAATGCAAATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCag                                           <  2:4001712/94‑1 (MQ=255)
                                                        cGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGt                                        >  2:3141226/1‑100 (MQ=255)
                                                        cGTTCACCATAGACCCATAAATAGACGCCTCCAATAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCagag                                         <  2:2874787/99‑1 (MQ=255)
                                                         gTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGt                                        >  1:1751549/1‑99 (MQ=255)
                                                           tCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGaa                                     <  2:333676/100‑1 (MQ=255)
                                                            cACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGaaa                                    <  1:3598453/100‑1 (MQ=255)
                                                            cACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGaaa                                    <  1:10231/100‑1 (MQ=255)
                                                             aagATAGACCGAGAAATAGACGCCGCCAACAATCCATATAATGGAGATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTg                                       <  2:2624920/94‑1 (MQ=255)
                                                             aCGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAg                                   <  1:1787333/100‑1 (MQ=255)
                                                             aCGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAg                                   <  1:1838853/100‑1 (MQ=255)
                                                              cGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGt                                  >  1:1163053/1‑100 (MQ=255)
                                                              cGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGt                                  >  1:3047198/1‑100 (MQ=255)
                                                                   gACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgct                             >  1:2118074/1‑100 (MQ=255)
                                                                   gACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgct                             >  2:2849625/1‑100 (MQ=255)
                                                                   gACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgct                             <  1:263787/100‑1 (MQ=255)
                                                                   gACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgct                             >  1:3943978/1‑100 (MQ=255)
                                                                    aCCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgcta                            <  2:4153451/100‑1 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           >  1:1843836/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           >  1:1906683/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           >  2:3629407/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           >  2:3054214/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           >  1:3491326/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           >  1:1650614/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           >  1:774236/1‑100 (MQ=255)
                                                                      ccAGAAATAGACGCCGCCAACAATCCATATAATGGACATCAAAAACGGCCAGAAGAAAACGAAACTCATCATCAGTTCACCAGAGTGAAAGtacgatacg                          <  2:2263794/100‑1 (MQ=255)
                                                                      ccAGAAATAGACGCCGCCAACAATCCATAAAATGGAGATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                           <  2:1863582/99‑1 (MQ=255)
                                                                       cAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGGAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctacgc                         >  2:3205703/1‑100 (MQ=255)
                                                                       cAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctacgc                         >  2:3413069/1‑100 (MQ=255)
                                                                       cAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctacgc                         >  2:1248609/1‑100 (MQ=255)
                                                                         gAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCAt                       >  1:1911560/1‑100 (MQ=255)
                                                                             tAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCGTAggg                  >  1:1290062/1‑100 (MQ=255)
                                                                             tAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggg                  >  2:159904/1‑100 (MQ=255)
                                                                              aGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggg                 <  2:1060241/100‑1 (MQ=255)
                                                                               gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                >  2:2759947/1‑99 (MQ=255)
                                                                               gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                >  1:482683/1‑99 (MQ=255)
                                                                               gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                >  1:1655185/1‑99 (MQ=255)
                                                                               gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                <  1:3494741/100‑2 (MQ=255)
                                                                                aCGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggggt               <  1:3876145/100‑3 (MQ=255)
                                                                                 cgccgcCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCAAAggggg                >  1:3024911/1‑97 (MQ=255)
                                                                                  gccgccAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggggtat             <  1:596575/100‑5 (MQ=255)
                                                                                   ccgccAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggggtatg            <  1:2772155/100‑6 (MQ=255)
                                                                                    cgccAACAATCCATATAATGGACATAATAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGc            <  1:147787/100‑1 (MQ=255)
                                                                                    cgccAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggggtatg            >  2:2728425/1‑94 (MQ=255)
                                                                                       caacaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa         <  1:240268/100‑1 (MQ=255)
                                                                                       caacaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa         <  2:3182314/100‑1 (MQ=255)
                                                                                       caacaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa         <  1:2309246/100‑1 (MQ=255)
                                                                                        aacaaTCCATATAATGGACATAAAAATCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa         <  2:849837/99‑1 (MQ=255)
                                                                                        aacaaTCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAac       <  1:2018979/100‑1 (MQ=255)
                                                                                         ncaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAac       >  1:2073701/2‑100 (MQ=255)
                                                                                         acaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAac       >  2:3664137/1‑100 (MQ=255)
                                                                                         acaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAac       <  1:2428966/100‑1 (MQ=255)
                                                                                           aaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAacac     >  2:2854682/1‑100 (MQ=255)
                                                                                           aaTCCACATAATTTACAGAAAAAACGGCTAGAAGAAAACTAACCTCATCATCAGTTCACCAGATTGAAAGTACGCTACGCATAGGGGTATGCATAACAcc    <  2:3434134/100‑1 (MQ=255)
                                                                                           aaTCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACAcc    <  1:3282188/100‑1 (MQ=255)
                                                                                           aaTCAACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACAcc    <  2:404847/100‑1 (MQ=255)
                                                                                              ccATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAg  >  2:1572375/1‑100 (MQ=255)
                                                                                              ccATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAg  >  1:3681888/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TAGACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAG  >  NC_000913/1087618‑1087810

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: