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.C98.6% 976.4 / 3.9 281coding (1772/1776 nt)yihKputative GTP‑binding factor
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base C (152/125);  total (152/129)
Fisher's exact test for biased strand distribution p-value = 4.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 303 total aligned reads displayed.

TATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATG  >  NC_000913/4057666‑4057845
                                                                                             |                                                                                        
tATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCa                                                                                           >  1:471230/1‑93 (MQ=255)
 aTCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAc                                                                                         >  1:1984557/1‑93 (MQ=255)
  tCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcc                                                                                        >  2:5940317/1‑92 (MQ=255)
   ccGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGATCTCTATCCGTATTCGTAAACGTCAcct                                                                                       >  1:6832673/1‑91 (MQ=255)
   ccGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCAcct                                                                                       >  1:100339/1‑91 (MQ=255)
        tGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACgg                                                                                  <  1:3063232/93‑1 (MQ=255)
        tGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACgg                                                                                  >  1:6434318/1‑93 (MQ=255)
        tGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACgg                                                                                  >  2:8667198/1‑93 (MQ=255)
         gACTCTGGAACAAGCTCTGGAGTTCATCGATGATGGCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGa                                                                                 >  2:5608298/1‑93 (MQ=255)
         gACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGa                                                                                 <  1:7289392/93‑1 (MQ=255)
           ctctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaa                                                                               >  1:7804420/1‑93 (MQ=255)
            tctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaaa                                                                              >  1:1920597/1‑93 (MQ=255)
            tctGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGaaaa                                                                              >  1:6669713/1‑93 (MQ=255)
               ggAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGa                                                                           >  2:2144960/1‑93 (MQ=255)
               ggAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGa                                                                           <  1:8365313/93‑1 (MQ=255)
                gAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAt                                                                          <  1:1754273/93‑1 (MQ=255)
                gAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAt                                                                          >  1:5589265/1‑93 (MQ=255)
                 aaCAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtc                                                                         <  2:5820617/93‑1 (MQ=255)
                  aCAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGAtcg                                                                        <  2:6986809/93‑1 (MQ=255)
                       ctctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                                   <  1:6374171/93‑1 (MQ=255)
                       ctctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgcc                                                                   >  1:1316947/1‑93 (MQ=255)
                        tctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccg                                                                  <  2:3568324/93‑1 (MQ=255)
                        tctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccg                                                                  >  1:3971007/1‑93 (MQ=255)
                         ctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccgc                                                                 >  2:8111032/1‑93 (MQ=255)
                         ctGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTcgccgc                                                                 >  2:3962923/1‑93 (MQ=255)
                            gAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCCCCGCGcc                                                              >  2:6211691/1‑93 (MQ=255)
                             aGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCa                                                             >  2:909998/1‑93 (MQ=255)
                               ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAc                                                           <  2:3703937/93‑1 (MQ=255)
                               ttCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAc                                                           <  1:2335845/93‑1 (MQ=255)
                               ttCATCGATGAGGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAAc                                                           <  1:2443146/93‑1 (MQ=255)
                                 cATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcg                                                         >  2:4796932/1‑93 (MQ=255)
                                   tCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACcgcg                                                       >  2:4905266/1‑93 (MQ=255)
                                     gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACTTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCa                                                     >  2:7707348/1‑93 (MQ=255)
                                     gatgatGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCa                                                     >  2:5325095/1‑93 (MQ=255)
                                      atgatgACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCAc                                                    >  1:7634143/1‑93 (MQ=255)
                                         atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGa                                                 >  1:5955225/1‑93 (MQ=255)
                                         atgaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCACCGCGCCAACCGCGCACCGa                                                 >  2:6540493/1‑93 (MQ=255)
                                           gaCGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGaaa                                               <  1:8581896/93‑1 (MQ=255)
                                              gAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGAc                                            <  2:1864084/93‑1 (MQ=255)
                                              gAACTGGTAGAAGTGACTCCGACATCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGAc                                            <  2:4768325/93‑1 (MQ=255)
                                                aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                                          <  1:3888747/93‑1 (MQ=255)
                                                aCTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGa                                          <  2:8782478/93‑1 (MQ=255)
                                                 cTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGAt                                         <  2:4205422/93‑1 (MQ=255)
                                                   ggTAGAAGTGAGTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTa                                       >  2:5057371/1‑93 (MQ=255)
                                                    gTAGAAGGTCCTCCCACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTaa                                      <  1:7848607/93‑1 (MQ=255)
                                                     tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAt                                     >  2:6768152/1‑93 (MQ=255)
                                                     tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAt                                     >  2:8221688/1‑93 (MQ=255)
                                                     tAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAt                                     <  2:1662843/93‑1 (MQ=255)
                                                       gAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAttt                                   >  1:1510570/1‑93 (MQ=255)
                                                       gAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAttt                                   >  2:4432016/1‑93 (MQ=255)
                                                       gAAGTGACCCCGGCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAAttt                                   <  2:8492777/93‑1 (MQ=255)
                                                        atgTGACTCCGACCTCTATCCGTTTTCGTAAACGTCACCTGACGGAATACGATCGTAGCCTCGCCAAGCGCGCACCGAAAGCCGGTTAAtttg                                  >  2:6487353/3‑92 (MQ=255)
                                                         aGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATGTCg                                 >  2:7855139/1‑93 (MQ=255)
                                                          gTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGt                                <  2:2921834/93‑1 (MQ=255)
                                                          gTGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAAAGATCGTCGCCGCGCCAACCACGCACCGAAAGACGCTTAATTTCGt                                >  2:8728646/1‑93 (MQ=255)
                                                           tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  2:8457380/93‑1 (MQ=255)
                                                           tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  2:2637534/93‑1 (MQ=255)
                                                           tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  1:1117704/93‑1 (MQ=255)
                                                           tGACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  1:8667198/93‑1 (MQ=255)
                                                           ggACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtt                               <  1:7059314/92‑1 (MQ=255)
                                                            gACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttt                              <  1:6132625/93‑1 (MQ=255)
                                                            gACTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACAGCGCACCGAAAGACGATTAATTTCGttt                              >  1:3348233/1‑93 (MQ=255)
                                                             aCTCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttt                             >  2:4720920/1‑93 (MQ=255)
                                                              cTCCGCCCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGttttt                            <  1:8089778/93‑1 (MQ=255)
                                                               tCCGACCTCTATCCGTATTCGTAAACGTCACCTGACGGTAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGtttttt                           >  1:4373247/1‑93 (MQ=255)
                                                                ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTa                          >  2:7421226/1‑93 (MQ=255)
                                                                ccGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTa                          <  1:6049819/93‑1 (MQ=255)
                                                                 cGACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAg                         <  2:1452784/93‑1 (MQ=255)
                                                                  gACCTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGt                        >  1:5674571/1‑93 (MQ=255)
                                                                    ccTCTATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCa                      >  1:3773611/1‑93 (MQ=255)
                                                                     ctctATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAt                     <  1:6902497/93‑1 (MQ=255)
                                                                        tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaa                  >  2:3778134/1‑93 (MQ=255)
                                                                        tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaa                  <  2:7309010/93‑1 (MQ=255)
                                                                        tATCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaa                  >  2:121076/1‑93 (MQ=255)
                                                                         aTCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaa                 <  2:3165236/93‑1 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                >  1:2158101/1‑93 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                >  1:5533996/1‑93 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                >  1:3200678/1‑93 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                >  1:2553565/1‑93 (MQ=255)
                                                                          tCCGTATTCGTAAACGTCACATGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATaaaaa                <  1:3561313/93‑1 (MQ=255)
                                                                            cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAc              >  2:1989600/1‑93 (MQ=255)
                                                                            cGTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAc              <  1:7212344/93‑1 (MQ=255)
                                                                             gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTTATAAAAAAcc             >  1:8199443/1‑93 (MQ=255)
                                                                             gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAcc             <  2:5704675/93‑1 (MQ=255)
                                                                             gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAcc             <  1:8814918/93‑1 (MQ=255)
                                                                             gTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAcc             <  1:8456801/93‑1 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACTATTAATTTCGTTTTTTAGTCATAAAAAACCt            >  2:1787001/1‑93 (MQ=255)
                                                                              tATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCt            >  2:667489/1‑93 (MQ=255)
                                                                               aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTg           <  1:7172618/93‑1 (MQ=255)
                                                                                ttCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGc          <  2:2801003/93‑1 (MQ=255)
                                                                                 tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGcc         >  2:3357647/1‑93 (MQ=255)
                                                                                 tCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGcc         >  2:2818863/1‑93 (MQ=255)
                                                                                 tCGTAAACGTCACCTGACGAAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGcc         <  2:806473/93‑1 (MQ=255)
                                                                                  cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCa        >  2:1025446/1‑93 (MQ=255)
                                                                                   gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAg       <  2:825394/93‑1 (MQ=255)
                                                                                     aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCg     <  1:4899952/93‑1 (MQ=255)
                                                                                     aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCg     >  1:7366864/1‑93 (MQ=255)
                                                                                      aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa    >  1:8733594/1‑93 (MQ=255)
                                                                                       aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGAt   <  1:5150350/93‑1 (MQ=255)
                                                                                        cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATg  <  1:547107/93‑1 (MQ=255)
                                                                                             |                                                                                        
TATCCGCATGACTCTGGAACAAGCTCTGGAGTTCATCGATGATGACGAACTGGTAGAAGTGACTCCGACCTCTATCCGTATTCGTAAACGTCACTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATG  >  NC_000913/4057666‑4057845

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