Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 887,392 T→A 17.8% intergenic (‑53/‑31) ybjJ ← / → rcdA putative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913887,3920TA17.8% 341.5 / 64.9 146intergenic (‑53/‑31)ybjJ/rcdAputative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator
Reads supporting (aligned to +/- strand):  ref base T (67/53);  new base A (14/12);  total (81/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

Only 100 of 152 total aligned reads displayed.

GCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAC  >  NC_000913/887293‑887489
                                                                                                   |                                                                                                 
gcCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAt                                                                                                   >  1:1344366/1‑100 (MQ=255)
 cccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGGAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                  >  1:2698777/1‑100 (MQ=255)
  ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAttt                                                                                                 >  2:183209/1‑100 (MQ=255)
  ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAttt                                                                                                 >  1:2483034/1‑100 (MQ=255)
   cAGGTTCGGCGTTTCAATGGAATAAGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACTACATTAGCAAGGTTGTGTACAt                                                                                                   <  2:2215544/97‑1 (MQ=255)
   cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGGTGTGTACATTTg                                                                                                >  1:3625409/1‑100 (MQ=255)
     ggTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                              <  1:3276601/100‑1 (MQ=255)
        tCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTacac                                                                                           >  2:2611087/1‑100 (MQ=255)
        tCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTacac                                                                                           <  1:3491141/100‑1 (MQ=255)
          ggCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaa                                                                                         <  1:3689052/100‑1 (MQ=255)
           gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                        <  2:679600/100‑2 (MQ=255)
           gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                        <  2:425743/100‑2 (MQ=255)
           gCGTTTCAATGCATTACGTAAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGTCAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                        <  1:989048/100‑2 (MQ=255)
            cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                       <  1:2160624/100‑3 (MQ=255)
             gTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacct                                                                                      <  2:574035/100‑4 (MQ=255)
             gTTTCAATGCATTACGTGAAAAATTTACGGTCATAGATCACGTCAAAAGAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                        <  2:570502/98‑2 (MQ=255)
              tttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTg                                                                                      >  1:105942/1‑99 (MQ=255)
               ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCa                                                                                    <  2:733133/100‑1 (MQ=255)
                tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                   >  2:2896018/1‑100 (MQ=255)
                 cAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCaga                                                                                  <  1:2814725/100‑1 (MQ=255)
                   aTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagaga                                                                                <  1:540422/100‑1 (MQ=255)
                    tGCATTACGTTAAGAGTTTTCGGTCAGATATCACGTCAAAATAAGAAGAGGTAAGACAACATTAGCAAGGGTGTGTACAAATGTACACaacctttc                                                                                   <  2:1315043/96‑7 (MQ=255)
                    attcgcccatGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                               <  1:837438/91‑1 (MQ=255)
                     gCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                              >  1:3210907/1‑100 (MQ=255)
                      cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                             <  2:1767789/100‑1 (MQ=255)
                       aTTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                             >  2:1083779/1‑99 (MQ=255)
                       aTTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                            >  1:2469710/1‑100 (MQ=255)
                           cGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAAt                                                                        >  1:815300/1‑100 (MQ=255)
                            gTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATg                                                                       <  2:3489946/100‑1 (MQ=255)
                            gTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATg                                                                       >  1:2240945/1‑100 (MQ=255)
                            gTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATg                                                                       <  2:3132169/100‑1 (MQ=255)
                             tGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGa                                                                      >  2:2171782/1‑100 (MQ=255)
                             tGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGa                                                                      >  2:2412907/1‑100 (MQ=255)
                              gaagagTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                     >  1:859635/1‑100 (MQ=255)
                              gaagaaTTTACGGTCATAGATCACGTCAACATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATg                                                                       <  2:3807637/98‑1 (MQ=255)
                              gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                     <  1:936423/100‑1 (MQ=255)
                              gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                     >  2:1429165/1‑100 (MQ=255)
                              gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                     >  1:3122989/1‑100 (MQ=255)
                              gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                     >  2:1086351/1‑100 (MQ=255)
                               aagtaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGAAAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACa                                                                    <  1:127152/100‑1 (MQ=255)
                               aagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGCGGGAAGACAACCTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGa                                                                      <  1:173800/98‑1 (MQ=255)
                               aagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACa                                                                    >  1:3317788/1‑100 (MQ=255)
                                agaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATGAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACAt                                                                   <  2:2447188/100‑1 (MQ=255)
                                 gaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATg                                                                  <  1:98221/100‑1 (MQ=255)
                                     ttACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTc                                                              <  1:1305101/100‑1 (MQ=255)
                                      tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcg                                                             >  1:540237/1‑100 (MQ=255)
                                      tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcg                                                             >  2:1914956/1‑100 (MQ=255)
                                       aCGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgc                                                            <  1:264340/100‑1 (MQ=255)
                                       aCGGTCATAGATCACGTCAAAATAAGAAGAGGGAAAACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGCAATGACATGCGTcgc                                                            <  2:3043908/100‑1 (MQ=255)
                                           tCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGATATGTACATTAGTACACAATTGCAGAGGGAGGAAACGACATGCGTCgcgc                                                          >  2:1635727/1‑98 (MQ=255)
                                            cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                       >  1:2459153/1‑100 (MQ=255)
                                             aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                      <  2:1993201/100‑1 (MQ=255)
                                             aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                      <  2:2556179/100‑1 (MQ=255)
                                              tAGACCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGACTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                     >  2:603374/1‑100 (MQ=255)
                                               aGATCACGTCAAAATAAGAAGTGGGAAGACAACATTAGCAAGGTTGTGTACATTGGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGa                                                    >  2:3238574/1‑100 (MQ=255)
                                                 aTCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATc                                                  >  1:799643/1‑100 (MQ=255)
                                                     cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                              >  2:75060/1‑100 (MQ=255)
                                                       tCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAAc                                            >  1:2533744/1‑100 (MQ=255)
                                                         aaaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACgg                                          >  1:629709/1‑100 (MQ=255)
                                                          aaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGAGCTAACGATCCGCAACGGc                                         >  1:317411/1‑100 (MQ=255)
                                                           naTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATAGCAGTGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAAAGGCg                                        >  2:1335529/2‑100 (MQ=255)
                                                            aTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCga                                       >  1:3702924/1‑100 (MQ=255)
                                                             tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgtg                                      >  2:2309148/1‑98 (MQ=255)
                                                             tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                      <  2:66128/100‑1 (MQ=255)
                                                             tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                        >  2:3601056/1‑98 (MQ=255)
                                                                gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaa                                   >  2:1093033/1‑100 (MQ=255)
                                                                gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaa                                   >  1:586317/1‑100 (MQ=255)
                                                                   gagGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaaa                                <  2:3448704/100‑1 (MQ=255)
                                                                     gggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAAtt                              >  2:920946/1‑100 (MQ=255)
                                                                      ggAAGACAACATTAGCAATGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTa                             <  2:1390290/100‑1 (MQ=255)
                                                                      ggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTa                             <  1:1086360/100‑1 (MQ=255)
                                                                       gAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTCACGATCCGCAACGGCGAGAAAAAATTa                             >  2:13728/1‑99 (MQ=255)
                                                                          gACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTACCGATCCGCAACGGCGAGAAAAAATTATCCa                         >  1:769199/1‑100 (MQ=255)
                                                                             aacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGc                      <  1:3072242/100‑1 (MQ=255)
                                                                              acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                     >  1:3780117/1‑100 (MQ=255)
                                                                               caTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCa                    <  1:1346262/100‑1 (MQ=255)
                                                                                aTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCAc                   >  1:2795964/1‑100 (MQ=255)
                                                                                   aGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGATGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                <  1:771790/100‑1 (MQ=255)
                                                                                   aGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                >  2:3667766/1‑100 (MQ=255)
                                                                                    gCAAGGTTGTGTACATTTGTACCCAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                >  2:2970154/1‑99 (MQ=255)
                                                                                    gCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg               >  1:283037/1‑100 (MQ=255)
                                                                                    ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGAGAGAAAAAATTATCCAGGCCACGCTg               >  2:3471684/7‑100 (MQ=38)
                                                                                     tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg              >  1:1902763/6‑100 (MQ=255)
                                                                                     cAAGGTTGTGTACATTTGTACACACTTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAACCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg              >  2:1232769/1‑100 (MQ=255)
                                                                                       caattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg            >  2:1097898/4‑100 (MQ=255)
                                                                                       caattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg            >  1:3415441/4‑100 (MQ=255)
                                                                                       caattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg            >  1:397961/4‑100 (MQ=255)
                                                                                         gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc          >  2:1648944/1‑100 (MQ=255)
                                                                                         gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc          <  1:2309163/100‑1 (MQ=255)
                                                                                         attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc          >  1:2094228/2‑100 (MQ=255)
                                                                                         attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc          >  2:1960497/2‑100 (MQ=255)
                                                                                           tgtgTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg        >  1:3393772/1‑100 (MQ=255)
                                                                                            gtgtACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt       >  2:3728302/1‑100 (MQ=255)
                                                                                             tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTg      >  1:1532414/1‑100 (MQ=255)
                                                                                               tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa    >  1:1395669/1‑100 (MQ=255)
                                                                                               tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa    <  1:1189855/100‑1 (MQ=255)
                                                                                                aCATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTAGCCAGGCCACGCTGGAGGCGGTGaaa   <  1:1111055/100‑1 (MQ=255)
                                                                                                 cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc  >  2:3746956/1‑100 (MQ=255)
                                                                                                 cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc  <  2:200402/100‑1 (MQ=255)
                                                                                                 cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc  <  2:682996/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAC  >  NC_000913/887293‑887489

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