Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,523 A→G L239P (CTG→CCG)  yieO ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,5230AG99.6% 1680.7 / ‑4.1 519L239P (CTG→CCG) yieOputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (270/246);  total (270/249)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 534 total aligned reads displayed.

ACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGA  >  NC_000913/3937433‑3937612
                                                                                          |                                                                                         
aCGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGt                                                                                         <  2:354777/93‑1 (MQ=255)
  gATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTaa                                                                                       <  1:9199067/93‑1 (MQ=255)
   aTACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAc                                                                                      >  2:6214889/1‑93 (MQ=255)
   aTACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAc                                                                                      >  1:4741759/1‑93 (MQ=255)
   aTACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAc                                                                                      >  2:3782108/1‑93 (MQ=255)
    tACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAccc                                                                                     >  1:9480618/1‑93 (MQ=255)
    tACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAcc                                                                                     >  2:7151217/1‑93 (MQ=255)
     aCCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAccc                                                                                    >  1:8444589/1‑93 (MQ=255)
     aCCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAAccc                                                                                    >  2:4755077/1‑93 (MQ=255)
       cGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGa                                                                                  <  1:6216919/93‑1 (MQ=255)
        gATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGAt                                                                                 >  2:2221827/1‑93 (MQ=255)
        gATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGAt                                                                                 >  1:6706341/1‑93 (MQ=255)
         aTCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATg                                                                                <  1:4432643/93‑1 (MQ=255)
          tCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGc                                                                               >  1:135490/1‑93 (MQ=255)
            gagaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTg                                                                             >  1:4210148/1‑93 (MQ=255)
             agaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTgg                                                                            <  1:3486193/93‑1 (MQ=255)
             agaAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTgg                                                                            >  2:8061688/1‑93 (MQ=255)
               aaaGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTg                                                                          >  2:3928745/1‑93 (MQ=255)
                   tGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACaa                                                                      >  1:1975100/1‑93 (MQ=255)
                    gCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAAt                                                                     >  1:2476169/1‑93 (MQ=255)
                     cGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAAtt                                                                    >  2:5765199/1‑93 (MQ=255)
                     cGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAAtt                                                                    <  2:9869467/93‑1 (MQ=255)
                      gggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTa                                                                   >  1:4047615/1‑93 (MQ=255)
                       ggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTAc                                                                  <  2:7848132/93‑1 (MQ=255)
                       ggTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTAc                                                                  <  2:1546143/93‑1 (MQ=255)
                         ttttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCg                                                                <  2:6706277/93‑1 (MQ=255)
                         ttttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCg                                                                >  2:8500776/1‑93 (MQ=255)
                          tttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGTGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGt                                                               >  1:5869857/1‑93 (MQ=255)
                          tttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGt                                                               >  2:7244787/1‑93 (MQ=255)
                           ttAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTc                                                              >  1:8925313/1‑93 (MQ=255)
                             aaataaatCTAATGAAATTAATGGTTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCaa                                                            <  1:6064034/93‑1 (MQ=255)
                             aaataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCaa                                                            <  2:6118744/93‑1 (MQ=255)
                             aaataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCaa                                                            <  1:7065902/93‑1 (MQ=255)
                             aaataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCaa                                                            <  1:9763675/93‑1 (MQ=255)
                             aaataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCaa                                                            >  1:2203946/1‑93 (MQ=255)
                              aataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAg                                                           >  1:6013938/1‑93 (MQ=255)
                              aataaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAg                                                           >  1:2739181/1‑93 (MQ=255)
                                 aaatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCa                                                        >  2:9353523/1‑93 (MQ=255)
                                  aatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCaa                                                       <  2:5168058/93‑1 (MQ=255)
                                  aatCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCaa                                                       >  1:7927689/1‑93 (MQ=255)
                                   atCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTAAAGGCAAt                                                      >  1:8737310/1‑93 (MQ=255)
                                   atCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCAAGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAAt                                                      <  2:7967057/93‑1 (MQ=255)
                                     cTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATcc                                                    <  1:299921/93‑1 (MQ=255)
                                     cTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATcc                                                    <  1:1164473/93‑1 (MQ=255)
                                      tAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCa                                                   >  2:5461663/1‑93 (MQ=255)
                                      tAATGAAAATAATGGGTTTGGCGTGCGTCGTGAATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCa                                                   <  2:1204497/93‑1 (MQ=255)
                                          gAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTg                                               <  2:5807216/93‑1 (MQ=255)
                                            aaTTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGc                                             >  2:1495163/1‑93 (MQ=255)
                                            aaTTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGc                                             >  2:3346979/1‑93 (MQ=255)
                                             aTTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGAAATCCAGCTGGCg                                            >  1:5019786/1‑93 (MQ=255)
                                              ttAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGa                                           >  1:9051100/1‑93 (MQ=255)
                                               tAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGAGACAATTACCGTCAAGGCAATCCAGCTGGCGAc                                          >  1:9676676/1‑93 (MQ=255)
                                                aaTGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACa                                         >  1:7683647/1‑93 (MQ=255)
                                                aaTGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACa                                         <  2:9704976/93‑1 (MQ=255)
                                                aaTGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACa                                         <  1:7393803/93‑1 (MQ=255)
                                                 aTGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACaa                                        <  2:2108788/93‑1 (MQ=255)
                                                 aTGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACaa                                        >  1:5576313/1‑93 (MQ=255)
                                                  tGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAAt                                       >  2:825139/1‑93 (MQ=255)
                                                     gTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCtt                                    >  1:6286552/1‑93 (MQ=255)
                                                     gTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCtt                                    >  1:3297614/1‑93 (MQ=255)
                                                      tttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttt                                   <  2:5668015/93‑1 (MQ=255)
                                                      tttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttt                                   <  1:10017660/93‑1 (MQ=255)
                                                      tttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttt                                   <  2:8679428/93‑1 (MQ=255)
                                                      tttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCttt                                   >  2:9032220/1‑93 (MQ=255)
                                                       ttGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCtttt                                  <  1:369928/93‑1 (MQ=255)
                                                        tGGCGTGCTTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTc                                 >  1:5715941/1‑93 (MQ=255)
                                                        tGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTc                                 >  1:5199787/1‑93 (MQ=255)
                                                         ggCGTGCGTCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTcc                                <  1:1679459/93‑1 (MQ=255)
                                                              gcgtCGTGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCTTTTCCCCGaa                           >  2:2841311/1‑93 (MQ=255)
                                                               cgtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCTTTTCCCCGAAt                          <  2:1855653/93‑1 (MQ=255)
                                                               cgtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAAt                          <  2:2483305/93‑1 (MQ=255)
                                                                gtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCTTTTCCCCGAATa                         >  2:490431/1‑93 (MQ=255)
                                                                gtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCTTTTCCCCGAATa                         <  1:4834864/93‑1 (MQ=255)
                                                                gtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATa                         <  1:2834020/93‑1 (MQ=255)
                                                                gtcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATa                         <  2:9218332/93‑1 (MQ=255)
                                                                 tcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAg                        >  1:6551510/1‑93 (MQ=255)
                                                                 tcgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAg                        <  2:6999856/93‑1 (MQ=255)
                                                                  cgtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGc                       >  2:1004959/1‑93 (MQ=255)
                                                                   gtGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGct                      >  2:7687859/1‑93 (MQ=255)
                                                                    tGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCTTTTCCCCGAATAGctc                     >  1:1547785/1‑93 (MQ=255)
                                                                    tGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGctc                     >  1:1480582/1‑93 (MQ=255)
                                                                    tGCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGctc                     <  2:2623811/93‑1 (MQ=255)
                                                                     gCATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGctct                    >  2:5025691/1‑93 (MQ=255)
                                                                      cATGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTa                   <  2:9475147/93‑1 (MQ=255)
                                                                        tGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                 >  2:451044/1‑93 (MQ=255)
                                                                        tGGAGAATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTAtt                 <  1:6992972/93‑1 (MQ=255)
                                                                          gagaATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTcc               >  2:7412596/1‑93 (MQ=255)
                                                                          gagaAGATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTcc               <  2:3279563/93‑1 (MQ=255)
                                                                            gaATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGc             >  1:4602221/1‑93 (MQ=255)
                                                                            gaATATAGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGc             >  2:41893/1‑93 (MQ=255)
                                                                               tataGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGTAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTg          <  1:7517519/93‑1 (MQ=255)
                                                                               tataGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTg          >  1:5091710/1‑93 (MQ=255)
                                                                               tataGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTg          <  2:3187022/93‑1 (MQ=255)
                                                                               tataGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTg          >  1:5348897/1‑93 (MQ=255)
                                                                                 taGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgag        <  2:6374751/93‑1 (MQ=255)
                                                                                 taGAGAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgag        >  2:8820145/1‑93 (MQ=255)
                                                                                  agagAAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgaga       <  2:8387206/93‑1 (MQ=255)
                                                                                   gagaAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaa      <  2:5392161/93‑1 (MQ=255)
                                                                                    agaAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaag     <  1:8156301/93‑1 (MQ=255)
                                                                                     gaAGCGGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTgagaaga    >  2:5968338/1‑93 (MQ=255)
                                                                                          |                                                                                         
ACGATACCGATCGAGAAAGTGCGGGTTTTAAATAAATCTAATGAAATTAATGGGTTTGGCGTGCGTCGTGCATGGAGAATATAGAGAAGCAGTAACCCGATGCTGGTGACAATTACCGTCAAGGCAATCCAGCTGGCGACAATCTTTTCCCCGAATAGCTCTATTCCGCTTGAGAAGAGA  >  NC_000913/3937433‑3937612

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