Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,057,759 +C coding (1773/1776 nt) yihK → putative GTP‑binding factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,057,7581.C99.5% 696.2 / ‑4.1 211coding (1772/1776 nt)yihKputative GTP‑binding factor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (120/90);  total (120/91)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 243 total aligned reads displayed.

TATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATG  >  NC_000913/4057666‑4057845
                                                                                              |                                                                                       
tATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCa                                                                                          >  2:1979873/1‑93 (MQ=255)
tATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCa                                                                                          >  1:2160901/1‑93 (MQ=255)
 aTCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTTGTAAACGTCAc                                                                                        <  2:4620016/93‑1 (MQ=255)
 aTCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAc                                                                                        >  1:1181030/1‑93 (MQ=255)
 aTCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAc                                                                                        >  1:1802712/1‑93 (MQ=255)
  tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                                       <  2:109603/93‑2 (MQ=255)
  tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                                       >  1:4550327/1‑92 (MQ=255)
  tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                                       >  2:4706196/1‑92 (MQ=255)
   ccGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcct                                                                                      >  2:60871/1‑91 (MQ=255)
    cGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctg                                                                                     <  2:1322909/93‑4 (MQ=255)
    cGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcctg                                                                                     >  2:3952805/1‑90 (MQ=255)
       aTGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACg                                                                                   >  2:2259977/1‑93 (MQ=255)
        tGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACgg                                                                                  >  2:4003940/1‑93 (MQ=255)
            tctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaaa                                                                              >  1:3071864/1‑93 (MQ=255)
            tctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaaa                                                                              >  2:135391/1‑93 (MQ=255)
             ctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACGGGTAGAAGTGACTCCGACCTCTATCCGTATGCGTAAACGTCACCTGACGGAAAAc                                                                             <  2:3270606/93‑1 (MQ=255)
                gAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAt                                                                          <  2:4729608/93‑1 (MQ=255)
                gAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAt                                                                          <  2:1904078/93‑1 (MQ=255)
                gAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAt                                                                          >  1:414182/1‑93 (MQ=255)
                 aaCAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                                         >  2:4656021/1‑93 (MQ=255)
                 aaCAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                                         <  1:1862653/93‑1 (MQ=255)
                 aaCAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                                         <  1:3714057/93‑1 (MQ=255)
                    aaGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                                      <  1:1009744/93‑1 (MQ=255)
                    aaGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                                      >  2:2804256/1‑93 (MQ=255)
                    aaGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcgtc                                                                      <  2:3267241/93‑1 (MQ=255)
                        tctGGTGGTCATCGATGATGACGAACCGGTAGAAGTGACGCCGACCTCTATCCGTAAGCGCAAACGCCACCTGACGGAAAACGAGCGAcgccg                                                                  >  2:4206800/1‑93 (MQ=255)
                           ggAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCgcgc                                                               <  2:2503967/93‑1 (MQ=255)
                           ggAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCgcgc                                                               <  2:1377944/93‑1 (MQ=255)
                              gTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCaa                                                            >  1:4487436/1‑93 (MQ=255)
                               ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAc                                                           >  1:2975886/1‑93 (MQ=255)
                                 cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcg                                                         >  1:773163/1‑93 (MQ=255)
                                  aTCGATGATGACGAACTGGTAGAAGCGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGACGCCGCGCCAACcgc                                                        >  1:2477854/1‑93 (MQ=255)
                                   tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgcg                                                       >  1:1263791/1‑93 (MQ=255)
                                     gatgatGACTACCTGCTAGACGTCACTCCGCCCTCTATCCGTATCCGTAAACGTCCCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCa                                                     <  2:257854/93‑1 (MQ=255)
                                     gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCa                                                     >  1:2228323/1‑93 (MQ=255)
                                     gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCa                                                     <  2:1148371/93‑1 (MQ=255)
                                     gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCa                                                     <  2:2189832/93‑1 (MQ=255)
                                      atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAc                                                    <  1:4536983/93‑1 (MQ=255)
                                       tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                                   <  2:4462287/93‑1 (MQ=255)
                                       tgatgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAcc                                                   >  2:4185736/1‑93 (MQ=255)
                                         atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGa                                                 >  2:3672841/1‑93 (MQ=255)
                                           gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaaa                                               >  1:2853954/1‑93 (MQ=255)
                                               aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                                           >  1:665105/1‑93 (MQ=255)
                                               aaCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACg                                           <  2:883811/93‑1 (MQ=255)
                                                aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                                          <  1:1140558/93‑1 (MQ=255)
                                                aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                                          <  1:3152260/93‑1 (MQ=255)
                                                 cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAGCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                                         <  2:3648005/93‑1 (MQ=255)
                                                 cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                                         <  1:3964142/93‑1 (MQ=255)
                                                   ggTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTa                                       >  1:1625624/1‑93 (MQ=255)
                                                    gTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTaa                                      >  2:1207680/1‑93 (MQ=255)
                                                     tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGAGCACCGAAAGACGATTAAt                                     >  1:4617718/1‑93 (MQ=255)
                                                      aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAtt                                    >  1:4387081/1‑93 (MQ=255)
                                                      aGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAtt                                    <  1:2000895/93‑1 (MQ=255)
                                                       gAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAttt                                   <  1:823348/93‑1 (MQ=255)
                                                        aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTc                                  <  1:681679/93‑1 (MQ=255)
                                                        aaGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTc                                  <  2:2179014/93‑1 (MQ=255)
                                                             aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttt                             >  2:936177/1‑93 (MQ=255)
                                                             aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttt                             <  1:1561506/93‑1 (MQ=255)
                                                              cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt                            >  1:1358294/1‑93 (MQ=255)
                                                              cTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt                            <  1:2063302/93‑1 (MQ=255)
                                                               tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttttt                           <  2:4293255/93‑1 (MQ=255)
                                                                 cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAg                         <  2:4050739/93‑1 (MQ=255)
                                                                  gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGt                        >  1:3794071/1‑93 (MQ=255)
                                                                  gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGt                        >  2:3784558/1‑93 (MQ=255)
                                                                   aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTc                       <  1:4316944/93‑1 (MQ=255)
                                                                   aCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTc                       >  2:3339314/1‑93 (MQ=255)
                                                                    ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa                      >  2:158259/1‑93 (MQ=255)
                                                                    ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa                      <  2:3714626/93‑1 (MQ=255)
                                                                    ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa                      <  1:3352041/93‑1 (MQ=255)
                                                                    ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa                      <  2:3118190/93‑1 (MQ=255)
                                                                     ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCGCCGAAAGACGATTAATTTCGTTTTTTAGTCAt                     >  1:881763/1‑93 (MQ=255)
                                                                     ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAt                     <  2:2779414/93‑1 (MQ=255)
                                                                     ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAt                     <  2:3487056/93‑1 (MQ=255)
                                                                      tctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATa                    >  2:2659733/1‑93 (MQ=255)
                                                                       ctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCGCCGAAAGACGATTAATTTCGTTTTTTAGTCATaa                   >  1:3561338/1‑93 (MQ=255)
                                                                        tATCCGTATTCGTAAAGGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaa                  >  2:4507563/1‑93 (MQ=255)
                                                                        tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaa                  >  2:1613370/1‑93 (MQ=255)
                                                                         aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaca                 >  1:4054704/1‑91 (MQ=255)
                                                                         aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaa                 <  1:2362576/93‑1 (MQ=255)
                                                                         aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaa                 <  1:3309010/93‑1 (MQ=255)
                                                                         aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaa                 <  2:1767689/93‑1 (MQ=255)
                                                                          tCCTTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                <  2:638819/93‑1 (MQ=255)
                                                                           ccGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaaa               <  2:1689624/93‑1 (MQ=255)
                                                                           ccGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaaa               <  2:4383399/93‑1 (MQ=255)
                                                                            cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAc              >  2:1860628/1‑93 (MQ=255)
                                                                            cGTATGCGTAAACGTCACCTGGCGGAAAACGATCGTCGGCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAc              <  2:3385097/93‑1 (MQ=255)
                                                                             gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAcc             >  1:3905847/1‑93 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCt            <  1:3399619/93‑1 (MQ=255)
                                                                               aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTg           >  1:2409698/1‑93 (MQ=255)
                                                                               aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTg           <  2:2109598/93‑1 (MQ=255)
                                                                               aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTg           >  1:1244594/1‑93 (MQ=255)
                                                                                ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACTATTAATTTCGTTTTTTAGTCATAAAAAACCTGc          <  2:4553696/93‑1 (MQ=255)
                                                                                ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGc          <  2:3730667/93‑1 (MQ=255)
                                                                                 tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGcc         <  2:356649/93‑1 (MQ=255)
                                                                                  cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCa        >  1:2565122/1‑93 (MQ=255)
                                                                                     aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCg     <  2:2341412/93‑1 (MQ=255)
                                                                                     aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCg     <  1:1782153/93‑1 (MQ=255)
                                                                                       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt   <  1:2196328/93‑1 (MQ=255)
                                                                                       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt   >  2:1258597/1‑93 (MQ=255)
                                                                                        cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATg  <  2:2489282/93‑1 (MQ=255)
                                                                                              |                                                                                       
TATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATG  >  NC_000913/4057666‑4057845

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