Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,238,726 A→G 17.7% I365V (ATT→GTT)  yjbH → DUF940 family extracellular polysaccharide protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,238,7260AG17.7% 218.7 / 38.7 103I365V (ATT→GTT) yjbHDUF940 family extracellular polysaccharide protein
Reads supporting (aligned to +/- strand):  ref base A (43/41);  new base G (12/6);  total (55/48)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 106 total aligned reads displayed.

ACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGA  >  NC_000913/4238627‑4238825
                                                                                                   |                                                                                                   
aCCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCa                                                                                                     >  1:2812496/1‑100 (MQ=255)
 ccGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCAt                                                                                                    >  1:1074440/1‑100 (MQ=255)
 ccGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCAt                                                                                                    >  2:2254190/1‑100 (MQ=255)
   ggCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTa                                                                                                  >  1:1663129/1‑100 (MQ=255)
   ggCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAAAGATCCGCATTa                                                                                                  >  1:1162346/1‑100 (MQ=255)
     cGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACg                                                                                                <  2:663827/100‑1 (MQ=255)
       aGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGa                                                                                              >  1:2835305/1‑100 (MQ=255)
       aGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGa                                                                                              >  1:1374090/1‑100 (MQ=255)
        gCAGGTGAAATATCGTGATTCGCGCGAAGGCATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGaa                                                                                             <  2:2477823/100‑1 (MQ=255)
           ggTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAAt                                                                                          >  2:1318613/1‑100 (MQ=255)
           ggTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAAt                                                                                          <  1:468365/100‑1 (MQ=255)
           ggTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAAt                                                                                          >  1:732547/1‑100 (MQ=255)
           ggTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAAt                                                                                          >  2:1370300/1‑100 (MQ=255)
           ggTGAAATATCGTGATTCGCGCGAAGGAATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGaaggc                                                                                          >  2:1133341/1‑97 (MQ=255)
             tGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCg                                                                                        <  1:1518483/100‑1 (MQ=255)
              gAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGc                                                                                       <  2:2472492/100‑1 (MQ=255)
               aaaTATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGGAAATCGcc                                                                                      <  1:2383707/100‑1 (MQ=255)
                aaTATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCt                                                                                     >  2:646396/1‑100 (MQ=255)
                 atatCGTGATTCGCGCGAAGGGATCATCCTTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCt                                                                                     <  1:2537407/99‑1 (MQ=255)
                   atCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTaa                                                                                  <  1:1388351/100‑1 (MQ=255)
                   atCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGAATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTaa                                                                                  >  2:1557659/1‑100 (MQ=255)
                    tCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAAc                                                                                 <  2:1895623/100‑1 (MQ=255)
                    tCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAAc                                                                                 <  1:1173689/100‑1 (MQ=255)
                      gTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACAt                                                                               <  2:446590/100‑1 (MQ=255)
                      gTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACAt                                                                               >  1:212373/1‑100 (MQ=255)
                       tGATTCGCGCGAAGGCATCATCCGTGCCAATCCGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGTAAAATCGCCTTAACATg                                                                              <  2:2467751/100‑1 (MQ=255)
                         aTTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATgac                                                                            >  2:1833915/1‑98 (MQ=255)
                         aTTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGcc                                                                            >  1:760653/1‑100 (MQ=255)
                         aTTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGcc                                                                            <  2:212371/100‑1 (MQ=255)
                         aTTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGcc                                                                            <  1:659857/100‑1 (MQ=255)
                           tCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGc                                                                          <  1:646400/100‑1 (MQ=255)
                            cgcgcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGCTCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCa                                                                         >  1:91095/1‑100 (MQ=255)
                              cgcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAgg                                                                       <  2:851161/100‑1 (MQ=255)
                               gcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAgg                                                                       >  2:2080753/1‑99 (MQ=255)
                               gcgAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGt                                                                      >  2:1258321/1‑100 (MQ=255)
                                 gAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTAAGAAAAATCGCCTTAAAAAGCCgaaagtaa                                                                    >  1:1337590/1‑93 (MQ=255)
                                    gggATCATCCGTGCGAATCGTATCGTGATGAACGATTTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTgacga                                                                 <  1:495424/100‑1 (MQ=255)
                                    gggATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTgacga                                                                 >  1:2161036/1‑100 (MQ=255)
                                    gggATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTgacga                                                                 >  1:2079969/1‑100 (MQ=255)
                                     ggATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTgacgac                                                                <  1:2860250/100‑1 (MQ=255)
                                      gATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGaagacg                                                               >  1:3009837/1‑100 (MQ=255)
                                       atcatcCGTGCGAATCGGATCGTGATGAACGATCTGCCTGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACgg                                                              <  1:1611299/100‑1 (MQ=255)
                                       atcatcCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACgg                                                              <  1:1610624/100‑1 (MQ=255)
                                        tcatcCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGa                                                             >  2:2470939/1‑100 (MQ=255)
                                         catcCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGaa                                                            >  1:3010913/1‑100 (MQ=255)
                                           tcCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAAc                                                          >  1:419519/1‑100 (MQ=255)
                                            ccGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGACGCAGGTGACGACGGAAAcc                                                         <  2:1854454/100‑1 (MQ=255)
                                              gTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGa                                                       <  1:2101728/100‑1 (MQ=255)
                                              gTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGa                                                       <  2:1695541/100‑1 (MQ=255)
                                                gCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATg                                                     <  2:2383619/100‑1 (MQ=255)
                                                 cGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGt                                                    >  2:2385417/1‑100 (MQ=255)
                                                   aaTCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGGATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCg                                                  <  2:250138/100‑1 (MQ=255)
                                                    aTCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGc                                                 >  1:2220786/1‑100 (MQ=255)
                                                       ggATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAg                                              >  1:2930059/1‑100 (MQ=255)
                                                       ggATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGAAGGAAACCGATGTCGCCa                                               >  1:136774/1‑99 (MQ=255)
                                                        gATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGc                                             >  1:634853/1‑100 (MQ=255)
                                                         aTCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGcc                                            <  2:405576/100‑1 (MQ=255)
                                                            gTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCAGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGa                                         >  1:1370625/1‑100 (MQ=255)
                                                              gatgaACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGaaa                                       >  1:2375298/1‑100 (MQ=255)
                                                               atgaACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGACAATCGCCTTAACATGCAGCTGGTGAAGACGGAAACCGATGTCGCCAGCCTCaaaa                                      >  2:2546260/1‑100 (MQ=255)
                                                                tgaACGATCTGCCGGATGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGaaaaa                                     >  1:1348745/1‑100 (MQ=255)
                                                                tgaACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGaaaaa                                     <  2:679420/100‑1 (MQ=255)
                                                                   ctgATCTGCCGGATGGGATCAACACGATCCGCTTTACGGAAAATCGCCTTAACATGTCGCAGGTGACGTCGGAAACCGATGTCGCCAGCCTGAAAAatcn                                  <  2:246851/98‑2 (MQ=255)
                                                                   aCGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatca                                  <  1:646219/100‑1 (MQ=255)
                                                                     naTCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatcaag                                >  1:246854/2‑98 (MQ=255)
                                                                     gATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatcatc                                >  1:1866197/1‑100 (MQ=255)
                                                                     gATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatcatc                                <  1:265677/100‑1 (MQ=255)
                                                                     gATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAatcatc                                <  1:488033/100‑1 (MQ=255)
                                                                       tCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCAtctc                              <  2:157462/100‑1 (MQ=255)
                                                                        cTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCg                             >  1:1671966/1‑100 (MQ=255)
                                                                         tGCCGGATGGGATCAAAACGATACGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCg                             >  2:1103505/1‑99 (MQ=255)
                                                                         tGCCGGATGGGATAAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATACtctc                              >  2:2309508/1‑98 (MQ=255)
                                                                          gCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGcc                           >  1:2513376/1‑100 (MQ=255)
                                                                          gCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGcc                           >  1:1472144/1‑100 (MQ=255)
                                                                           ccGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCg                          <  1:2801561/100‑1 (MQ=255)
                                                                               aTGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGgaga                      <  2:1663109/100‑1 (MQ=255)
                                                                               aTGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCATGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCg                          <  2:184968/96‑1 (MQ=255)
                                                                                tGGGATCAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGagag                     <  1:2521010/100‑1 (MQ=255)
                                                                                   gATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCg                  >  2:311326/1‑100 (MQ=255)
                                                                                   gATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCg                  >  2:1769137/1‑100 (MQ=255)
                                                                                   gATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGACAGCCTGAAAAATAATCTCGCCGGAGAGCCg                  >  1:94366/1‑100 (MQ=255)
                                                                                    aTCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGt                 >  2:1318742/1‑100 (MQ=255)
                                                                                     tCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGtt                <  1:833045/100‑1 (MQ=255)
                                                                                      cAAAACGATCCGCGTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTg               >  1:1883885/1‑100 (MQ=255)
                                                                                      cAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTg               <  1:2682614/100‑1 (MQ=255)
                                                                                      cAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTg               <  2:2952174/100‑1 (MQ=255)
                                                                                       aaaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTgg              >  1:698179/1‑100 (MQ=255)
                                                                                       aaaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTgg              <  1:2969473/100‑1 (MQ=255)
                                                                                       aaaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTg               >  1:222531/1‑99 (MQ=255)
                                                                                        aaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTggg             <  2:1543708/100‑1 (MQ=255)
                                                                                        aaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTggg             <  1:1667893/100‑1 (MQ=255)
                                                                                        aaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTggg             <  1:1470107/100‑1 (MQ=255)
                                                                                        aaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTggg             >  2:2163430/1‑100 (MQ=255)
                                                                                        aaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATATCGCCAGCCTGAAAAATCATCTCTACAGAAAGCCTTTgg              >  1:2166001/1‑99 (MQ=255)
                                                                                        aaaCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCCATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTggg             <  1:1745706/100‑1 (MQ=255)
                                                                                          aCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGcc           <  1:1947513/100‑1 (MQ=255)
                                                                                          aCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGcc           <  1:1360564/100‑1 (MQ=255)
                                                                                          aCGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGcc           <  1:850935/100‑1 (MQ=255)
                                                                                                   aTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAacga  <  1:2716560/100‑1 (MQ=255)
                                                                                                   aTTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAacg   <  2:336581/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
ACCGGCGAGCAGGTGAAATATCGTGATTCGCGCGAAGGGATCATCCGTGCGAATCGGATCGTGATGAACGATCTGCCGGATGGGATCAAAACGATCCGCATTACGGAAAATCGCCTTAACATGCCGCAGGTGACGACGGAAACCGATGTCGCCAGCCTGAAAAATCATCTCGCCGGAGAGCCGTTGGGCCACGAAACGA  >  NC_000913/4238627‑4238825

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