Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,049,986 G→T 15.4% 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.4% 455.5 / 64.0 175intergenic (+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 (75/73);  new base T (15/12);  total (90/85)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-02

Only 100 of 192 total aligned reads displayed.

AAAGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAATT  >  NC_000913/4049887‑4050085
                                                                                                   |                                                                                                   
aaaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                     >  2:2558488/1‑100 (MQ=255)
aaaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTAACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcag                                                                                                     >  1:1330881/1‑100 (MQ=255)
 aaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagt                                                                                                    <  1:1074262/100‑1 (MQ=255)
 aaGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCcagtcagt                                                                                                    <  2:1880399/100‑1 (MQ=255)
  aGTTAGTCGCGTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTa                                                                                                   >  1:1204699/1‑100 (MQ=255)
          gcgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgg                                                                                           <  2:3310293/100‑1 (MQ=255)
           cgTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTggg                                                                                          <  2:846355/100‑1 (MQ=255)
            gTAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTgggg                                                                                         <  1:1913746/100‑1 (MQ=255)
            gTAGGGTACACAGGTAAGATGTTCTATCGTTCAGACCTTTTACTTCAGGTATTCGGATTTGGTTGATTATTTTAGCCTCCATAGTCATTACTGAGTgggg                                                                                         <  2:3457667/100‑1 (MQ=255)
             tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                        <  1:41184/100‑1 (MQ=255)
             tAGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGc                                                                                        >  1:3116252/1‑100 (MQ=255)
              aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCATTACTGACTGGGGc                                                                                        <  2:3676497/99‑1 (MQ=255)
              aGGGTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCg                                                                                       <  2:304511/100‑1 (MQ=255)
               gggTACAGATGTAAGATGTTCTATATTTCAGTCTTTTTACGTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCTTAGTCATTACTGACTGGGGCgg                                                                                      <  2:3840796/100‑2 (MQ=255)
               gggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGAATTTAGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGt                                                                                      >  1:132496/1‑100 (MQ=255)
                ggTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtt                                                                                     >  1:2143348/1‑100 (MQ=255)
                 gTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGCTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATAACTGGGGc                                                                                        >  1:896626/1‑96 (MQ=255)
                 gTACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAAGCGGATTTGGCTGAATATTTTAGCCGCCCTAGTCATTACTGACTGGGGCggct                                                                                    <  2:2876879/100‑4 (MQ=38)
                  tACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttt                                                                                   >  1:3233577/1‑100 (MQ=255)
                  tACAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttt                                                                                   <  1:2195574/100‑1 (MQ=255)
                   aCAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGttttt                                                                                  <  1:2859900/100‑1 (MQ=255)
                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGTCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                 <  2:3582016/100‑1 (MQ=255)
                    cAGAGGTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGtttttt                                                                                 <  2:2227619/100‑1 (MQ=255)
                      gagGTAAGATGTTCTATCTTTCAGACCTTTTACTTCATGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTAt                                                                               <  2:394359/100‑1 (MQ=255)
                        ggTAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTg                                                                             >  2:1160771/1‑100 (MQ=255)
                         gTATGATGTTCTATCGTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTgg                                                                            <  1:3874928/100‑1 (MQ=255)
                          tAAGATGTTCTATCTTTCAGACCTTTTATTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                           <  1:2709246/100‑1 (MQ=255)
                          tAAGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTggg                                                                           <  1:3277413/100‑1 (MQ=255)
                           aaGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc                                                                          >  2:1577735/1‑100 (MQ=255)
                           aaGATGTTCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGc                                                                          <  1:1985358/100‑1 (MQ=255)
                                 ttCTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgaaaga                                                                    >  2:913970/1‑100 (MQ=255)
                                   cTATCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCgtaag                                                                     >  2:1409257/1‑97 (MQ=255)
                                     aTCTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGAGGCGTTTTTTATTGGGCGaaagaaaa                                                                 >  1:2585501/1‑99 (MQ=255)
                                       cTTTCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGAt                                                              >  2:3604559/1‑100 (MQ=255)
                                        tttCAGACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATc                                                             >  2:652297/1‑100 (MQ=255)
                                             gACCTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTaa                                                        >  1:3844216/1‑100 (MQ=255)
                                               ccTTTTACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATg                                                      >  2:3631663/1‑100 (MQ=255)
                                                 ttttACTTCATGTAATCGGATTTGGCTGAATATTTTAGCCGCCACAGTCAGTAATGACTTGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGcc                                                    <  1:3845116/100‑1 (MQ=255)
                                                   ttACTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCGCAGTCAGTTATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTg                                                  <  1:1116620/100‑1 (MQ=255)
                                                     aCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGAt                                                <  2:3199346/100‑1 (MQ=255)
                                                     aCTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATACGTAATGcat                                                   >  2:3516769/1‑95 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTATCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATg                                               <  1:14947/100‑1 (MQ=255)
                                                      cTTCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGATTTTTATTGGGAGAAAGAAAAGATCCGTAATACCTGa                                                 >  1:2180614/1‑98 (MQ=255)
                                                       ttCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgc                                              <  2:3065552/100‑1 (MQ=255)
                                                       ttCACGTAATCGAATTTGGCTGAATAATTTAGCCGACCCAGTAAGTAATGACTGGGGCGTTTTTTATTGGGCGGAAGAAAAGATCCGTAATGCCTGAt                                                >  2:3214223/1‑98 (MQ=255)
                                                        tCACTTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTCTTTTAGTGGGCGTAACAAAAGATCCGTAATg                                                      >  2:2810733/1‑91 (MQ=255)
                                                        tCACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcg                                             >  2:3293025/1‑100 (MQ=255)
                                                         cACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcgc                                            >  2:2692641/1‑100 (MQ=255)
                                                         cACGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTAGGCGAAAGAAAAGATCCGTAATGCCTGATgcgc                                            >  1:1256405/1‑100 (MQ=255)
                                                         cACGGAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATgcgc                                            <  2:2014534/100‑1 (MQ=255)
                                                          aCGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCt                                           <  1:3718753/100‑1 (MQ=255)
                                                           cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTa                                          >  2:2178991/1‑100 (MQ=255)
                                                           cGTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGCCGTAAGAAAAGATCCGTAATGCCTGATGCGCTa                                          <  2:2578381/100‑1 (MQ=255)
                                                            gTAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTAt                                         >  2:1450611/1‑100 (MQ=255)
                                                             tAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                        <  2:1226219/100‑1 (MQ=255)
                                                             tAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                        <  2:3908599/100‑1 (MQ=255)
                                                             tAATCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATg                                        <  1:505014/100‑1 (MQ=255)
                                                              aaTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGt                                       >  2:3870411/1‑100 (MQ=255)
                                                              aaTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGt                                       <  2:2364039/100‑1 (MQ=255)
                                                               aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGtt                                      >  2:1482497/1‑100 (MQ=255)
                                                               aTCGGATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGtt                                      >  1:3919458/1‑100 (MQ=255)
                                                                   gATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                  >  2:1939544/1‑100 (MQ=255)
                                                                   gATTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATc                                  <  2:561274/100‑1 (MQ=255)
                                                                    aTTTGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCAGTAATGCCTGATGCGCTATGTTTATCa                                 >  2:44106/1‑100 (MQ=255)
                                                                     tttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAg                                <  1:1773839/100‑1 (MQ=255)
                                                                      ttGGCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAgg                               <  1:2068643/100‑1 (MQ=255)
                                                                        ggCTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGcc                             <  2:177640/100‑1 (MQ=255)
                                                                          cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGGGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCa                            <  1:1685300/99‑1 (MQ=255)
                                                                          cTGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCaa                           >  1:2367653/1‑100 (MQ=255)
                                                                           tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAAc                          >  2:1474778/1‑100 (MQ=255)
                                                                           tGAATATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAAc                          <  1:3516808/100‑1 (MQ=255)
                                                                             aaTATTTTAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACgg                        <  1:3293059/100‑1 (MQ=255)
                                                                                  tttAGCCGCCCCAGTCAGTAATGACTGGGGCGTTTTTTTTTTTTTTATACAAAAGATACGTAATGATTGATGCGATATGTTCATAAggta                             >  2:2649916/1‑88 (MQ=255)
                                                                                     aGCCGCCCCAGTCAGTAATGACTGGGGTGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTg                <  1:1158533/100‑1 (MQ=255)
                                                                                       ccgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa              <  1:2616772/100‑1 (MQ=255)
                                                                                        cgccCCAGTCATTACTGACTGGGGCGTTTTTTATTTGGCTAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa             >  2:3900978/1‑100 (MQ=39)
                                                                                        cgccCCAGTCATTACTGACTGGGGCGTTTTTTATTTGGCGTAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTa              >  2:466548/1‑99 (MQ=38)
                                                                                        cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa             >  2:1344339/1‑100 (MQ=255)
                                                                                        cgccCCAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa             <  2:874357/100‑1 (MQ=255)
                                                                                        cgccCCAGTCAGTAATGACTGGGACGATTTATATTGTTCGAAATAAAAAAACCTTAATGCCTGATAAGCTATGAATATCAGGCATACGGTAGTATTGTaa             >  2:3120030/1‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTTTTTTTAAGAAAAGAACCGTAATGCCTGATTCGCTATGTTTATCAGGCCAACGCTAGACTTGTa              >  2:1760554/1‑98 (MQ=19)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt            >  2:3402140/1‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt            >  2:1117499/1‑100 (MQ=255)
                                                                                         gccCCAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAAt            >  1:3767589/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTaa             >  2:3772691/1‑98 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:3365226/1‑100 (MQ=255)
                                                                                          ccccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCACTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:1774482/1‑100 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:1387741/1‑100 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           <  2:2353060/100‑1 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:652307/1‑100 (MQ=255)
                                                                                          ccccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATc           >  1:308231/1‑100 (MQ=255)
                                                                                           cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCt          >  1:1644220/1‑100 (MQ=255)
                                                                                           cccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCt          >  1:968738/1‑100 (MQ=255)
                                                                                            ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCt          >  2:588960/1‑99 (MQ=255)
                                                                                            ccAGTCATTACTGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa         >  2:780146/1‑100 (MQ=255)
                                                                                            ccAGTCAGTAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTa         >  2:3281174/1‑100 (MQ=255)
                                                                                              agtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTAtt       >  1:3628867/1‑100 (MQ=255)
                                                                                               gtcagtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTg      <  1:3006130/100‑1 (MQ=255)
                                                                                                  agtAATTATTTGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCGTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt   <  1:225628/100‑1 (MQ=255)
                                                                                                  agtAATGACTGGGGCGTTTTTTATTGGGCGAAAGAAAAGATCCGTAATGCCTGATGCGCTATGTTTATCAGGCCAACGGTAGAATTGTAATCTATTGAAt   >  1:102423/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 < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: