Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 15.6% intergenic (+234/+147) polA → / ← yihA 5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease/cell division GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,049,9860GT15.6% 470.8 / 69.7 168intergenic (+234/+147)polA/yihA5' to 3' DNA polymerase and 3' to 5'/5' to 3' exonuclease/cell division GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (71/70);  new base T (13/13);  total (85/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

Only 100 of 185 total aligned reads displayed.

AAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATT  >  NC_000913/4049887‑4050085
                                                                                                   |                                                                                                   
 aaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagt                                                                                                    >  1:3555236/1‑100 (MQ=255)
    ttAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAAt                                                                                                 <  2:882691/100‑1 (MQ=255)
      aGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGa                                                                                               <  2:2341370/100‑1 (MQ=255)
        tCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACt                                                                                             <  2:716189/100‑1 (MQ=255)
        tCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACt                                                                                             >  2:2900096/1‑100 (MQ=255)
         cgcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                                                                                            >  2:1287820/1‑100 (MQ=255)
         cgcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTGTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                                                                                            <  2:2762930/100‑1 (MQ=255)
         cgcgTAGGGTACAGAGGTAAGATGACCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                                                                                            >  2:1163862/1‑100 (MQ=255)
          gcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                           >  1:1325082/1‑100 (MQ=255)
           cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                          <  1:3668757/100‑1 (MQ=255)
           cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                          <  2:3139711/100‑1 (MQ=255)
           cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                          <  1:2180752/100‑1 (MQ=255)
            gatctgTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggg                                                                                         <  2:3629669/95‑1 (MQ=255)
            gTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgggg                                                                                         <  1:9408/100‑1 (MQ=255)
            gTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgggg                                                                                         <  1:3309629/100‑1 (MQ=255)
             tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                        <  1:3600171/100‑1 (MQ=255)
              aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCg                                                                                       <  1:2889374/100‑1 (MQ=255)
               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCgg                                                                                      <  2:2140584/100‑2 (MQ=255)
                 gTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTCGGCTGAATATTTTCGCCGCCCCAGTCACTAATGCCTGGGGCGtt                                                                                     >  1:1399103/1‑99 (MQ=255)
                 gTACAGAGGTAAGATGTTCTATCCTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggggagatc                                                                                    >  1:3629697/1‑95 (MQ=255)
                   aCAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGttttt                                                                                  <  2:2599000/100‑1 (MQ=255)
                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                 >  1:1954985/1‑100 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                               >  1:463511/1‑100 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                               >  1:2279853/1‑100 (MQ=255)
                       agGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAtt                                                                              <  2:2400427/100‑1 (MQ=255)
                        ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                             <  2:1971130/100‑1 (MQ=255)
                          tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                           <  2:2666877/100‑1 (MQ=255)
                          tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                           <  1:1663359/100‑1 (MQ=255)
                            aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCg                                                                         <  1:1622732/100‑1 (MQ=255)
                            aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCg                                                                         <  2:515515/100‑1 (MQ=255)
                              aTGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaa                                                                       >  2:1006846/1‑100 (MQ=255)
                               tGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaa                                                                      <  1:3359937/100‑1 (MQ=255)
                                gTTCTATCTTTCGGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaag                                                                     <  1:2466172/100‑1 (MQ=255)
                                 ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaaga                                                                    >  1:28888/1‑100 (MQ=255)
                                   cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGTCgaaagaaa                                                                  >  2:3066327/1‑100 (MQ=255)
                                    tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                 <  1:1715915/100‑1 (MQ=255)
                                    tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                 <  2:3784180/100‑1 (MQ=255)
                                    tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                 <  2:2992350/100‑1 (MQ=255)
                                     aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTTTACGaaagaaaag                                                                >  2:1306467/1‑100 (MQ=255)
                                       cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACGAGGGCGTTTTTGATTGGGCGGAAGAACAGAt                                                              >  2:43835/1‑100 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                             <  2:3542499/100‑1 (MQ=255)
                                         ttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                                            >  1:74956/1‑100 (MQ=255)
                                         ttCAGACCTTTTACTTCACGTAATCGGAGTTGGCTGAATATTTTAGCCGCCCGAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                             <  2:2469794/99‑1 (MQ=255)
                                         tactcattttttACTTCACGTAATCGGATTTGTCTGAATATTTTTCCCGTTCCAGTCAGTAATGACTTATGCGTTTTTTAATTGCTTaaagaaaan                                                                <  2:2233944/88‑2 (MQ=255)
                                           cAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGt                                                          <  2:1657290/100‑1 (MQ=255)
                                            aGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTa                                                         <  1:1866372/100‑1 (MQ=255)
                                             gACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                                        >  2:509821/1‑100 (MQ=255)
                                               ccTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATg                                                      <  1:1041044/100‑1 (MQ=255)
                                                  tttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCt                                                   <  1:3940373/100‑1 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTTTTTTTCGAAAGAAAAGATCCGTAATGCCTGATg                                               >  2:453792/1‑100 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                               <  1:450141/100‑1 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                               >  1:342651/1‑100 (MQ=255)
                                                       ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgc                                              <  1:1006852/100‑1 (MQ=255)
                                                       ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgc                                              >  1:252139/1‑100 (MQ=255)
                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                             <  2:602569/100‑1 (MQ=255)
                                                         cACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcgc                                            >  1:2562645/1‑100 (MQ=255)
                                                          aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                           <  1:3723473/100‑1 (MQ=255)
                                                          aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                           >  1:3569903/1‑100 (MQ=255)
                                                           cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTa                                          >  1:1297761/1‑100 (MQ=255)
                                                            gTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTAt                                         >  2:48196/1‑100 (MQ=255)
                                                            gTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTAt                                         >  1:2255963/1‑100 (MQ=255)
                                                             tAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                        <  2:2795055/100‑1 (MQ=255)
                                                                 cGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTa                                    <  1:3066346/100‑1 (MQ=255)
                                                                   gATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                  >  1:860632/1‑100 (MQ=255)
                                                                    aTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                 >  2:3684556/1‑100 (MQ=255)
                                                                        ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGTGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcn                             <  2:3450371/100‑2 (MQ=255)
                                                                          cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                           >  1:786058/1‑100 (MQ=255)
                                                                           tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAAc                          >  1:2599502/1‑100 (MQ=255)
                                                                            gAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACg                         >  2:820123/1‑100 (MQ=255)
                                                                            gAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACg                         >  1:1795248/1‑100 (MQ=255)
                                                                              atatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACgg                        <  2:1812693/99‑1 (MQ=255)
                                                                                atTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAg                     >  1:3153673/1‑100 (MQ=255)
                                                                                atTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAg                     >  1:173416/1‑100 (MQ=255)
                                                                                atTTTAGCCGCCCCAGCCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAg                     <  2:2307335/100‑1 (MQ=255)
                                                                                 ttttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGa                    >  2:1199789/1‑100 (MQ=255)
                                                                                 ttttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGa                    >  1:2145619/1‑100 (MQ=255)
                                                                                   ttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                  <  1:1199798/100‑1 (MQ=255)
                                                                                   ttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                  <  1:2227059/100‑1 (MQ=255)
                                                                                    tAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt                 >  1:2905901/1‑100 (MQ=255)
                                                                                     aGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                <  2:2735706/100‑1 (MQ=255)
                                                                                     aGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                <  1:958817/100‑1 (MQ=255)
                                                                                     aGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                <  2:1491314/100‑1 (MQ=255)
                                                                                      gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGt               >  2:943255/1‑100 (MQ=255)
                                                                                      gacgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                >  2:3943599/3‑99 (MQ=255)
                                                                                        cgccCCAGTCATTACTGACTGGGGCGTCTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa             >  2:3679699/1‑100 (MQ=39)
                                                                                        cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                <  2:3868968/97‑1 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTTTTCTTAATAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt            >  2:1809133/1‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt            >  1:2656091/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:1523328/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  2:3052914/1‑100 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:3510549/1‑100 (MQ=255)
                                                                                            ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa         >  1:388387/1‑100 (MQ=255)
                                                                                            ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa         >  1:3713632/1‑100 (MQ=255)
                                                                                             cagtcattACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAt        >  1:1556320/1‑100 (MQ=255)
                                                                                               gtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg      <  2:2828831/100‑1 (MQ=255)
                                                                                                 cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGCATTGTGATccactgca    >  2:2392949/1‑93 (MQ=255)
                                                                                                 cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa    >  1:759157/1‑100 (MQ=255)
                                                                                                   gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa    <  2:1559798/98‑1 (MQ=255)
                                                                                                   gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAtt  >  2:1849883/1‑100 (MQ=255)
                                                                                                   gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAtt  >  1:240898/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATT  >  NC_000913/4049887‑4050085

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: