Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 611,741 A→G intergenic (‑24/‑297) fepA ← / → fes FepA/enterochelin esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913611,7410AG88.7% 494.5 / 49.2 194intergenic (‑24/‑297)fepA/fesFepA/enterochelin esterase
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base G (85/87);  total (97/97)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 196 total aligned reads displayed.

TTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTG  >  NC_000913/611649‑611832
                                                                                            |                                                                                           
ttGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCa                                                                                             <  2:1837883/93‑1 (MQ=255)
    gCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGaa                                                                                         <  2:407534/93‑1 (MQ=255)
    gCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGaa                                                                                         <  2:5305498/93‑1 (MQ=255)
    gCCTGCGCTACCCCATAAAGCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGaa                                                                                         <  1:741062/93‑1 (MQ=255)
     ccTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAAt                                                                                        <  1:7119594/93‑1 (MQ=255)
      cTGCGCTACCCCATGAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAAtt                                                                                       >  2:5183523/1‑93 (MQ=255)
      cTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAAtt                                                                                       >  1:5857604/1‑93 (MQ=255)
       tGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTg                                                                                      <  1:5487075/93‑1 (MQ=255)
        gcgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGc                                                                                     <  1:4534741/93‑1 (MQ=255)
         cgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCa                                                                                    <  2:6132373/93‑1 (MQ=255)
         cgcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCa                                                                                    >  2:1922581/1‑93 (MQ=255)
          gcTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCaa                                                                                   >  1:4849766/1‑93 (MQ=255)
           cTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAAc                                                                                  <  2:3469613/93‑1 (MQ=255)
            tACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACt                                                                                 <  1:3317336/93‑1 (MQ=255)
             aCCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTg                                                                                <  1:1723197/93‑1 (MQ=255)
             aCCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTg                                                                                >  2:3240481/1‑93 (MQ=255)
             aCCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTg                                                                                <  1:7211309/93‑1 (MQ=255)
               cccATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTc                                                                              >  2:6248984/1‑93 (MQ=255)
               cccATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTc                                                                              <  1:936279/93‑1 (MQ=255)
                ccATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCg                                                                             >  2:4693035/1‑93 (MQ=255)
                 cATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCgg                                                                            <  1:3373982/93‑1 (MQ=255)
                  aTAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCggg                                                                           >  2:4411297/1‑93 (MQ=255)
                   tAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGc                                                                          >  1:5117224/1‑93 (MQ=255)
                    aaaTCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCa                                                                         <  1:6134259/93‑1 (MQ=255)
                     aaTCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCAt                                                                        >  2:1014979/1‑93 (MQ=255)
                     aaTCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCAt                                                                        <  2:6070119/93‑1 (MQ=255)
                        ccccAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGgt                                                                     >  1:5611287/1‑93 (MQ=255)
                           cAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGgtcgt                                                                  <  1:1088953/93‑1 (MQ=255)
                             gATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGtca                                                                >  1:986453/1‑93 (MQ=255)
                                tGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGtcatca                                                             >  2:913291/1‑93 (MQ=255)
                                 gACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCaa                                                            >  2:4179259/1‑93 (MQ=255)
                                 gACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCaa                                                            <  2:3739297/93‑1 (MQ=255)
                                   ccAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                          <  1:2569288/93‑1 (MQ=255)
                                   ccAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAaca                                                          <  1:6874290/93‑1 (MQ=255)
                                    caacaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAacac                                                         >  2:7121394/1‑93 (MQ=255)
                                      acaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACacga                                                       <  2:2600849/93‑1 (MQ=255)
                                      acaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACacga                                                       >  1:2267207/1‑93 (MQ=255)
                                      acaaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACacga                                                       >  1:6475428/1‑93 (MQ=255)
                                       caaGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACacgac                                                      >  2:3437880/1‑93 (MQ=255)
                                         aGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGc                                                    >  1:3133111/1‑93 (MQ=255)
                                         aGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGc                                                    <  1:4410722/93‑1 (MQ=255)
                                          ggCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                   <  2:3620509/93‑1 (MQ=255)
                                          ggCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                   <  2:1442812/93‑1 (MQ=255)
                                          ggCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                   <  1:4980491/93‑1 (MQ=255)
                                          ggCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTACTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCa                                                   >  2:5006057/1‑93 (MQ=255)
                                             cAGGGAATGAATCTTTTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATcc                                                >  2:5695529/1‑93 (MQ=255)
                                             cAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATcc                                                <  2:4846811/93‑1 (MQ=255)
                                             cAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATcc                                                <  1:1447727/93‑1 (MQ=255)
                                             cAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATcc                                                >  2:1110482/1‑93 (MQ=255)
                                               gggAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCg                                              >  1:1224039/1‑93 (MQ=255)
                                                ggAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGc                                             >  2:1234933/1‑93 (MQ=255)
                                                ggAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGc                                             <  1:5720428/93‑1 (MQ=255)
                                                ggAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGc                                             >  1:6011160/1‑93 (MQ=255)
                                                 gaatgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCt                                            >  2:6606397/1‑93 (MQ=255)
                                                  aatgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTa                                           >  2:790042/1‑93 (MQ=255)
                                                   atgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAc                                          >  1:3606039/1‑93 (MQ=255)
                                                   atgaatCTTATTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAc                                          <  1:1883671/93‑1 (MQ=255)
                                                    tgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                         >  2:6018887/1‑93 (MQ=255)
                                                    tgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                         >  2:5563506/1‑93 (MQ=255)
                                                    tgaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTAcc                                         >  2:3181094/1‑93 (MQ=255)
                                                     gaatCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACcg                                        >  1:1629478/1‑93 (MQ=255)
                                                     gaatCTTCCTGTTCATTGTTTTATTTCTGCATTTTTTCCGCGAATTTCAACTTTCGGGCATGGGCGTCATCAAAACGACCCATCCCGCTAccc                                        >  1:2383364/1‑92 (MQ=255)
                                                           tcttGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGaaaa                                  >  1:1064457/1‑93 (MQ=255)
                                                             ttGTTCATTTTTTTATTCCTGCATTTTCGCCGCGAATTGCAACTGTCGGGCATGATCAACACCAACACGACGCAACCCTCCACCGCGAAAAcc                                <  1:4626099/93‑1 (MQ=255)
                                                             ttGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAAcc                                >  1:4260055/1‑93 (MQ=255)
                                                              tGTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCt                               <  1:7254891/93‑1 (MQ=255)
                                                               gTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                              <  1:5703112/93‑1 (MQ=255)
                                                               gTTCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCtt                              >  1:4016922/1‑93 (MQ=255)
                                                                 tCATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTg                            >  1:6505517/1‑93 (MQ=255)
                                                                  cATTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGa                           <  2:7241561/93‑1 (MQ=255)
                                                                   aTTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGAt                          <  2:2666509/93‑1 (MQ=255)
                                                                   aTTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGAt                          <  1:5137340/93‑1 (MQ=255)
                                                                   aTTGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGAt                          >  2:2071677/1‑93 (MQ=255)
                                                                     tGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATcc                        <  1:302603/93‑1 (MQ=255)
                                                                     tGTTTTATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATcc                        <  1:5894192/93‑1 (MQ=255)
                                                                       ttttATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTg                      >  2:1176524/1‑93 (MQ=255)
                                                                        tttATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGa                     >  1:528804/1‑93 (MQ=255)
                                                                        tttATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGa                     >  2:1675619/1‑93 (MQ=255)
                                                                         ttATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaa                    >  2:4217101/1‑93 (MQ=255)
                                                                         ttATTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGaa                    >  2:5107017/1‑93 (MQ=255)
                                                                           aTTCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                  <  1:1829305/93‑1 (MQ=255)
                                                                           aTTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAg                  <  2:1855653/93‑1 (MQ=255)
                                                                            ttCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGa                 >  2:3890840/1‑93 (MQ=255)
                                                                            ttCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGa                 <  1:7380725/93‑1 (MQ=255)
                                                                            ttCCTGCATTTTTGCCGCGAATTGCAACTGTCGGGAATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGCTCCTGAAAGa                 >  2:2503128/1‑93 (MQ=255)
                                                                            ttCCCGTATTTTGGCCGCGAATTCCAATTGTCGGGCATGTTGGTATTCACCACGACGCACCCCTTTCCGCGAAAACCTTTGATCCTGAAAGac                <  2:188658/93‑1 (MQ=255)
                                                                                tGCATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg             >  1:3552545/1‑93 (MQ=255)
                                                                                tGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACg             >  1:786077/1‑93 (MQ=255)
                                                                                  cATTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgca           <  1:5374621/93‑1 (MQ=255)
                                                                                   aTTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcag          >  2:2927648/1‑93 (MQ=255)
                                                                                   aTTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcag          <  2:224172/93‑1 (MQ=255)
                                                                                   aTTTTTGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcag          >  1:7521090/1‑93 (MQ=255)
                                                                                    tttttGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt         >  2:5936120/1‑93 (MQ=255)
                                                                                    tttttGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt         <  1:4071259/93‑1 (MQ=255)
                                                                                    tttttGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt         >  1:5733221/1‑93 (MQ=255)
                                                                                    tttttGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt         <  2:4615907/93‑1 (MQ=255)
                                                                                    tttttGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagt         <  2:5853364/93‑1 (MQ=255)
                                                                                     ttttGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtg        >  1:6784542/1‑93 (MQ=255)
                                                                                       ttGCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtgca      <  2:770425/93‑1 (MQ=255)
                                                                                         gCCGCGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACgcagtgcagt    >  2:3931879/1‑93 (MQ=255)
                                                                                            |                                                                                           
TTGTGCCTGCGCTACCCCATAAATCCCCAGATTGACCAACAAGGCCAGGGAATGAATCTTCTTGTTCATTGTTTTATTCCTGCATTTTTGCCACGAATTGCAACTGTCGGGCATGGTCGTCATCAACACGACGCATCCCGCTACCGCGAAAACCTTTGATCCTGAAAGACACGCAGTGCAGTTG  >  NC_000913/611649‑611832

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