Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,087,715 C→T 46.9% 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,7150CT46.9% 38.7 / 167.5 132M44I (ATG→ATApgaCbiofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase
Reads supporting (aligned to +/- strand):  ref base C (39/31);  new base T (36/26);  total (75/57)
Fisher's exact test for biased strand distribution p-value = 8.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

TAGACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGAC  >  NC_000913/1087618‑1087812
                                                                                                 |                                                                                                  
 aGACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCAtata                                                                                                 <  1:341796/100‑1 (MQ=255)
  gACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATaa                                                                                                >  2:1962454/1‑100 (MQ=255)
     ggATTATCTTTCAACTGGGGAGCGGGTGCGTTTGCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATg                                                                                              <  1:1063611/99‑1 (MQ=255)
      gATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGa                                                                                            <  2:2015793/100‑1 (MQ=255)
      gATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGa                                                                                            <  1:1111482/100‑1 (MQ=255)
       aTTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGAc                                                                                           <  2:265941/100‑1 (MQ=255)
        ttATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGAc                                                                                           >  1:172631/1‑99 (MQ=255)
        ttATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCAATCCCAGTGGCGTTCATGATAGACCCAGAACTAGACGCCGCCGACAATCCATATAATGGGCa                                                                                          >  1:1033243/1‑100 (MQ=255)
         tATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACAt                                                                                         >  2:379039/1‑100 (MQ=255)
            cTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATaaa                                                                                      >  2:577904/1‑100 (MQ=255)
            cTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATaaa                                                                                      >  2:1382281/1‑100 (MQ=255)
               tCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGACCAGTGGCATTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATaaaaaa                                                                                   <  2:2768961/100‑1 (MQ=255)
                cAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAAc                                                                                  >  2:3094714/1‑100 (MQ=255)
                 aaCTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACg                                                                                 <  1:1600132/100‑1 (MQ=255)
                  aCTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACgg                                                                                <  2:2955135/100‑1 (MQ=255)
                   cTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGc                                                                               >  1:2865661/1‑100 (MQ=255)
                    tGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGcc                                                                              <  2:952884/100‑1 (MQ=255)
                    tGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGcc                                                                              <  2:190968/100‑1 (MQ=255)
                     ggggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCa                                                                             <  2:2669960/100‑1 (MQ=255)
                     ggggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCa                                                                             >  1:3291530/1‑100 (MQ=255)
                     ggggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCa                                                                             <  1:4265299/100‑1 (MQ=255)
                     ggggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCAACATACTGGACATAAAAAACGGcc                                                                              >  1:1892400/1‑99 (MQ=255)
                      gggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCag                                                                            >  1:3835419/1‑100 (MQ=255)
                      gggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCag                                                                            >  1:1470526/1‑100 (MQ=255)
                      gggAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCag                                                                            <  1:1167078/100‑1 (MQ=255)
                      gggAGCTGGTGCGTTTGCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACgg                                                                                <  1:2447714/96‑1 (MQ=255)
                         aGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCagaag                                                                         <  2:1470530/100‑1 (MQ=255)
                         aGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCagaag                                                                         >  2:1398506/1‑100 (MQ=255)
                          gCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCagaaga                                                                        <  2:4135838/100‑1 (MQ=255)
                          gCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCagaaga                                                                        <  2:4136043/100‑1 (MQ=255)
                            tGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGaaa                                                                      >  2:4044949/1‑100 (MQ=255)
                            tGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGaaa                                                                      >  2:2534747/1‑100 (MQ=255)
                            tGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGaaa                                                                      >  2:505574/1‑100 (MQ=255)
                             ggTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGaaaa                                                                     >  2:2314034/1‑100 (MQ=255)
                             ggTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGaaaa                                                                     >  2:31639/1‑100 (MQ=255)
                              gTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAAc                                                                    <  1:294453/100‑1 (MQ=255)
                              gTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAAc                                                                    >  1:47861/1‑100 (MQ=255)
                              gTGCGTTTTCTCCCAACGGCCAGTAGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAAACCACATAATGGACATAAAAAACGGCCAGAAGAAAAc                                                                    <  1:2430043/100‑1 (MQ=255)
                              gTGCGTTTTCCCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAAc                                                                    <  2:753907/100‑1 (MQ=255)
                               tGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACACAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACg                                                                   >  1:3534096/1‑100 (MQ=255)
                               tGCGTATGCTCCCCACGGCCAGTGGCGTTCACGATAGAACCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACg                                                                   <  1:2442358/100‑1 (MQ=255)
                                gCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGa                                                                  >  1:647575/1‑100 (MQ=255)
                                gCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGa                                                                  >  1:1511162/1‑100 (MQ=255)
                                 cGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGaa                                                                 <  2:2391881/100‑1 (MQ=255)
                                 cGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGaa                                                                 <  2:393310/100‑1 (MQ=255)
                                  gTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAAc                                                                >  2:346422/1‑100 (MQ=255)
                                  gTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAAc                                                                >  1:2006327/1‑100 (MQ=255)
                                     ttCTCCCCACGGCCAGTGGCGTTCGCGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtc                                                             <  1:495696/100‑1 (MQ=255)
                                       ctcCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcat                                                           <  1:1975948/100‑1 (MQ=255)
                                         ccccACGGCCAGTGGCGTTCACGATAGACCCAGAAATGGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatca                                                         >  1:2269847/1‑100 (MQ=255)
                                         ccccACGGCCAGTGGCGTTCACGATAGACCCAGAAATCGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatca                                                         >  1:2276494/1‑100 (MQ=255)
                                         ccccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatca                                                         >  1:2627958/1‑100 (MQ=255)
                                          cccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatcat                                                        >  2:79730/1‑100 (MQ=255)
                                          cccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatcat                                                        >  2:3795528/1‑100 (MQ=255)
                                          cccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatcat                                                        >  2:2437540/1‑100 (MQ=255)
                                           ccACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCtcatcatc                                                       <  2:2079843/100‑1 (MQ=255)
                                                 ccAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCa                                                 <  1:1684113/100‑1 (MQ=255)
                                                  cAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCAc                                                >  1:1144368/1‑100 (MQ=255)
                                                  cAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCAc                                                <  2:875523/100‑1 (MQ=255)
                                                   aGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCAcc                                               <  2:269755/100‑1 (MQ=255)
                                                   aGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCAcc                                               >  1:320773/1‑100 (MQ=255)
                                                   aGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCAcc                                               <  1:344443/100‑1 (MQ=255)
                                                    gTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCa                                              >  1:1838333/1‑100 (MQ=255)
                                                          tgcACGATAGACCTAGAAATAGACGCCGCCAGCAATCCACATGATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGa                                        >  2:2887993/3‑100 (MQ=255)
                                                            cACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGaaa                                      <  1:2174089/100‑1 (MQ=255)
                                                            cACGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGaaa                                      >  1:4204222/1‑100 (MQ=255)
                                                            cACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGaaa                                      >  2:823115/1‑100 (MQ=255)
                                                             aCGATAGACCCAGAAATAGACGCCGCCAACAATCCACAAAAAAAAAAAAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAg                                     <  1:31639/100‑1 (MQ=255)
                                                              cGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGt                                    >  1:2390001/1‑100 (MQ=255)
                                                              cGATAGACCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGt                                    <  2:1205166/100‑1 (MQ=255)
                                                               gATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGta                                   >  2:3555174/1‑100 (MQ=255)
                                                                aTAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCTGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtac                                  >  1:2730383/1‑100 (MQ=255)
                                                                 tAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacg                                 <  2:2314931/100‑1 (MQ=255)
                                                                   gACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgct                               <  2:2964186/100‑1 (MQ=255)
                                                                   gACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgct                               >  2:3347832/1‑100 (MQ=255)
                                                                   gACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGAGAAAGtacgct                               >  1:2867851/1‑100 (MQ=255)
                                                                    aCCCAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgcta                              >  1:1330380/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                             >  1:4178234/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                             >  1:2587863/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCATATAATGGACATAAAAAACCGCCAGAAGAAACCGAAACTCATCATCAGTTCACCAGCGTGAACGtac                                  >  1:1746095/1‑95 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                             >  1:388956/1‑100 (MQ=255)
                                                                     cccAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctac                             <  2:1912395/100‑1 (MQ=255)
                                                                        aaaaaTAGACGCCGCCAACAATCCACATAATGAACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCa                          <  1:357906/98‑1 (MQ=255)
                                                                        aGAAATAGACGCCGCCAACAATCCACATAATGGACATATAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGtacgctacgc                           <  1:3959584/99‑1 (MQ=255)
                                                                        aGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCa                          <  1:1398502/100‑1 (MQ=255)
                                                                           aaTAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAg                      >  2:3705383/1‑100 (MQ=255)
                                                                             tAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggg                    <  1:917014/100‑1 (MQ=255)
                                                                             tAGACGCCGCAAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggg                    <  2:294228/100‑1 (MQ=255)
                                                                              aGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggg                   >  2:1613150/1‑100 (MQ=255)
                                                                               gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                  >  2:621582/1‑99 (MQ=255)
                                                                               gACGCCGCCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAggggg                  <  2:2813719/100‑2 (MQ=255)
                                                                               gACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGt                  >  1:1874086/1‑100 (MQ=255)
                                                                               gACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGt                  <  2:1421346/100‑1 (MQ=255)
                                                                                 cgccgcCAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAgggggta                >  2:1747128/1‑97 (MQ=255)
                                                                                    cgccAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCa             >  2:2338504/1‑100 (MQ=255)
                                                                                     gccAACAATCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCa             <  2:3914014/100‑1 (MQ=255)
                                                                                       caacaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATa           >  2:327646/1‑100 (MQ=255)
                                                                                       caacaaTCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATaa          >  1:4227160/1‑100 (MQ=255)
                                                                                         acaaTCCATATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGGTATGCATAac         >  1:3426382/1‑100 (MQ=255)
                                                                                         acaaTCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAaca        >  1:1455610/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
TAGACGGATTATCTTTCAACTGGGGAGCTGGTGCGTTTTCTCCCCACGGCCAGTGGCGTTCACGATAGACCCAGAAATAGACGCCGCCAACAATCCACATAATGGACATAAAAAACGGCCAGAAGAAAACGAACCTCATCATCAGTTCACCAGAGTGAAAGTACGCTACGCATAGGGGTATGCATAACACCAGAC  >  NC_000913/1087618‑1087812

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: