Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 887,393 T→A 18.3% 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,3930TA18.3% 245.2 / 44.9 107intergenic (‑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 (44/43);  new base A (11/9);  total (55/52)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

Only 100 of 109 total aligned reads displayed.

CCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACT  >  NC_000913/887295‑887490
                                                                                                  |                                                                                                 
ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAttt                                                                                                  >  1:2137279/1‑100 (MQ=255)
  aGGTTGCGCGTTTCAATGCATGACGGGAAGAATTTACGGTCATAGATCACGGCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                   <  1:2381077/97‑1 (MQ=255)
   ggTTCGGCTTTTCAATGCATTACGTGAAGAATTTGCATTCATAGATCGCGTCAAAATAAGAAGAGGGAAGTCAACATTAGCAAGGTTGTGTACATTTGTa                                                                                               <  1:2337280/100‑1 (MQ=255)
   ggTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                               <  1:2172313/100‑1 (MQ=255)
   ggTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                               <  1:549666/100‑1 (MQ=255)
   ggTTCGGCGTTTCAATGAATTACGTGAAGAAGTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                               <  1:439435/100‑1 (MQ=255)
      tCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTacac                                                                                            <  2:1221496/100‑1 (MQ=255)
       cGGCGTTTCAATGCTTTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAcaca                                                                                           <  1:1764957/100‑1 (MQ=255)
        ggCGTTACAATGCATTACGGAAACAATTTACGGTCATAGATCACGTCAAAATAAGGAGAGGGAAGACAACATTAGCAAGGTCGTGTACATAGGTACACaa                                                                                          >  1:1585272/1‑100 (MQ=255)
         gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAGGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  2:632496/100‑2 (MQ=255)
         gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAAt                                                                                         <  1:1830583/100‑1 (MQ=255)
         gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  2:492328/100‑2 (MQ=255)
         gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  2:2110815/100‑2 (MQ=255)
         gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  1:1046987/100‑2 (MQ=255)
         gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  2:553399/100‑2 (MQ=255)
          cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAAtt                                                                                        >  2:2018758/1‑100 (MQ=255)
          cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                        <  2:1679113/100‑3 (MQ=255)
           gTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTCGCAAGGTTGTGTACATTTGTACACAAtt                                                                                        >  1:1594697/1‑99 (MQ=255)
            tttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGc                                                                                      >  2:2517729/1‑100 (MQ=255)
             ttCAATGCATTACGTGAAGACTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCa                                                                                     >  1:1585907/1‑100 (MQ=255)
             ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTTCACAATTg                                                                                       >  1:1327893/1‑98 (MQ=255)
             ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCa                                                                                     <  1:1733908/100‑1 (MQ=255)
             ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaaccttg                                                                                     <  2:59947/100‑6 (MQ=255)
              tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    <  2:865616/100‑1 (MQ=255)
              tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    >  2:858511/1‑100 (MQ=255)
              tCAAGGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    <  2:648495/100‑1 (MQ=255)
               cAATGCATTACGTGAAAAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacct                                                                                       <  1:2421286/96‑4 (MQ=255)
                aaTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagag                                                                                  <  2:1694655/100‑1 (MQ=255)
                aaTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagag                                                                                  <  2:456004/100‑1 (MQ=255)
                  tGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                                >  2:44078/1‑100 (MQ=255)
                   gCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                               <  2:1879434/100‑1 (MQ=255)
                   gCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                               >  1:145064/1‑100 (MQ=255)
                    cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                              <  2:1327771/100‑1 (MQ=255)
                     aTTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                             <  1:343843/100‑1 (MQ=255)
                        aCGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAATATGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                             <  1:1790181/97‑1 (MQ=255)
                        aCGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGaaa                                                                          >  1:1396246/1‑100 (MQ=255)
                        aCGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGaaa                                                                          >  2:526142/1‑100 (MQ=255)
                         cGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTTCACAATTGCAGAGAGAGGAAAt                                                                         <  1:1150573/100‑1 (MQ=255)
                           tGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGa                                                                       >  2:2481393/1‑100 (MQ=255)
                           ngaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGa                                                                       >  1:1289833/2‑100 (MQ=255)
                            gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                      >  2:484877/1‑100 (MQ=255)
                             aagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACa                                                                     >  1:1616941/1‑100 (MQ=255)
                              agaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACAt                                                                    <  1:2488745/100‑1 (MQ=255)
                                 aTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCg                                                                 >  1:940817/1‑100 (MQ=255)
                                  tttACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCgn                                                                <  2:72509/100‑2 (MQ=255)
                                      cGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAAGATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgc                                                             >  2:26105/1‑99 (MQ=255)
                                       ggTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCgcgc                                                           <  1:587961/100‑1 (MQ=255)
                                          cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                        >  1:824315/1‑100 (MQ=255)
                                          cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                        >  2:211274/1‑100 (MQ=255)
                                           aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                       <  2:1177141/100‑1 (MQ=255)
                                           aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                       <  1:1089108/100‑1 (MQ=255)
                                           aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                       >  2:2440011/1‑100 (MQ=255)
                                            tAGATCACGTCAACATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                      >  1:668959/1‑100 (MQ=255)
                                            tAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                      >  2:1462894/1‑100 (MQ=255)
                                             aGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGa                                                     >  2:257497/1‑100 (MQ=255)
                                              gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                    >  2:852756/1‑100 (MQ=255)
                                              gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                    >  2:962349/1‑100 (MQ=255)
                                                 cACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                 <  1:429676/100‑1 (MQ=255)
                                                   cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                               >  1:849302/1‑100 (MQ=255)
                                                    gTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                               >  2:729870/1‑99 (MQ=255)
                                                       aaaaTAAGAAGGGGGAATACAACATTAGCCAGGTTGTGGACAGTTGTACACAAGTGCAGAGAGAGGAAAGGACATGCGTCGCGCTAACGATCCGCAACg                                            <  1:2517710/99‑1 (MQ=255)
                                                        aaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                          <  1:93821/100‑1 (MQ=255)
                                                          aTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCga                                        >  2:1084937/1‑100 (MQ=255)
                                                           tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                       <  2:1562378/100‑1 (MQ=255)
                                                           tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                       >  1:955254/1‑100 (MQ=255)
                                                             agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                          >  1:1948536/1‑95 (MQ=255)
                                                             agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaa                                     >  2:61721/1‑100 (MQ=255)
                                                              gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaa                                    >  2:1610090/1‑100 (MQ=255)
                                                                  agGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAAt                                <  2:2535274/100‑1 (MQ=255)
                                                                   gggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGCGAGAGGAAATGACATGCGTCGCGATAACGATCCGCAACGGCGAGAAAAAAt                                >  1:981066/1‑99 (MQ=255)
                                                                   gggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAAtt                               >  1:690985/1‑100 (MQ=255)
                                                                    ggAAGAAAAGATGAGAAAGGTGGTGTAGATTTGTACACAATTGCAGTGAGAGGAAAGGACATGAATCGCGCTAACGATTTGCAACAGAGAAAAAAAATTa                              <  1:2137233/100‑1 (MQ=255)
                                                                        gACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCa                          <  1:993542/100‑1 (MQ=255)
                                                                         acaacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                         >  2:1926914/1‑100 (MQ=255)
                                                                         acaacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                         <  2:51858/100‑1 (MQ=255)
                                                                          caacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAgg                        >  1:1089824/1‑100 (MQ=255)
                                                                          caacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAgg                        >  2:1717573/1‑100 (MQ=255)
                                                                            acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                      <  1:2039282/100‑1 (MQ=255)
                                                                                ctctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGc                  >  2:1329831/9‑100 (MQ=38)
                                                                                 tctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                 >  1:1400243/8‑100 (MQ=38)
                                                                                  gCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                <  1:343622/100‑1 (MQ=255)
                                                                                  gCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                >  2:2126368/1‑100 (MQ=255)
                                                                                  ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                >  2:894065/7‑100 (MQ=39)
                                                                                   tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCTAGAAAAAATTATCCAGGCCAc                    >  1:2565983/6‑95 (MQ=38)
                                                                                   tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg               >  2:1456080/6‑100 (MQ=255)
                                                                                   tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg               >  1:1802354/6‑100 (MQ=255)
                                                                                    gcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGa              >  2:469423/5‑100 (MQ=255)
                                                                                     aGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg             <  2:1845587/100‑1 (MQ=255)
                                                                                       attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc           >  1:2560806/2‑100 (MQ=255)
                                                                                         tgtgTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg         >  1:1388331/1‑100 (MQ=255)
                                                                                          gtgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt        <  2:1711957/100‑1 (MQ=255)
                                                                                          gtgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt        <  2:655368/100‑1 (MQ=255)
                                                                                          gtgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt        <  1:350483/100‑1 (MQ=255)
                                                                                           tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTg       <  1:2134306/100‑1 (MQ=255)
                                                                                            gtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa      <  1:37259/100‑1 (MQ=255)
                                                                                            gtACATTTGTACACAATTGCAGAGAGAGGAAATGACACGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa      >  1:1378601/1‑100 (MQ=255)
                                                                                             tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa     <  2:1856053/100‑1 (MQ=255)
                                                                                              aCATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaaa    >  1:1537930/1‑100 (MQ=255)
                                                                                              aCATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaaa    >  2:807330/1‑100 (MQ=255)
                                                                                                aTTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACt  >  2:670643/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
CCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACT  >  NC_000913/887295‑887490

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: