Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,087,795:1 +G 100% coding (52/1326 nt) 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,7911.G76.7% 393.2 / 134.1 150coding (56/1326 nt)pgaCbiofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase
Reads supporting (aligned to +/- strand):  ref base . (18/17);  new base G (64/50);  total (83/67)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 163 total aligned reads displayed.

TAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATAAGGTC  >  NC_000913/1087695‑1087889
                                                                                                 |                                                                                                  
tAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCGTAggg                                                                                                 >  1:1290062/1‑100 (MQ=255)
  gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                                                                                               >  1:715298/1‑99 (MQ=255)
  gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                                                                                               <  1:3494741/100‑2 (MQ=255)
  gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                                                                                               >  1:482683/1‑99 (MQ=255)
  gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                                                                                               >  2:2759947/1‑99 (MQ=255)
    cgccgcCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATACGGGGTa                                                                                              >  2:1588145/1‑100 (MQ=255)
      ccgccAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggggtatg                                                                                           <  1:2772155/100‑6 (MQ=255)
       cgccAACAATCCATATAATGGACATAATAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGc                                                                                           <  1:147787/100‑1 (MQ=255)
       cgccAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggggtatg                                                                                           >  2:2728425/1‑94 (MQ=255)
        gccAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCa                                                                                          <  1:1248615/100‑1 (MQ=255)
          caacaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa                                                                                        <  2:3182314/100‑1 (MQ=255)
          caacaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa                                                                                        <  1:2309246/100‑1 (MQ=255)
          caacaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa                                                                                        <  1:3866605/100‑1 (MQ=255)
           aacaaTCCATATAATGGACATAAAAATCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa                                                                                        <  2:849837/99‑1 (MQ=255)
            ncaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAac                                                                                      >  1:2073701/2‑100 (MQ=255)
            acaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAac                                                                                      <  1:2428966/100‑1 (MQ=255)
              aaTCCACATAATTTACAGAAAAAACGGCTAGAAGAAAACTAACCTCATCATCAGTTCACCAGATTGAAAGTACGCTACGCATAGGGGTATGCATAACAcc                                                                                   <  2:3434134/100‑1 (MQ=255)
              aaTCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACAcc                                                                                   <  1:3282188/100‑1 (MQ=255)
              aaTCAACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACAcc                                                                                   <  2:404847/100‑1 (MQ=255)
               aTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACAcc                                                                                   >  1:3296643/1‑100 (MQ=255)
                 ccATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAg                                                                                 >  1:3681888/1‑100 (MQ=255)
                 ccATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAg                                                                                 >  2:1572375/1‑100 (MQ=255)
                  cataTAATGTACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGGAAGTACGCTATGTCTAGGGGGTCTGCATAaca                                                                                     >  2:1849982/4‑95 (MQ=255)
                  cataTAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGAATGCATAACAACAGa                                                                                >  2:3593453/4‑100 (MQ=255)
                  cataTAATGGACATAAAAAACGGCCAGAAGAAAACGAACATCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGa                                                                                >  2:1433469/4‑100 (MQ=255)
                   ataTAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGAc                                                                               >  1:3746857/3‑100 (MQ=255)
                     aTAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACan                                                                             <  2:2073690/100‑2 (MQ=255)
                      tAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACata                                                                            >  2:333781/1‑100 (MQ=255)
                      tAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACata                                                                            >  1:2904959/1‑100 (MQ=255)
                      tAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACata                                                                            >  1:2504845/1‑100 (MQ=255)
                      tAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACataa                                                                           >  2:711998/1‑100 (MQ=255)
                       aaTGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACataa                                                                           <  2:774232/100‑1 (MQ=255)
                        aTGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACataat                                                                          >  2:2998460/1‑100 (MQ=255)
                        aTGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACataat                                                                          <  1:2035352/100‑1 (MQ=255)
                           gACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAAtata                                                                       >  1:1545984/1‑100 (MQ=255)
                                aaaaaaCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATaaaaaa                                                                  >  1:2554188/1‑100 (MQ=255)
                                aaaaaaCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAac                                                                 >  1:1983135/1‑100 (MQ=255)
                                 aaaaaCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAac                                                                 >  2:2586726/1‑100 (MQ=255)
                                 aaaaaCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAac                                                                 <  1:3719365/100‑1 (MQ=255)
                                 aaaaaCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAac                                                                 <  1:2117267/100‑1 (MQ=255)
                                  aaaaCGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAacga                                                               <  1:3191515/100‑1 (MQ=255)
                                     aaaaaaaaaagaaAACGAACCTAATCATCAGTTCACCAGAGTAAAAGTACGATACGCATAGGGGGTATTCATAACACAAAAAAAAAAAAAAAAAacgata                                                             <  2:2592386/92‑1 (MQ=255)
                                      cGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAacgatac                                                            <  1:3818628/100‑1 (MQ=255)
                                        gCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAacgatacga                                                          >  1:4114326/1‑100 (MQ=255)
                                        gCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAacgatacga                                                          >  2:1156372/1‑100 (MQ=255)
                                         ccAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAacgatacgat                                                         >  1:3049658/1‑100 (MQ=255)
                                           agaagaAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGc                                                       <  2:4038179/100‑1 (MQ=255)
                                             aagaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGa                                                     >  2:51261/1‑100 (MQ=255)
                                             aagaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGAt                                                    <  1:712002/100‑1 (MQ=255)
                                              agaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATCATATAAAAAACGATACGATGCgct                                                    >  2:3794615/1‑98 (MQ=255)
                                              agaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGAt                                                    >  1:3106408/1‑100 (MQ=255)
                                               gaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGAtt                                                   >  1:3993515/1‑100 (MQ=255)
                                               gaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGAtt                                                   <  2:1253610/100‑1 (MQ=255)
                                               gaaAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGAtt                                                   >  1:738647/1‑100 (MQ=255)
                                                aaaaCGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTa                                                  >  2:1179148/1‑100 (MQ=255)
                                                aaaaCGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTa                                                  <  1:340284/100‑1 (MQ=255)
                                                 aaaCGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTaa                                                 >  2:2797367/1‑100 (MQ=255)
                                                 aaaCGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAAt                                                <  2:1728816/100‑1 (MQ=255)
                                                 aaaCGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAAt                                                >  2:2951523/1‑100 (MQ=255)
                                                   aCGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATc                                               >  2:1143133/1‑100 (MQ=255)
                                                        ccTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCAtttt                                          >  2:2024221/1‑99 (MQ=255)
                                                          tcatcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAGTCATTTTTCg                                       <  1:3453826/100‑1 (MQ=255)
                                                           catcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCg                                       <  1:868162/100‑1 (MQ=255)
                                                            atcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCgg                                      >  1:4107990/1‑100 (MQ=255)
                                                            atcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCgg                                      <  1:2980279/100‑1 (MQ=255)
                                                            atcatcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCgg                                      >  1:2304274/1‑100 (MQ=255)
                                                            atcatcaGTTCACCAGAGTGAAACTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCg                                       <  2:1822855/99‑1 (MQ=255)
                                                              catcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGAt                                    >  1:3993364/1‑100 (MQ=255)
                                                               atcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAGCACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATa                                   >  1:2949365/1‑100 (MQ=255)
                                                               atcaGTTCACCAGAATGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATa                                   <  2:3350140/100‑1 (MQ=255)
                                                                tcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATAc                                  >  2:2881183/1‑100 (MQ=255)
                                                                tcaGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATAcc                                 >  1:1161917/1‑100 (MQ=255)
                                                                 caGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATAcc                                 <  1:1988481/100‑1 (MQ=255)
                                                                  aGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCa                                >  1:4091272/1‑100 (MQ=255)
                                                                  aGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAg                               >  1:4141021/1‑100 (MQ=255)
                                                                   gTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAg                               <  1:2967607/100‑1 (MQ=255)
                                                                     tCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGc                             >  1:1095890/1‑100 (MQ=255)
                                                                     tCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGc                             <  1:1826775/100‑1 (MQ=255)
                                                                       aCCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTg                           >  1:540660/1‑100 (MQ=255)
                                                                       aCCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTg                           <  1:975006/100‑1 (MQ=255)
                                                                        ccAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCt                            >  2:3833111/1‑98 (MQ=255)
                                                                         cAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGtt                         <  1:238148/100‑1 (MQ=255)
                                                                            agTGAAAGTACGCTACGAATAGGGGGTATGCATAACTTCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTg                        <  2:1040266/98‑1 (MQ=255)
                                                                              tGAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGaaaa                    <  1:2362314/100‑1 (MQ=255)
                                                                               gAAAGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAAc                   <  1:623999/100‑1 (MQ=255)
                                                                                 aaGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTc                 <  2:2030053/100‑1 (MQ=255)
                                                                                 aaGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCa                <  2:4127085/100‑1 (MQ=255)
                                                                                 aaGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCa                <  1:1239006/100‑1 (MQ=255)
                                                                                  aGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCa                <  1:2199597/100‑1 (MQ=255)
                                                                                  aGTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCa                <  2:110256/100‑1 (MQ=255)
                                                                                  aGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAg               <  1:1489764/100‑1 (MQ=255)
                                                                                   gTACGCTACGCATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAg               <  2:910154/100‑1 (MQ=255)
                                                                                     acgctacgcATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGa             >  2:4143131/1‑100 (MQ=255)
                                                                                      cgctacgcATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGAc            <  2:3047180/100‑1 (MQ=255)
                                                                                       gctacgcATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCTGATACCAGTCTGTTGAAGACTCAGGACg           >  2:286913/1‑100 (MQ=255)
                                                                                           cgcATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATa       <  2:2763173/100‑1 (MQ=255)
                                                                                           cgcATAGGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATa       >  2:3184484/1‑100 (MQ=255)
                                                                                           cgcATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATaa      <  1:2850391/100‑1 (MQ=255)
                                                                                             cATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATAAgg    >  1:3691720/1‑100 (MQ=255)
                                                                                               tAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATAAGGTc  >  1:1834942/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
TAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGACATAATATAAAAAACGATACGATGCGATTAATCATTTTTCGGATACCAGGCTGTTGAAAACTCAGGACGAATAAGGTC  >  NC_000913/1087695‑1087889

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: