Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,689,553 A→G intergenic (+169/+57) ydgA → / ← uidC hypothetical protein/membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,689,5530AG99.6% 701.6 / NA 223intergenic (+169/+57)ydgA/uidChypothetical protein/membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (99/123);  total (99/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

GCTTATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTCAT  >  NC_000913/1689463‑1689642
                                                                                          |                                                                                         
gCTTATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGcg                                                                                         <  1:3733286/93‑1 (MQ=255)
 cTTATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAgcgc                                                                                        >  2:6700017/1‑93 (MQ=255)
  ttATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGcgcg                                                                                       >  1:1731962/1‑93 (MQ=255)
  ttAGCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGcgcg                                                                                       >  2:8928425/1‑93 (MQ=255)
   tATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAgcgcgc                                                                                      <  1:4065734/93‑1 (MQ=255)
    aTCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCa                                                                                     <  1:1380029/93‑1 (MQ=255)
     tCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCAc                                                                                    <  2:9503858/93‑1 (MQ=255)
     tCAGCCCTACCTTTTCCTGCAACACTTTGAATTTAGGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCGGCATCCGGCATAAAAAACGCGCAc                                                                                    <  2:4562823/93‑1 (MQ=25)
      cAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACt                                                                                   >  1:715441/1‑93 (MQ=255)
      cAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACt                                                                                   >  1:8111777/1‑93 (MQ=255)
       aGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACtt                                                                                  <  1:7930552/93‑1 (MQ=255)
       aGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACtt                                                                                  >  2:7101286/1‑93 (MQ=255)
        ggCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACttt                                                                                 >  2:3609064/1‑93 (MQ=255)
         gCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTg                                                                                >  1:5937681/1‑93 (MQ=255)
          ccTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGt                                                                               >  1:8999804/1‑93 (MQ=255)
          ccTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGt                                                                               <  2:2902522/93‑1 (MQ=255)
          ccTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGt                                                                               <  1:4226115/93‑1 (MQ=255)
            tACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTca                                                                             <  2:2814007/93‑1 (MQ=255)
            tACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTca                                                                             >  1:1208081/1‑93 (MQ=255)
            tACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTca                                                                             >  2:4623394/1‑93 (MQ=255)
            tACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCAATTTGTca                                                                             >  2:808796/1‑93 (MQ=255)
            tACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGATGGATAAGGCGTTTACGCCGACTCCGGCATACAAAAAGCGCACTTTGtta                                                                             >  1:7027261/1‑91 (MQ=25)
            ctacTTTTTCTGCCACCCCTTTGATTTATGAGTTTTTGGAGGCTGGGAAAAGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTca                                                                             <  1:4350147/90‑1 (MQ=12)
             aCCTTTTCCTGCAACACTTTGCATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaa                                                                            >  1:1979737/1‑93 (MQ=255)
              ccTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaac                                                                           <  2:975739/93‑1 (MQ=255)
               cTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaaca                                                                          >  1:915669/1‑93 (MQ=255)
                ttttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaacaa                                                                         <  2:1718417/93‑1 (MQ=255)
                ttttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaacaa                                                                         >  1:5842465/1‑93 (MQ=255)
                ttttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaacaa                                                                         <  1:164525/93‑1 (MQ=255)
                ttttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTcaacaa                                                                         <  2:5972495/93‑1 (MQ=255)
                  ttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATc                                                                       <  2:6543012/93‑1 (MQ=255)
                  ttCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATc                                                                       <  2:3785417/93‑1 (MQ=255)
                  ttCCTGCAACACTTTGAATTTATGAGTGTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCTGCATAAAAAGCGCGCACTTTGTCAACAATc                                                                       <  2:8650444/93‑1 (MQ=255)
                      tGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaa                                                                   >  2:2264694/1‑93 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  <  2:1942812/93‑1 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  <  2:5837675/93‑1 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  <  1:4478099/93‑1 (MQ=255)
                       gCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGaaa                                                                  >  2:8303311/1‑93 (MQ=255)
                        cAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAAc                                                                 <  1:3864828/93‑1 (MQ=255)
                        cAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGAATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAAc                                                                 <  1:4604274/93‑1 (MQ=255)
                            acTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCg                                                             <  1:6415480/93‑1 (MQ=255)
                            acTTCGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCg                                                             >  2:1898095/1‑93 (MQ=255)
                              tttGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCCCGCACTTTGTCAACAATCCGAAACGCCGga                                                           <  2:122972/93‑1 (MQ=255)
                                tGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGgaga                                                         <  2:7358021/93‑1 (MQ=255)
                                tGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGgaga                                                         <  2:8245408/93‑1 (MQ=255)
                                tGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGgaga                                                         >  1:9171613/1‑93 (MQ=255)
                                 gAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGAt                                                        <  2:4121415/93‑1 (MQ=255)
                                  aaTTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGAtt                                                       <  2:5685192/93‑1 (MQ=255)
                                     ttATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGAttttt                                                    >  1:2159653/1‑93 (MQ=255)
                                     ttATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGAttttt                                                    <  2:3525458/93‑1 (MQ=255)
                                       aTGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtct                                                  >  1:9995756/1‑93 (MQ=255)
                                       aTGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtct                                                  <  1:7599630/93‑1 (MQ=255)
                                          aGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtctctc                                               >  1:7513759/1‑93 (MQ=255)
                                          aGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtctctc                                               <  1:3461024/93‑1 (MQ=255)
                                          aGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTtctctc                                               <  2:9450900/93‑1 (MQ=255)
                                           gTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTcc                                              <  1:6196673/93‑1 (MQ=255)
                                           gTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTcc                                              >  2:936639/1‑93 (MQ=255)
                                            tttttGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCg                                             <  2:3209732/93‑1 (MQ=255)
                                              tttGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCTGGAGATTTTTCTCTCCGGc                                           >  2:7756311/1‑93 (MQ=255)
                                              tttGTAGCCGGGAGAGGGCTTTCACGCCGCATCCGGAAAAAAAGGCGCGCTTTTTGCCAACAATCTAAACCGCCGGATATTTTTCTCTCCGGc                                           <  2:5443921/93‑1 (MQ=12)
                                               ttGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCg                                          <  1:2336276/93‑1 (MQ=255)
                                                 gTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGtt                                        <  1:2665459/93‑1 (MQ=255)
                                                     gCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAttt                                    >  2:8780299/1‑93 (MQ=255)
                                                      cTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttt                                   <  1:1017004/93‑1 (MQ=255)
                                                      cTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttt                                   <  1:8771808/93‑1 (MQ=255)
                                                        ggATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttttt                                 >  2:634888/1‑93 (MQ=255)
                                                        ggATAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTAtttttt                                 >  2:6192561/1‑93 (MQ=255)
                                                         gATAGGGCGTTTACGCGCCATCCGGCATAAAAAGCCCCCACTTGGCCAACAATCGGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTa                                <  2:2012762/93‑1 (MQ=37)
                                                           tAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACt                              >  1:5939377/1‑93 (MQ=255)
                                                           tAAGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACt                              <  1:6944610/93‑1 (MQ=255)
                                                             aGGCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTc                            <  1:1003423/93‑1 (MQ=255)
                                                              ggCGTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCa                           >  1:8364059/1‑93 (MQ=255)
                                                                 gTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCa                        <  2:7796049/93‑1 (MQ=255)
                                                                 gTTTACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCa                        >  1:6770320/1‑93 (MQ=255)
                                                                  tttACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCAt                       <  1:8407618/93‑1 (MQ=255)
                                                                    tACGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATaa                     <  2:9806498/93‑1 (MQ=255)
                                                                     aCGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATaaa                    <  2:4744376/93‑1 (MQ=255)
                                                                     aCGCCGCATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATaaa                    <  1:2987780/93‑1 (MQ=255)
                                                                      cgccgcATCCGGCATAAAAAGCGCGCACTTTGTCAACGATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAg                   >  1:5267505/1‑93 (MQ=255)
                                                                      cgccgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAg                   <  1:6879771/93‑1 (MQ=255)
                                                                       gccgcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCCTACAGt                  >  1:6521698/1‑93 (MQ=255)
                                                                          gcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGTAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCATCATAAAGTCAt               >  1:5707494/1‑93 (MQ=255)
                                                                          gcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCAt               <  2:7220056/93‑1 (MQ=255)
                                                                          gcATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCAt               >  1:3161869/1‑93 (MQ=255)
                                                                           cATCTGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATa              <  1:1080361/93‑1 (MQ=255)
                                                                           cATCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATa              >  1:5044810/1‑93 (MQ=255)
                                                                            aTCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATAc             >  2:9189434/1‑93 (MQ=255)
                                                                            aTCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATAc             <  1:2016514/93‑1 (MQ=255)
                                                                             tCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACt            >  2:1840646/1‑93 (MQ=255)
                                                                             tCCGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACt            >  1:4727463/1‑93 (MQ=255)
                                                                               cGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGa          <  1:7267008/93‑1 (MQ=255)
                                                                               cGGCATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGa          <  2:2144147/93‑1 (MQ=255)
                                                                                  cATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATa       <  1:5163151/93‑1 (MQ=255)
                                                                                  cATAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATa       <  1:9858906/93‑1 (MQ=255)
                                                                                    tAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACt     >  2:4416292/1‑93 (MQ=255)
                                                                                    tAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACt     >  2:953114/1‑93 (MQ=255)
                                                                                    tAAAAAGCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACt     >  2:162161/1‑93 (MQ=255)
                                                                                     aaaaagCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTc    >  1:2745877/1‑93 (MQ=255)
                                                                                      aaaagCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTCa   <  2:4985398/93‑1 (MQ=255)
                                                                                       aaagCGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTCAt  >  2:6800401/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCTTATCAGGCCTACCTTTTCCTGCAACACTTTGAATTTATGAGTTTTTGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCATAAAAAACGCGCACTTTGTCAACAATCTGAAACGCCGGAGATTTTTCTCTCCGGCGTTATTTTTTACTTCAGCATAAAGTCATACTGAATACTCAT  >  NC_000913/1689463‑1689642

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