Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 903,254 A→G V146A (GTC→GCC)  ybjP ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913903,2540AG100.0% 730.8 / NA 226V146A (GTC→GCC) ybjPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (96/130);  total (96/130)

Only 100 of 235 total aligned reads displayed.

TTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACGG  >  NC_000913/903163‑903343
                                                                                            |                                                                                         
  gATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCa                                                                                       >  1:5848468/1‑93 (MQ=255)
    tCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAac                                                                                     >  2:3853537/1‑93 (MQ=255)
      aaGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAacac                                                                                   >  2:1097603/1‑93 (MQ=255)
      aaGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACCCTACCGCCAAGTAACGCACATCGCCGCATGCCCAacac                                                                                   >  2:142986/1‑93 (MQ=255)
      aaGGCGTGAGCGATTCTCGATGGACTGACGGAGCTTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAacac                                                                                   <  2:2382970/93‑1 (MQ=255)
       aGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACt                                                                                  >  1:4208360/1‑93 (MQ=255)
       aGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACt                                                                                  >  2:3025805/1‑93 (MQ=255)
       aGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACt                                                                                  >  2:653021/1‑93 (MQ=255)
       aGGCGTTAGCGATTCTCGAGGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACt                                                                                  <  2:3683718/93‑1 (MQ=255)
        ggCGTTAGTGATCGGCGAGGGACTGCCGGATCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTg                                                                                 <  2:3582122/93‑1 (MQ=255)
        ggCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTg                                                                                 >  2:2743346/1‑93 (MQ=255)
        ggCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTg                                                                                 >  2:4419684/1‑93 (MQ=255)
         gctttAGCGATTCTCGAGTGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTgg                                                                                <  2:3370871/90‑1 (MQ=255)
           gTTAGCTATCCTCGATGGACTGGCGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGcc                                                                              <  1:110249/93‑1 (MQ=255)
            ttAGGGATTCTCGATGGAGTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                                             <  1:653382/93‑1 (MQ=255)
            ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                                             >  2:3457342/1‑93 (MQ=255)
            ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                                             >  2:457017/1‑93 (MQ=255)
            tgaGCGATTCTCGATGGACGGACGGAGCGTCCCCGCCGGGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                                             <  2:6002056/91‑1 (MQ=255)
              aGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTc                                                                           >  1:2170143/1‑93 (MQ=255)
               gCGATTCTCGATGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                                          <  1:3615348/93‑1 (MQ=255)
               gCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                                          >  2:4843870/1‑93 (MQ=255)
               gCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                                          <  1:3831259/93‑1 (MQ=255)
               gCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                                          >  1:353037/1‑93 (MQ=255)
               gCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                                          >  2:783165/1‑93 (MQ=255)
               gCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                                          <  1:2983551/93‑1 (MQ=255)
                cGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCt                                                                         >  2:4671777/1‑93 (MQ=255)
                 gATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCtt                                                                        >  1:4875034/1‑92 (MQ=255)
                 gATTCTCGAGGGACTGAGGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTg                                                                        <  1:595025/93‑1 (MQ=255)
                  aTTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGa                                                                       >  1:5604173/1‑93 (MQ=255)
                  aTTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGa                                                                       >  2:5716120/1‑93 (MQ=255)
                   ttCTCGATGGAGTGACGGAGCTTACCCGCTGTGGCATGAACGCTCCCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                                                      <  1:694292/93‑1 (MQ=255)
                   ttCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                                                      <  1:5640182/93‑1 (MQ=255)
                      tcGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                                                   <  1:6277327/93‑1 (MQ=255)
                      tcGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                                                   >  2:1649671/1‑93 (MQ=255)
                      tcGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                                                   <  2:4353736/93‑1 (MQ=255)
                           ggACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAac                                                              >  2:250189/1‑93 (MQ=255)
                            gACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAaca                                                             >  2:2248876/1‑93 (MQ=255)
                              cTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACt                                                           <  2:1057232/93‑1 (MQ=255)
                              cTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACt                                                           >  2:3252539/1‑93 (MQ=255)
                                 aCGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCa                                                        <  1:2372354/93‑1 (MQ=255)
                                 aCGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCa                                                        <  2:5597654/93‑1 (MQ=255)
                                   ggAGCGTGCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATc                                                      <  2:1515979/93‑1 (MQ=255)
                                   ggAGCGTCCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATc                                                      <  1:4047841/93‑1 (MQ=255)
                                    gAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATcc                                                     >  2:4366410/1‑93 (MQ=255)
                                    gAGCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATcc                                                     <  2:264028/93‑1 (MQ=255)
                                      gCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTg                                                   >  1:4353148/1‑93 (MQ=255)
                                      gCGTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTg                                                   >  1:2772570/1‑93 (MQ=255)
                                        gtccccGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGcc                                                 <  2:2408535/90‑1 (MQ=255)
                                        gTACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGcc                                                 >  1:3863856/1‑93 (MQ=255)
                                         tACCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCa                                                >  2:5813005/1‑93 (MQ=255)
                                          aCCCGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCaa                                               >  1:5020860/1‑93 (MQ=255)
                                           cccGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCTCAGCCCTGGCCTTCCTGAACCATTAACACTTCATCCTGCCacc                                              >  1:5745883/1‑91 (MQ=255)
                                           cccGCCGTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAc                                              <  2:539701/93‑1 (MQ=255)
                                             cgccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAccc                                            >  1:5447866/1‑93 (MQ=255)
                                              gccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCt                                           <  1:5825809/93‑1 (MQ=255)
                                              gccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCt                                           >  1:1458319/1‑93 (MQ=255)
                                               ccgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTg                                          >  2:6346830/1‑93 (MQ=255)
                                                cgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTgg                                         <  2:5533862/93‑1 (MQ=255)
                                                cgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTgg                                         >  2:5991489/1‑93 (MQ=255)
                                                cgTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTgg                                         >  1:2030242/1‑93 (MQ=255)
                                                catGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTgg                                         <  1:1898351/91‑1 (MQ=255)
                                                 gTGGCATGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGt                                        >  2:905785/1‑93 (MQ=255)
                                                      aTGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCAccc                                   <  1:5041721/93‑1 (MQ=255)
                                                      aTGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCAccc                                   >  1:2932137/1‑93 (MQ=255)
                                                       tcaaCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCt                                  <  2:3966394/91‑1 (MQ=255)
                                                       tGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGTCAACCCTGGTCACCCt                                  >  1:1770124/1‑93 (MQ=255)
                                                       tGAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCt                                  <  1:5293981/93‑1 (MQ=255)
                                                        gAACGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTg                                 >  2:3300394/1‑93 (MQ=255)
                                                         aaCGCTCCCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACACTGGTCACCCTGt                                <  2:3735548/93‑1 (MQ=255)
                                                         aaCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt                                >  1:139614/1‑93 (MQ=255)
                                                         aaCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt                                <  1:3691147/93‑1 (MQ=255)
                                                         aaCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt                                <  2:1808435/93‑1 (MQ=255)
                                                          aCGCTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGtt                               <  2:1701147/93‑1 (MQ=255)
                                                             cTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTga                            <  2:3911385/93‑1 (MQ=255)
                                                             cTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTga                            <  2:2914483/93‑1 (MQ=255)
                                                             cTACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTga                            <  1:790693/93‑1 (MQ=255)
                                                              tACCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTgag                           <  2:4133/93‑1 (MQ=255)
                                                               aCCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTgaga                          <  2:3649062/93‑1 (MQ=255)
                                                               aCCGCCAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTgaga                          <  2:6365413/93‑1 (MQ=255)
                                                                   gccaaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATcc                       <  1:2226550/90‑1 (MQ=255)
                                                                    ccaaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCa                      >  1:95385/3‑93 (MQ=255)
                                                                     caaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAc                     >  1:6587587/2‑93 (MQ=255)
                                                                      aaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAcg                    >  2:2694690/1‑93 (MQ=255)
                                                                       aGTATCTCACATCGTCGATGTCCCAACACTGGCTTTCCTGAATCATTAACACTTCATCCTGCCACCCCTTGTCACCCTGTTTGAGATCCAcgc                   <  2:1542750/93‑1 (MQ=255)
                                                                          aaCGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCa                >  1:5513421/1‑93 (MQ=255)
                                                                            cGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAAt              <  2:5238889/93‑1 (MQ=255)
                                                                             gCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATg             <  2:6649557/93‑1 (MQ=255)
                                                                             gCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATg             <  1:5708513/93‑1 (MQ=255)
                                                                             gCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATg             >  2:6270497/1‑93 (MQ=255)
                                                                               acaTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGa           >  1:378563/1‑93 (MQ=255)
                                                                               acaTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCTGGTCACCCTGTTTGAGATCCACGCGCAATGGaa          >  1:5172249/1‑93 (MQ=255)
                                                                                 aTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAAt         >  2:5503006/1‑93 (MQ=255)
                                                                                  tcgtcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATg        >  1:5388507/1‑93 (MQ=255)
                                                                                  tcgtcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATg        >  2:405065/1‑93 (MQ=255)
                                                                                     tcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTa     <  1:3117200/93‑1 (MQ=255)
                                                                                     tcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTa     >  2:3087287/1‑93 (MQ=255)
                                                                                      cgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTAc    <  2:836623/93‑1 (MQ=255)
                                                                                      cgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTAc    <  1:4557411/93‑1 (MQ=255)
                                                                                       gATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACg   >  2:4936444/1‑93 (MQ=255)
                                                                                         aTGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACgg  <  2:6023033/93‑1 (MQ=255)
                                                                                            |                                                                                         
TTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACGG  >  NC_000913/903163‑903343

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