Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 17.5% 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,9860GT17.5% 387.0 / 72.3 149intergenic (+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 (69/54);  new base T (16/10);  total (85/64)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 164 total aligned reads displayed.

AAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAT  >  NC_000913/4049887‑4050084
                                                                                                   |                                                                                                  
aaaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                    >  2:1218542/1‑100 (MQ=255)
aaaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                    >  1:2399645/1‑100 (MQ=255)
aaaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                    >  2:489981/1‑100 (MQ=255)
 aaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagt                                                                                                   >  2:2280674/1‑100 (MQ=255)
 aaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagt                                                                                                   <  2:379913/100‑1 (MQ=255)
       gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                             >  1:221835/1‑100 (MQ=255)
       gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                             >  1:1121731/1‑100 (MQ=255)
         cgcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTg                                                                                           <  2:1757457/100‑1 (MQ=255)
           cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTggg                                                                                         >  2:1324694/1‑100 (MQ=255)
             tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                       <  1:3229150/100‑1 (MQ=255)
             tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                       <  2:2958359/100‑1 (MQ=255)
              aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCg                                                                                      <  1:2358876/100‑1 (MQ=255)
              aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTAATGACTGGGGCg                                                                                      <  1:1781400/100‑1 (MQ=255)
              aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCg                                                                                      <  2:361837/100‑1 (MQ=255)
              aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCg                                                                                      >  1:1946085/1‑100 (MQ=255)
               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCg                                                                                      <  1:3093174/99‑1 (MQ=255)
               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGt                                                                                     >  1:1241904/1‑100 (MQ=255)
               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCAAGTAATCGGATTTGGCTGAATATTTTAGCCGCCGCAGTCAGTAATGACTGGGGCGt                                                                                     >  2:1389159/1‑100 (MQ=255)
                ggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCggc                                                                                    <  2:1907963/100‑3 (MQ=255)
                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                <  2:2879919/100‑1 (MQ=255)
                     agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTATCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTa                                                                               <  1:988535/100‑1 (MQ=255)
                     agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGCggctaaaa                                                                               <  2:3094352/100‑8 (MQ=38)
                     agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTa                                                                               <  1:1060204/100‑1 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                              <  1:1435916/100‑1 (MQ=255)
                        ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                            <  2:2712830/100‑1 (MQ=255)
                        ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                            >  2:469457/1‑100 (MQ=255)
                          tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                          >  1:224053/1‑100 (MQ=255)
                          tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                          <  2:1857210/100‑1 (MQ=255)
                          tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                          >  2:488066/1‑100 (MQ=255)
                           aaGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc                                                                         <  1:1825343/100‑1 (MQ=255)
                           aaGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc                                                                         <  2:353293/100‑1 (MQ=255)
                            aGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCg                                                                        <  2:774327/100‑1 (MQ=255)
                              aTGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaa                                                                      >  1:701186/1‑100 (MQ=255)
                              aTGTTCTATCTTTCAGAACTTTTACTTCACGTAATAGGATTTGGCTGAATATTTTAGTCGCCCCAGTCAGTAATGACTGGGGCGTCTTTTATTAGGc                                                                         >  1:2844235/1‑97 (MQ=255)
                               tGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTGTATTGGGCgaaa                                                                     >  1:1707012/1‑100 (MQ=255)
                                gTTCTATCTTTCAGACCTTTTACTTCACGTAATCGTATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGTGCgaaag                                                                    <  2:1630827/100‑1 (MQ=255)
                                gTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaag                                                                    <  2:736621/100‑1 (MQ=255)
                                  tCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaa                                                                  >  2:881168/1‑100 (MQ=255)
                                    tATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaa                                                                >  2:1133056/1‑100 (MQ=255)
                                     aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaag                                                               <  2:850942/100‑1 (MQ=255)
                                       cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                             <  1:2639467/100‑1 (MQ=255)
                                       cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCACCAGTCAGTAATGAATGGGAAGTTTTTTCTTGGGAGcaagacaaga                                                              >  2:281169/1‑99 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaaga                                                              >  2:1970601/1‑98 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                            >  1:433424/1‑100 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                            >  2:2434405/1‑100 (MQ=255)
                                         ttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                            >  1:2679767/1‑99 (MQ=255)
                                          tCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCg                                                          <  2:1083565/100‑1 (MQ=255)
                                            aGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGACAGAAAAGATCCGTa                                                        >  1:370383/1‑100 (MQ=255)
                                             gACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGAGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                                       >  2:1385737/1‑100 (MQ=255)
                                                      cTTCGCGTAATCGGATTTGGCTGCATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCTAAAGAAAAGATCCGTAATGCCTGATg                                              <  1:1043768/100‑1 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                              <  2:1638385/100‑1 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                              >  2:2416180/1‑100 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                              >  1:503225/1‑100 (MQ=255)
                                                       ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgc                                             <  2:3176508/100‑1 (MQ=255)
                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                            <  2:2327334/100‑1 (MQ=255)
                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                            >  1:1493420/1‑100 (MQ=255)
                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                            <  2:1121734/100‑1 (MQ=255)
                                                          aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                          >  2:2804506/1‑100 (MQ=255)
                                                               aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGtt                                     <  2:2896269/100‑1 (MQ=255)
                                                                tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTTTTTTTTTTTAGAAAAGATCCGTAATGCCTGATGCGCTATGttt                                    >  1:1384194/1‑100 (MQ=255)
                                                                tCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTTTTTTTTTGAAAAGATCAGTAACATCTGATGCGCTATGttt                                    >  1:1475540/1‑100 (MQ=255)
                                                                  ggATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTTGTTTTTTTAAAAGATATGTAATGCCTGATATGCTATTTTTAt                                  >  1:273561/1‑100 (MQ=255)
                                                                   gATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                 <  1:1686225/100‑1 (MQ=255)
                                                                   gATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                 <  1:526385/100‑1 (MQ=255)
                                                                    aTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAGTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                <  2:1114647/100‑1 (MQ=255)
                                                                       tGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGc                             <  2:1151056/100‑1 (MQ=255)
                                                                        ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                            <  2:1932521/100‑1 (MQ=255)
                                                                        ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                            >  2:2596578/1‑100 (MQ=255)
                                                                         gCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCa                           <  2:1621565/100‑1 (MQ=255)
                                                                          cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                          >  2:2450670/1‑100 (MQ=255)
                                                                            gAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACg                        <  1:1740045/100‑1 (MQ=255)
                                                                             aaTATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACgg                       >  2:2817767/1‑100 (MQ=255)
                                                                              atatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGt                      <  1:1866063/100‑1 (MQ=255)
                                                                               tatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTa                     <  1:1549693/100‑1 (MQ=255)
                                                                               tatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTa                     >  2:3231031/1‑100 (MQ=255)
                                                                                ataaaaaacgccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAg                    >  2:3089024/9‑100 (MQ=37)
                                                                                ataaaaaacgccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAg                    >  2:3015215/9‑100 (MQ=37)
                                                                                       ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa             <  2:1001958/100‑1 (MQ=255)
                                                                                        cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTTTTTTTAATTTAATATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa            >  1:1789783/1‑100 (MQ=255)
                                                                                        cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa            >  2:2125031/1‑100 (MQ=255)
                                                                                        cgccCCAGTCAGTAACGACTGGGGCGTTTTTTATTTGTCGATAGCAATGTCCCGCGACGCCTGATTCCGTACGATTATCAgcc                             >  1:3337537/1‑81 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGTGCTATGTTTATCAGGCCAACGGTAGAATTGTaa            >  1:2126023/1‑99 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt           >  2:1701110/1‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt           >  2:852135/1‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt           >  2:1882994/1‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt           >  2:2847084/1‑100 (MQ=255)
                                                                                         gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGTGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt           >  1:1839928/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc          >  2:2697828/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc          >  2:1970180/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc          >  2:2994899/1‑100 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc          >  2:710607/1‑100 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTAGGCGAAAGAAAAGATCAGTAATGACTGATGCGCTATGTTTATCAGTACAACGCTAGACTTGTAATc          >  1:1166712/1‑100 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTGTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt           <  1:2246261/99‑1 (MQ=255)
                                                                                           cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTActct         >  1:2026169/1‑100 (MQ=255)
                                                                                            ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa        >  1:2874770/1‑100 (MQ=255)
                                                                                            ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa        >  1:1950814/1‑100 (MQ=255)
                                                                                             cagtcagtAACGATTTGGGCTTTTTTTATTGGGCAAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAATGGTAGAATTGTAATCTAt       <  1:3228454/100‑1 (MQ=255)
                                                                                              agtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt      >  2:267984/1‑100 (MQ=255)
                                                                                                 cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa   <  2:320786/100‑1 (MQ=255)
                                                                                                  agtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt  <  2:255459/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAT  >  NC_000913/4049887‑4050084

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