Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 887,393 T→A 17.8% intergenic (‑54/‑30) 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,3930TA17.8% 334.4 / 65.3 146intergenic (‑54/‑30)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 = 8.92e-02

Only 100 of 150 total aligned reads displayed.

CCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTT  >  NC_000913/887294‑887492
                                                                                                   |                                                                                                   
cccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGGAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                     >  1:2698777/1‑100 (MQ=255)
 ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAttt                                                                                                    >  2:183209/1‑100 (MQ=255)
  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                   <  1:2906716/100‑1 (MQ=255)
  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                   <  1:202593/100‑1 (MQ=255)
  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                   <  1:1660925/100‑1 (MQ=255)
    ggTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                 <  1:3276601/100‑1 (MQ=255)
       tCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTacac                                                                                              <  1:3491141/100‑1 (MQ=255)
        cGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAcaca                                                                                             >  1:3448730/1‑100 (MQ=255)
         ggCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaa                                                                                            <  1:3689052/100‑1 (MQ=255)
          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                           <  2:679600/100‑2 (MQ=255)
          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                           <  2:1133738/100‑2 (MQ=255)
          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                           <  1:2826365/100‑2 (MQ=255)
          gCGTTTCAATGCATTACGTAAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGTCAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                           <  1:989048/100‑2 (MQ=255)
           cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                          <  1:2160624/100‑3 (MQ=255)
            gTTTCAATGCATTACGTGAAAAATTTACGGTCATAGATCACGTCAAAAGAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                           <  2:570502/98‑2 (MQ=255)
             tttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTg                                                                                         >  1:105942/1‑99 (MQ=255)
             tttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGc                                                                                        <  1:3634865/100‑1 (MQ=255)
              ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCa                                                                                       <  2:733133/100‑1 (MQ=255)
               taaaTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaaccttgc                                                                                      <  2:2452414/98‑7 (MQ=39)
               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                      >  2:2896018/1‑100 (MQ=255)
                cAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCaga                                                                                     <  1:2814725/100‑1 (MQ=255)
                  aTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagaga                                                                                   >  2:3668233/1‑100 (MQ=255)
                   tGCATTACGTTAAGAGTTTTCGGTCAGATATCACGTCAAAATAAGAAGAGGTAAGACAACATTAGCAAGGGTGTGTACAAATGTACACaacctttc                                                                                      <  2:1315043/96‑7 (MQ=255)
                   attcgcccatGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                                  <  1:837438/91‑1 (MQ=255)
                     cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGCAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                                >  2:811642/1‑100 (MQ=255)
                     cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                                <  2:1767789/100‑1 (MQ=255)
                      aTTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                               >  1:2469710/1‑100 (MQ=255)
                          cGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAAt                                                                           >  1:815300/1‑100 (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                                                                        >  2:1429165/1‑100 (MQ=255)
                             gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                        >  2:1092019/1‑100 (MQ=255)
                             gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                        >  2:1086351/1‑100 (MQ=255)
                             gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                        <  1:936423/100‑1 (MQ=255)
                              aagtaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGAAAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACa                                                                       <  1:127152/100‑1 (MQ=255)
                              aagaaTTTACGGTCATAGATCACGTAAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                        <  2:2613688/99‑1 (MQ=255)
                               agaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACAt                                                                      >  2:3138329/1‑100 (MQ=255)
                                 aaTTTACGGTCATAGATCACGGCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATg                                                                     <  2:3029030/99‑1 (MQ=255)
                                    ttACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTc                                                                 <  1:1305101/100‑1 (MQ=255)
                                     tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcg                                                                >  2:1914956/1‑100 (MQ=255)
                                     tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcg                                                                >  1:540237/1‑100 (MQ=255)
                                      aCGGTCATAGATCACGTCAAAATAAGAAGAGGGAAAACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGCAATGACATGCGTcgc                                                               <  2:3043908/100‑1 (MQ=255)
                                         gTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCt                                                            >  2:2265081/1‑100 (MQ=255)
                                          tCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGATATGTACATTAGTACACAATTGCAGAGGGAGGAAACGACATGCGTCgcgc                                                             >  2:1635727/1‑98 (MQ=255)
                                            aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                         <  2:2556179/100‑1 (MQ=255)
                                            aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                         <  2:1993201/100‑1 (MQ=255)
                                             ttgATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                        >  1:2593725/3‑100 (MQ=255)
                                             tAGACCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGACTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                        >  2:603374/1‑100 (MQ=255)
                                              aGATCACGTCAAAATAAGAAGTGGGAAGACAACATTAGCAAGGTTGTGTACATTGGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGa                                                       >  2:3238574/1‑100 (MQ=255)
                                               gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                      >  1:3776647/1‑100 (MQ=255)
                                               gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGAAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                      >  2:449099/1‑100 (MQ=255)
                                                aTCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATc                                                     >  1:799643/1‑100 (MQ=255)
                                                   aCGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGc                                                  <  2:3776619/100‑1 (MQ=255)
                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                                 >  1:95930/1‑100 (MQ=255)
                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                                 >  2:2435675/1‑100 (MQ=255)
                                                     gTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCaa                                                >  1:2753950/1‑100 (MQ=255)
                                                      tCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAAc                                               >  1:2533744/1‑100 (MQ=255)
                                                        aaaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACgg                                             >  1:629709/1‑100 (MQ=255)
                                                         aaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                            >  1:398392/1‑100 (MQ=255)
                                                         aaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGAGCTAACGATCCGCAACGGc                                            >  1:317411/1‑100 (MQ=255)
                                                            taaggaaggggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                         <  2:730884/92‑1 (MQ=255)
                                                            tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                         >  1:2472746/1‑100 (MQ=255)
                                                            tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                           >  2:3601056/1‑98 (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)
                                                                           caacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATTCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAGTTATCCAgg                          >  2:602887/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)
                                                                                 tAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAACATTATCCAGGCCACGc                    >  2:1049948/1‑100 (MQ=255)
                                                                                 ctctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGc                    >  1:2166338/9‑100 (MQ=38)
                                                                                  aGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGATGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                   <  1:771790/100‑1 (MQ=255)
                                                                                   ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                  >  2:1326045/7‑100 (MQ=39)
                                                                                   ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGAGAGAAAAAATTATCCAGGCCACGCTg                  >  2:3471684/7‑100 (MQ=38)
                                                                                    tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg                 >  1:1902763/6‑100 (MQ=255)
                                                                                      caattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg               >  2:1097898/4‑100 (MQ=255)
                                                                                        gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc             >  2:1648944/1‑100 (MQ=255)
                                                                                        gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc             <  1:2309163/100‑1 (MQ=255)
                                                                                        attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGCAAAAATTATCCAGGCCACGCTGGAGGc             >  2:674292/2‑100 (MQ=255)
                                                                                        attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc             >  1:3371045/2‑100 (MQ=255)
                                                                                        attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc             >  2:1960497/2‑100 (MQ=255)
                                                                                          tgtgTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg           >  1:3393772/1‑100 (MQ=255)
                                                                                           gtgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt          <  1:3260937/100‑1 (MQ=255)
                                                                                           gtgtACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt          >  2:3728302/1‑100 (MQ=255)
                                                                                            tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTg         >  1:1532414/1‑100 (MQ=255)
                                                                                              tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa       >  2:3494444/1‑100 (MQ=255)
                                                                                              tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa       <  1:1189855/100‑1 (MQ=255)
                                                                                              tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa       <  1:2390831/100‑1 (MQ=255)
                                                                                              tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa       <  1:3168219/100‑1 (MQ=255)
                                                                                               aCATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTAGCCAGGCCACGCTGGAGGCGGTGaaa      <  1:1111055/100‑1 (MQ=255)
                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc     <  1:1359268/100‑1 (MQ=255)
                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc     <  2:200402/100‑1 (MQ=255)
                                                                                                   ttGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACttt  <  2:1344357/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTT  >  NC_000913/887294‑887492

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: