Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,438,341 A→G V151A (GTG→GCG)  ytfL ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,438,3410AG99.7% 1081.6 / ‑5.5 343V151A (GTG→GCG) ytfLputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (173/169);  total (174/169)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 361 total aligned reads displayed.

TAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTC  >  NC_000913/4438251‑4438430
                                                                                          |                                                                                         
tAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGcc                                                                                         >  1:7724988/1‑93 (MQ=255)
tAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGcc                                                                                         <  1:858139/93‑1 (MQ=255)
tAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGcc                                                                                         <  2:1329427/93‑1 (MQ=255)
  gATGTCATCAGAAGTGATGTCATCTTTACGTACCATTTGCAGTTTGAAACTACGGAAGATTCTGTTCGCCAGGCAGTTGAAGAACCACACCAg                                                                                       >  2:2919468/1‑93 (MQ=255)
  gATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAg                                                                                       >  1:3144329/1‑93 (MQ=255)
  gATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAg                                                                                       <  2:4883646/93‑1 (MQ=255)
   aTGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTCCGCCAGGCCGTTGAAGAACCACGCCAGc                                                                                      <  2:88913/93‑1 (MQ=255)
    tGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCg                                                                                     >  1:3084932/1‑93 (MQ=255)
     gTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCgg                                                                                    <  1:6631992/93‑1 (MQ=255)
      tcatcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCggg                                                                                   >  2:2595959/1‑93 (MQ=255)
        atcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGt                                                                                 >  2:4548763/1‑93 (MQ=255)
        atcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGt                                                                                 <  1:4116596/93‑1 (MQ=255)
        atcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTGCGCCAGGCCGTTGAGTAACCACGCCCTCgggtg                                                                                 >  1:1560423/1‑91 (MQ=255)
         tcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTg                                                                                >  2:6522616/1‑93 (MQ=255)
         tcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTg                                                                                <  1:743011/93‑1 (MQ=255)
          caGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGc                                                                               <  1:7499313/93‑1 (MQ=255)
            gAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCaa                                                                             >  1:5423831/1‑93 (MQ=255)
            gAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCaa                                                                             >  1:989516/1‑93 (MQ=255)
             aaGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCaaa                                                                            <  2:7247603/93‑1 (MQ=255)
              aGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAAc                                                                           <  1:4438922/93‑1 (MQ=255)
                 gATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTa                                                                        >  1:5734485/1‑93 (MQ=255)
                  aTGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTAc                                                                       <  1:3632438/93‑1 (MQ=255)
                   tGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACa                                                                      >  1:6962064/1‑93 (MQ=255)
                    gtaaTCTTTACGTACAATGAGCAGTTTGAAAATTCGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAg                                                                     <  2:6232130/90‑1 (MQ=255)
                    gTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAg                                                                     <  1:4909397/93‑1 (MQ=255)
                       aTCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCa                                                                  >  1:5123043/1‑93 (MQ=255)
                        tCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAg                                                                 <  1:6314702/93‑1 (MQ=255)
                        tCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAg                                                                 >  1:2895111/1‑93 (MQ=255)
                        ccTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAg                                                                 >  2:7022127/2‑93 (MQ=255)
                         cTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGa                                                                <  1:67476/93‑1 (MQ=255)
                          tttACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGaa                                                               <  1:2521080/93‑1 (MQ=255)
                          tttACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGaa                                                               >  2:6116697/1‑93 (MQ=255)
                          tttACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGaa                                                               >  1:6176623/1‑93 (MQ=255)
                          tttACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGaa                                                               <  2:5556247/93‑1 (MQ=255)
                          tttACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGaa                                                               <  1:1876437/93‑1 (MQ=255)
                            tACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAgc                                                             >  2:2431776/1‑93 (MQ=255)
                             aCGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAgcg                                                            >  1:5618793/1‑93 (MQ=255)
                               gTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCa                                                          >  2:5135699/1‑93 (MQ=255)
                               gTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCa                                                          >  1:4615796/1‑93 (MQ=255)
                               gTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCa                                                          <  1:7339086/93‑1 (MQ=255)
                               gTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCa                                                          >  2:5790391/1‑93 (MQ=255)
                                tACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCAt                                                         >  2:7478051/1‑93 (MQ=255)
                                 aCCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATc                                                        <  2:5253328/93‑1 (MQ=255)
                                   cATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCgg                                                      <  2:6798331/93‑1 (MQ=255)
                                   cATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCgg                                                      >  1:2071477/1‑93 (MQ=255)
                                     ttGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGt                                                    <  2:7054509/93‑1 (MQ=255)
                                     ttGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGt                                                    >  2:157795/1‑93 (MQ=255)
                                        gCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTtga                                                 <  1:2320356/93‑1 (MQ=255)
                                        gCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTtga                                                 >  1:5806719/1‑93 (MQ=255)
                                           gTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTtgatga                                              <  1:376226/93‑1 (MQ=255)
                                            tttGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTgatgat                                             >  1:1073185/1‑93 (MQ=255)
                                            tttGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTgatgat                                             >  1:3945881/1‑93 (MQ=255)
                                             ttGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATa                                            <  1:4765614/93‑1 (MQ=255)
                                             ttGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATa                                            <  1:6518202/93‑1 (MQ=255)
                                               gAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACg                                          >  2:767553/1‑93 (MQ=255)
                                                aaaaTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGc                                         <  2:687113/93‑1 (MQ=255)
                                                 aaaTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCa                                        >  2:1430842/1‑93 (MQ=255)
                                                 aaaTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCa                                        <  2:4112649/93‑1 (MQ=255)
                                                  aaTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCaa                                       <  2:2310626/93‑1 (MQ=255)
                                                       ggAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCca                                  >  2:6366405/1‑93 (MQ=255)
                                                        gAAGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCcac                                 >  2:4725584/1‑93 (MQ=255)
                                                        gAAGATTATGTCCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCcac                                 >  2:6712943/1‑93 (MQ=255)
                                                         aaGATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAACcaca                                >  1:3422382/1‑93 (MQ=255)
                                                           gATTATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGc                              >  2:3910554/1‑93 (MQ=255)
                                                             ttATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCtt                            >  2:1778087/1‑93 (MQ=255)
                                                             ttATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCtt                            >  2:5385809/1‑93 (MQ=255)
                                                              tATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTc                           <  2:4415213/93‑1 (MQ=255)
                                                              tATGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTc                           <  1:924574/93‑1 (MQ=255)
                                                              tATGTTCGCCAGGCCGTGGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTc                           <  2:4268540/93‑1 (MQ=255)
                                                               aTGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCt                          <  2:3650739/93‑1 (MQ=255)
                                                                tGTTCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTg                         <  2:3438173/93‑1 (MQ=255)
                                                                tGTTCGCCAGGCCGGTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTg                         <  1:1700721/93‑1 (MQ=255)
                                                                  ttCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGgc                       <  1:6896530/93‑1 (MQ=255)
                                                                   tCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGgcg                      <  1:2769506/93‑1 (MQ=255)
                                                                   tCGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGgcg                      <  1:1444344/93‑1 (MQ=255)
                                                                    cGCCAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGgcgc                     >  1:261880/1‑93 (MQ=255)
                                                                      ccAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGc                              <  1:2200674/82‑1 (MQ=255)
                                                                      ccAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGc                              >  2:2200673/1‑82 (MQ=255)
                                                                       cAGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAAt                  <  2:2069268/93‑1 (MQ=255)
                                                                        aGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATc                 >  2:2702129/1‑93 (MQ=255)
                                                                        aGGCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATc                 >  2:2768676/1‑93 (MQ=255)
                                                                         ggCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCa                <  2:1267322/93‑1 (MQ=255)
                                                                         ggCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCa                >  2:2118792/1‑93 (MQ=255)
                                                                         ggCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCa                <  1:4310365/93‑1 (MQ=255)
                                                                          gCCGTTGAAGAACCACGCCAGCGGGGTGCAAGCGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCAt               <  2:7852756/93‑1 (MQ=255)
                                                                          gCCGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCAt               >  2:5869505/1‑93 (MQ=255)
                                                                            cGTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAc             <  2:4310601/93‑1 (MQ=255)
                                                                             gTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAcc            <  2:4194973/93‑1 (MQ=255)
                                                                             gTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAcc            >  2:3916767/1‑93 (MQ=255)
                                                                             gTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAcc            >  1:7821859/1‑93 (MQ=255)
                                                                             gTTGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATAcc            >  2:3059295/1‑93 (MQ=255)
                                                                               tGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGa          <  1:1164435/93‑1 (MQ=255)
                                                                               tGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGa          <  1:6105171/93‑1 (MQ=255)
                                                                               tGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGa          <  2:351412/93‑1 (MQ=255)
                                                                               tGAAGAACCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGa          <  1:2415237/93‑1 (MQ=255)
                                                                                gaagaaCCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGAt         <  1:3292780/93‑1 (MQ=255)
                                                                                gaagaaCCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGAt         <  1:2041585/93‑1 (MQ=255)
                                                                                  agaaCCACGCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGc       <  1:5278818/93‑1 (MQ=255)
                                                                                   gaaCCACGCCAGCGGGGTGCAACCGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCg      <  2:7286586/93‑1 (MQ=255)
                                                                                       cacgCCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTc  <  1:1397470/93‑1 (MQ=255)
                                                                                          |                                                                                         
TAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTC  >  NC_000913/4438251‑4438430

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