Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 778,306 T→C G189G (GGT→GGC)  tolB → Tol‑Pal system periplasmic protein TolB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913778,3060TC99.7% 1342.2 / ‑4.4 328G189G (GGT→GGC) tolBTol‑Pal system periplasmic protein TolB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (179/147);  total (180/148)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

Only 100 of 328 total aligned reads displayed.

ACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCG  >  NC_000913/778208‑778403
                                                                                                  |                                                                                                 
aCCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTa                                                                                                 >  1:1455954/1‑101 (MQ=255)
  cGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACa                                                                                               >  1:1455958/1‑101 (MQ=255)
   ggCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACaa                                                                                              <  2:1455432/101‑1 (MQ=255)
    gCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAAc                                                                                             <  2:1455381/101‑1 (MQ=255)
    gCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAAc                                                                                             >  2:1455963/1‑101 (MQ=255)
      aTTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCa                                                                                           >  2:1455968/1‑101 (MQ=255)
       ttAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAg                                                                                          <  1:1455506/101‑1 (MQ=255)
       ttAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAg                                                                                          >  2:1455972/1‑101 (MQ=255)
        tAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGAGCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGt                                                                                         >  1:1455973/1‑101 (MQ=255)
         aaaGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGtt                                                                                        >  2:1455980/1‑101 (MQ=255)
         aaaGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGtt                                                                                        >  1:1455979/1‑101 (MQ=255)
          aaGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGttt                                                                                       <  1:1455644/101‑1 (MQ=255)
          aaGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGttt                                                                                       <  1:1455337/101‑1 (MQ=255)
           aGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTg                                                                                      <  1:1455540/101‑1 (MQ=255)
           aGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTg                                                                                      >  1:1455983/1‑101 (MQ=255)
             gTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTc                                                                                    <  1:1455632/101‑1 (MQ=255)
                cGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgtt                                                                                 >  2:1455988/1‑101 (MQ=255)
                cGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgtt                                                                                 <  2:1455577/101‑1 (MQ=255)
                 gTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTa                                                                                                 >  1:1455992/1‑84 (MQ=255)
                   tCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttcac                                                                              >  2:1455993/1‑101 (MQ=255)
                     cGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttcaccg                                                                            >  2:1455998/1‑101 (MQ=255)
                      gTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTcgttcaccgt                                                                           <  1:1455902/101‑1 (MQ=255)
                      gTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttcaccgt                                                                           >  1:1455999/1‑101 (MQ=255)
                        aCTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttcaccgttc                                                                         <  2:1455487/101‑1 (MQ=255)
                        aCTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttcaccgttc                                                                         >  2:1456004/1‑101 (MQ=255)
                        aCTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttcaccgttc                                                                         >  1:1456006/1‑101 (MQ=255)
                         cTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttcaccgttca                                                                        <  2:1455520/101‑1 (MQ=255)
                           cGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCgttcaccgttcacc                                                                      <  2:1455482/101‑1 (MQ=255)
                            gTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGttcaccgttcaccg                                                                     >  1:1456009/1‑101 (MQ=255)
                             tATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGc                                                                    >  2:1456016/1‑101 (MQ=255)
                             tATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGc                                                                    >  1:1456014/1‑101 (MQ=255)
                             tATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGc                                                                    >  1:1456018/1‑101 (MQ=255)
                              aTTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCCTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCa                                                                   <  2:1455639/101‑1 (MQ=255)
                               ttGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAg                                                                  >  2:1456022/1‑101 (MQ=255)
                                 gCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGcc                                                                >  1:1456027/1‑101 (MQ=255)
                                 gCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGcc                                                                >  2:1456028/1‑101 (MQ=255)
                                    tACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCt                                                             <  2:1455757/101‑1 (MQ=255)
                                    tACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCt                                                             >  1:1456031/1‑101 (MQ=255)
                                    tACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCt                                                             <  1:1455649/101‑1 (MQ=255)
                                      cGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGa                                                           <  1:1455486/101‑1 (MQ=255)
                                      cGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGa                                                           <  2:1455401/101‑1 (MQ=255)
                                        ttgttCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATg                                                         <  2:1455508/101‑1 (MQ=255)
                                         tgttCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGt                                                        <  1:1455406/101‑1 (MQ=255)
                                         tgttCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGt                                                        >  2:1456038/1‑101 (MQ=255)
                                          gttCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGcc                                                                >  1:1456040/1‑92 (MQ=255)
                                            tCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGtctc                                                     <  1:1455128/101‑1 (MQ=255)
                                            tCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGtctc                                                     <  2:1455606/101‑1 (MQ=255)
                                            tCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGtctc                                                     >  1:1456044/1‑101 (MQ=255)
                                              aGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCg                                                   <  2:1455596/101‑1 (MQ=255)
                                              aGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCg                                                   >  1:1456047/1‑101 (MQ=255)
                                               gACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCgg                                                  >  1:1456051/1‑101 (MQ=255)
                                                  cAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGttt                                                                                       >  2:1456054/1‑61 (MQ=255)
                                                   aaCGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTg                                              >  1:1456055/1‑101 (MQ=255)
                                                    aCGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTgg                                             <  2:1455417/101‑1 (MQ=255)
                                                    aCGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTgg                                             <  2:1455285/101‑1 (MQ=255)
                                                      ggcggTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTcgttca                                                                               <  1:1455771/65‑1 (MQ=255)
                                                      ggcggTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTc                                           >  2:1456058/1‑101 (MQ=255)
                                                        cggTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCAc                                         >  2:1456063/1‑101 (MQ=255)
                                                           tCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAg                                      >  1:1456068/1‑101 (MQ=255)
                                                           tCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAg                                      >  2:1456067/1‑101 (MQ=255)
                                                           tCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAg                                      <  2:1455464/101‑1 (MQ=255)
                                                             aGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGAc                                    <  2:1455661/101‑1 (MQ=255)
                                                                tCCCGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGt                                 >  2:1456072/1‑101 (MQ=255)
                                                                 cccGTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGtt                                >  2:1456075/1‑101 (MQ=255)
                                                                   cGTATGAACTGAGCGTATCTGACTATGACGGATACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGAt                                                          >  2:1456079/1‑73 (MQ=255)
                                                                    gTATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGc                                                 >  1:1456083/1‑81 (MQ=255)
                                                                     tATGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCt                              >  2:1456084/1‑99 (MQ=255)
                                                                      aTGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTaaa                           <  1:1455459/101‑1 (MQ=255)
                                                                      aTGAACTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTaaa                           >  1:1456089/1‑101 (MQ=255)
                                                                         aaCTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTg                        >  2:1456094/1‑101 (MQ=255)
                                                                         aaCTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTg                        <  1:1455605/101‑1 (MQ=255)
                                                                         aaCTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTg                        >  1:1456093/1‑101 (MQ=255)
                                                                          aCTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGAc                                    >  1:1456098/1‑88 (MQ=255)
                                                                          aCTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTgg                       <  2:1455446/101‑1 (MQ=255)
                                                                           cTGCGCGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGc                      <  2:1455919/101‑1 (MQ=255)
                                                                             gcgcGTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCtt                    >  1:1456103/1‑101 (MQ=255)
                                                                              cgcgTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCgg                                                  >  2:1456105/1‑70 (MQ=255)
                                                                                cgTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtg                 <  2:1455686/101‑1 (MQ=255)
                                                                                 gTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtgt                >  2:1456110/1‑101 (MQ=255)
                                                                                 gTATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtgt                <  1:1455663/101‑1 (MQ=255)
                                                                                  tATCTGACTATGACGGCTACCACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTg                        >  2:1456113/1‑92 (MQ=255)
                                                                                  tATCTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTAtgtg               <  2:1455725/101‑1 (MQ=255)
                                                                                     cTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGAcc            <  2:1455535/101‑1 (MQ=255)
                                                                                     cTGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGAcc            >  2:1456116/1‑101 (MQ=255)
                                                                                      tGACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCt           <  1:1455722/101‑1 (MQ=255)
                                                                                       gACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCtt          >  1:1456118/1‑101 (MQ=255)
                                                                                       gACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCtt          <  1:1455691/101‑1 (MQ=255)
                                                                                       gACTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCtt          <  1:1455528/101‑1 (MQ=255)
                                                                                        aCTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTc         >  2:1456121/1‑101 (MQ=255)
                                                                                         cTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTc         <  2:1455500/100‑1 (MQ=255)
                                                                                         cTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCg        <  1:1455215/101‑1 (MQ=255)
                                                                                         cTATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCg        <  2:1455799/101‑1 (MQ=255)
                                                                                          tATGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGa       <  1:1455423/101‑1 (MQ=255)
                                                                                           aTGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGaa      >  2:1456126/1‑101 (MQ=255)
                                                                                            tGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGaaa     <  2:1455556/101‑1 (MQ=255)
                                                                                            tGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGaaa     >  1:1456128/1‑101 (MQ=255)
                                                                                            tGACGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGaaa     >  2:1456130/1‑101 (MQ=255)
                                                                                              aCGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGc   <  2:1455620/101‑1 (MQ=255)
                                                                                              aCGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGc   <  1:1455281/101‑1 (MQ=255)
                                                                                               cGGCTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCg  >  1:1456135/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
ACCGGCATTAAAGGTGCGTTCCGTACTCGTATTGCCTACGTTGTTCAGACCAACGGCGGTCAGTTCCCGTATGAACTGCGCGTATCTGACTATGACGGTTACAACCAGTTTGTCGTTCACCGTTCACCGCAGCCGCTGATGTCTCCGGCGTGGTCACCAGACGGTTCTAAACTGGCTTATGTGACCTTCGAAAGCG  >  NC_000913/778208‑778403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑