Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 15.8% 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.8% 365.7 / 52.0 138intergenic (+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 (65/50);  new base T (15/7);  total (81/57)
Fisher's exact test for biased strand distribution p-value = 3.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 149 total aligned reads displayed.

AAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATT  >  NC_000913/4049887‑4050085
                                                                                                   |                                                                                                   
aaaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                     >  1:815101/1‑100 (MQ=255)
   gTTAGTCGCGTATGGTACAGAGGTAAGATGTTCTATCTTTCAGATCTTTTACTTCACGTAATCGGAATCTGCTGAATATTTTAGCCCCCCCATTCAGTaa                                                                                                  >  2:608715/1‑100 (MQ=255)
    ttAGTCTCGTAGGGTTCATAGGTAAGATATTCTATTTTTCAGCCCTGTTACTTCACGTAATCGGATTTGGCTGAATATTTTATCCGCCCcagtcagt                                                                                                    <  2:3534227/97‑1 (MQ=255)
       gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGAc                                                                                              >  1:241954/1‑100 (MQ=255)
       gTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATCTTAGCCGCCCCAGTCAGTAATGAc                                                                                              >  1:437137/1‑100 (MQ=255)
           cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTggg                                                                                          <  1:2408187/100‑1 (MQ=255)
            gTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgggg                                                                                         <  1:135610/100‑1 (MQ=255)
            gTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTgggg                                                                                         <  1:478826/100‑1 (MQ=255)
            gTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggg                                                                                         <  2:1087448/100‑1 (MQ=255)
               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGt                                                                                      >  1:959463/1‑100 (MQ=255)
               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGt                                                                                      >  2:3104/1‑100 (MQ=255)
                    cAGAGGTAAGATGTTCTATTTTCCAGACCTTTTGCTTTACGTAATCGGATTTGGCTGAAGATTTTATCCGCCCCACTCGGTAATGACTGGGGCGcttttt                                                                                 <  2:2046969/100‑1 (MQ=255)
                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                 >  2:5363/1‑100 (MQ=255)
                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                 <  2:1601840/100‑1 (MQ=255)
                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                 >  2:3286223/1‑100 (MQ=255)
                     agagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTa                                                                                <  1:3272313/100‑1 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                               >  1:2322087/1‑100 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                               >  2:1055746/1‑100 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                               >  1:1169720/1‑100 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTCGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACGGGAGCGtttttca                                                                                >  1:1878065/1‑97 (MQ=255)
                        gttagattGTCCTATCTTTCGAACTTTTTCGTTCCCGTAATGGGCTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                             <  1:1691560/92‑1 (MQ=255)
                        ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                             <  2:694530/100‑1 (MQ=255)
                        ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                             <  2:2056025/100‑1 (MQ=255)
                         gTAAGATGTTCTTTCTTTTAGACCTTTTACTTCACGGACTCGGATTTGGCTGAATATTTTAGCCGCCCCCGTCAGTTATGACTGTGGCGTTTTTTATTg                                                                             <  2:1959226/99‑1 (MQ=255)
                         gTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATCTTCGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                            >  2:2554324/1‑100 (MQ=255)
                           aataTGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc                                                                          <  1:2145234/97‑1 (MQ=255)
                            aGATGTAATTTCTTCCAGACCTTTTACTTCAGGTAATTTGATTTCTCTGAATATTTTAGCCGTTCAATTCAATAATGTCTAGGGCGTTTTTTAt                                                                               >  2:3221942/1‑94 (MQ=255)
                             gATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCga                                                                        >  2:1415013/1‑100 (MQ=255)
                             gATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCga                                                                        >  1:1892067/1‑100 (MQ=255)
                                gTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaag                                                                     >  1:780920/1‑100 (MQ=255)
                                 ttCTATCTTTTAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaaga                                                                    <  2:2385967/100‑1 (MQ=255)
                                  tCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaa                                                                   >  1:2793624/1‑100 (MQ=255)
                                  tCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaagaa                                                                   >  1:3021918/1‑100 (MQ=255)
                                     aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaag                                                                <  1:2179900/100‑1 (MQ=255)
                                     aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGaaagaaaag                                                                <  1:1253042/100‑1 (MQ=255)
                                       cTTTCAGCCCTATTACTTCACGTCATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATAACTCGGGCGTTTTTTATTGAGCGAAAGAAAAGAt                                                              >  2:1914376/1‑100 (MQ=255)
                                       cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                              <  2:470272/100‑1 (MQ=255)
                                        tttCAGTCCTTTTACTTCACGTAATCGGATTTGTCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCTTTTTTTATTGGGCGAAAGAAAAGATc                                                             <  2:2579459/100‑1 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                              >  2:362215/1‑99 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                             >  1:231740/1‑100 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                             <  1:1339017/100‑1 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                             <  1:456314/100‑1 (MQ=255)
                                         ttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATcc                                                            >  1:1533471/1‑100 (MQ=255)
                                          tCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCg                                                           <  2:1679990/100‑1 (MQ=255)
                                           cAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGt                                                          <  1:1953768/100‑1 (MQ=255)
                                            aGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTCTTCGGCGAAAGAAAAGATCAGt                                                          >  2:3351040/1‑99 (MQ=255)
                                                cTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGTCGAAAGAAAAGATCCGTAATGc                                                     >  2:3200133/1‑100 (MQ=255)
                                                 ttttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGcc                                                    <  1:2906758/100‑1 (MQ=255)
                                                     aCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGAt                                                >  2:1268369/1‑100 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                               >  2:1293757/1‑100 (MQ=255)
                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATACGTAATGCCTGATgcg                                             >  2:886836/1‑100 (MQ=255)
                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGGAATGACTGGGGCGTTTTTTATTGGGCGAAAGAATAGATCCTTAATGCCTGATgcg                                             <  2:1353815/100‑1 (MQ=255)
                                                          aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                           <  2:2302910/100‑1 (MQ=255)
                                                            gTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTAt                                         >  1:3342238/1‑100 (MQ=255)
                                                             tAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                        <  2:1531297/100‑1 (MQ=255)
                                                               aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGtt                                      >  1:1744969/1‑100 (MQ=255)
                                                                 cGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTa                                    >  2:1955169/1‑100 (MQ=255)
                                                                 cGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTa                                    <  2:475784/100‑1 (MQ=255)
                                                                   gATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                  <  2:3354060/100‑1 (MQ=255)
                                                                    aTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCa                                 <  1:270318/100‑1 (MQ=255)
                                                                     tttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                <  1:3614957/100‑1 (MQ=255)
                                                                     tttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTGTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                >  1:1075502/1‑100 (MQ=255)
                                                                      ttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAgg                               <  1:3124544/100‑1 (MQ=255)
                                                                      ttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAgg                               <  2:3011245/100‑1 (MQ=255)
                                                                      ttGGATGAATATTTTAGATGAACCAGTCAGTAATGACTGAGGTGTTTGTTATTGGGCGCAAGAAAAGATCGGTAATGCCCGATGAGCTATGTTTATCAgg                               <  2:2171949/100‑1 (MQ=255)
                                                                       tGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGc                              <  2:818203/100‑1 (MQ=255)
                                                                        ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                             >  1:1623136/1‑100 (MQ=255)
                                                                          cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGATTATCAGGCCaa                           >  2:1428321/1‑100 (MQ=255)
                                                                           tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                           >  2:2785203/1‑99 (MQ=255)
                                                                           tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAAc                          <  2:1884729/100‑1 (MQ=255)
                                                                              atatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGt                       <  1:1115249/100‑1 (MQ=255)
                                                                               tatTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTa                      >  1:1732207/1‑100 (MQ=255)
                                                                                atTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAg                     <  2:1323970/100‑1 (MQ=255)
                                                                                 ttttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGa                    >  2:1516033/1‑100 (MQ=255)
                                                                                  tttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGaa                   <  2:71602/100‑1 (MQ=255)
                                                                                   ttAGCCCCCACAGTCAGTAATGACTGGGGCAATTTTTACCGCCAGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGTCAACGGTAGaa                   <  2:3087906/99‑1 (MQ=255)
                                                                                   aaaaacgccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAt                  >  1:1318812/6‑100 (MQ=38)
                                                                                    aaaacgccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAAtt                 >  1:39104/5‑100 (MQ=38)
                                                                                      gccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGt               <  2:2734754/100‑1 (MQ=255)
                                                                                       ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa              >  1:647285/1‑100 (MQ=255)
                                                                                       acgccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa              >  1:2498686/2‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:2047488/1‑100 (MQ=255)
                                                                                         gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTTTTTTTTATTAAAGATCCGTAATGCCTGATGCTCTATTTTTATCAGGCTCACGGTAGAATTGTTAt            >  2:3103000/1‑100 (MQ=255)
                                                                                         gccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt            >  2:1056582/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  2:2746008/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  2:22227/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:2981023/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:169044/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTACAATTGTAATc           >  1:2672519/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCAGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:26018/1‑100 (MQ=255)
                                                                                            ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa         >  1:3510048/1‑100 (MQ=255)
                                                                                            ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa         >  1:3146945/1‑100 (MQ=255)
                                                                                            ccAGTCATAACTGCCTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa         >  2:2923789/1‑100 (MQ=255)
                                                                                             cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAt        <  1:1405365/100‑1 (MQ=255)
                                                                                             cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAt        >  2:2690291/1‑100 (MQ=255)
                                                                                             cagtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAt        <  1:1041762/100‑1 (MQ=255)
                                                                                              agtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt       <  2:2611849/100‑1 (MQ=255)
                                                                                                 cagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGaa    >  2:3305153/1‑100 (MQ=255)
                                                                                                  agtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt   >  1:320792/1‑100 (MQ=255)
                                                                                                   gtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAtt  <  1:2800655/100‑1 (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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: