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% 401.0 / 74.0 168intergenic (‑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 (80/57);  new base A (14/17);  total (94/74)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

CCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACT  >  NC_000913/887294‑887490
                                                                                                   |                                                                                                 
cccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                   <  1:2880250/100‑1 (MQ=255)
cccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                   <  1:1328959/100‑1 (MQ=255)
  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 <  2:2858562/100‑1 (MQ=255)
  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 >  2:4131010/1‑100 (MQ=255)
  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 >  1:346832/1‑100 (MQ=255)
  cAGGTTCGGCGTGTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 <  2:2414717/100‑1 (MQ=255)
   aGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTCGCAAGGTTGTGTACATTTGt                                                                                                <  1:2699245/100‑1 (MQ=255)
    ggTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                               >  1:3701178/1‑100 (MQ=255)
    ggTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                               >  2:2666193/1‑100 (MQ=255)
    ggTTCGGCGTTTCAAGGCATTAGGTGAAGAATTTACGGGCATAGATCAAGTCAAAATAAGAAGAGAAAAGACAGAAGTAGCAAGGGTGGGTAGAtgtg                                                                                                 >  1:3673830/1‑98 (MQ=255)
     gTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTac                                                                                              >  1:2017557/1‑100 (MQ=255)
     gTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTac                                                                                              <  1:1532441/100‑1 (MQ=255)
      ntCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTaca                                                                                             >  2:3477154/2‑100 (MQ=255)
         ggCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACaa                                                                                          >  2:836606/1‑100 (MQ=255)
         ggCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACaa                                                                                          >  1:295993/1‑100 (MQ=255)
         ggCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaa                                                                                          <  2:1904765/100‑1 (MQ=255)
          gggTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  2:2280531/98‑2 (MQ=255)
          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGTTTGTGTACAAATGTACACaac                                                                                         <  1:959889/100‑2 (MQ=255)
          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  1:2953059/100‑2 (MQ=255)
          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  1:2850050/100‑2 (MQ=255)
            gTTTCAATGCATTACGTGAAGAATTTTCGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacct                                                                                       <  1:3908426/100‑4 (MQ=255)
            gTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTg                                                                                       <  1:490411/100‑1 (MQ=255)
            gTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTg                                                                                       <  2:2945433/100‑1 (MQ=255)
             tttCAATGCATTACGTGAAGAATTTACGGTCATAGATCATGGCAAACTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacct                                                                                       <  1:3605552/99‑4 (MQ=255)
             tttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGc                                                                                      <  1:4367428/100‑1 (MQ=255)
             tttCAATGCATTACGTGAAGAATTTACGGTCAAAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGc                                                                                      >  2:4191540/1‑100 (MQ=255)
              ttCAATGCATTACGTGAAGAAGTTACGGTCATATATCACGTCAAAATAAGAAGATGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                        <  1:613755/97‑3 (MQ=255)
               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCan                                                                                    <  1:283234/100‑2 (MQ=255)
               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    <  2:124909/100‑1 (MQ=255)
               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    <  2:4180097/100‑1 (MQ=255)
               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    <  1:620164/100‑1 (MQ=255)
                cAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCaga                                                                                   >  1:3798884/1‑100 (MQ=255)
                cAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    >  2:850201/1‑99 (MQ=255)
                cAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaaccttgct                                                                                   <  1:3584854/100‑8 (MQ=38)
                cAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaaccttgct                                                                                   <  2:3802457/100‑8 (MQ=38)
                 aaTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagag                                                                                  >  1:1705141/1‑100 (MQ=255)
                     cATTACGTGCAGAATTTACGGTGATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATATGTACCCAATTCCagagaga                                                                               >  1:248040/1‑99 (MQ=255)
                     cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                              <  2:399755/100‑1 (MQ=255)
                       ttACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGa                                                                            <  2:2086586/100‑1 (MQ=255)
                           gTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATg                                                                        >  2:1516766/1‑100 (MQ=255)
                            tGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGa                                                                       >  2:3632417/1‑100 (MQ=255)
                            tGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGa                                                                       >  1:4204261/1‑100 (MQ=255)
                             gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                      <  2:600197/100‑1 (MQ=255)
                             gaagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGAc                                                                      >  2:1360545/1‑100 (MQ=255)
                                gaaCTTACGGTCATAGATCACGTCAAAAGAATAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATg                                                                   <  2:4398419/100‑1 (MQ=255)
                                 aaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGc                                                                  <  2:1057820/100‑1 (MQ=255)
                                  aTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGGCAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCg                                                                 <  2:3516637/100‑1 (MQ=255)
                                    ttACGGACATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACCATTACAGAGAGAGGAAATGAAATGCGTc                                                               >  1:4393366/1‑100 (MQ=255)
                                       cGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgcg                                                            >  2:3527923/1‑100 (MQ=255)
                                       cGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgcg                                                            >  2:1293821/1‑100 (MQ=255)
                                         gTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCt                                                          <  2:1380453/100‑1 (MQ=255)
                                           cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGGGGAAATGACATGCGTCGCGCTaa                                                        >  1:3406248/1‑100 (MQ=255)
                                           cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                        >  1:1120830/1‑100 (MQ=255)
                                            aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                       <  2:3548679/100‑1 (MQ=255)
                                             tAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                      >  2:909503/1‑100 (MQ=255)
                                             tAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGAAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGCAATCACATGCGTCGCGCTAAc                                                       >  1:1316810/1‑99 (MQ=255)
                                               gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                    >  1:836914/1‑100 (MQ=255)
                                                aTCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                      >  1:4156903/1‑97 (MQ=255)
                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                               >  2:1264151/1‑100 (MQ=255)
                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                               >  1:2398474/1‑100 (MQ=255)
                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCa                                               <  1:187634/100‑1 (MQ=255)
                                                     gTCCAAATATGAAAAGGGTAGACAACATTAGCAAGGATTTGTACATTTGAACACAACTTCAGATTGAGGAAATGACAAACGTCGCGCTAActaacc                                                  >  1:125071/1‑91 (MQ=255)
                                                      tCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAAc                                             <  1:2229660/100‑1 (MQ=255)
                                                       cAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACg                                            >  1:563072/1‑100 (MQ=255)
                                                           aTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCga                                        <  2:2652467/100‑1 (MQ=255)
                                                            tAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgag                                       <  2:2186434/100‑1 (MQ=255)
                                                              agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaa                                     >  2:1478421/1‑100 (MQ=255)
                                                              agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaa                                     >  1:857420/1‑100 (MQ=255)
                                                               gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGAAACGGCGAGaaa                                    >  1:350517/1‑100 (MQ=255)
                                                                      taaGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCATAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGAGAGAAAAAATTAt                             >  1:1911775/2‑100 (MQ=255)
                                                                          acaacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                         <  1:2114995/100‑1 (MQ=255)
                                                                             acaTTAGCAAGGTTGTGTACATTTGTACTCAATTGCAGAGAGAGGAAATGACATGCGTAGCGCTAACGATCCGACACGGCGAGAAACAATTATCCAGGcc                      >  1:2667315/1‑100 (MQ=255)
                                                                             acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                      >  1:1054641/1‑100 (MQ=255)
                                                                             acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                      <  1:2245620/100‑1 (MQ=255)
                                                                             acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                      <  2:3192854/100‑1 (MQ=255)
                                                                             acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                      >  2:3035246/1‑100 (MQ=255)
                                                                              caTTAGCTAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCa                     >  2:2311741/1‑100 (MQ=255)
                                                                                ttAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACg                   >  2:1753631/1‑100 (MQ=255)
                                                                                  aGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                 >  1:3065298/1‑100 (MQ=255)
                                                                                   gCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                >  1:4134440/1‑100 (MQ=255)
                                                                                   ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCAGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                >  1:3433165/7‑100 (MQ=38)
                                                                                    tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGTGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg               >  2:2763975/6‑100 (MQ=255)
                                                                                    cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg               <  2:3905285/100‑1 (MQ=255)
                                                                                     gcaattGTGTACAAATGTACACAATTGCAGAGAGAGTAAATGACAGGCGTCGCGCTAACGTTCCGCAAAGGCGAGAAAAAATTATCCAGGCTACGCTGGa              >  1:2278852/5‑100 (MQ=255)
                                                                                      caattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg             >  2:859280/4‑100 (MQ=255)
                                                                                      aGGTTGTGTACATTTGTACACAATAGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATGCGCAACGGCGTGAAAAAATTATCCAGGCCAGGCTgg               <  1:908955/98‑1 (MQ=255)
                                                                                       ggTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAgg            >  2:3872646/1‑100 (MQ=255)
                                                                                        gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc           <  2:3352217/100‑1 (MQ=255)
                                                                                        attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc           >  1:2319999/2‑100 (MQ=255)
                                                                                          tgtgTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg         >  1:3210835/1‑100 (MQ=255)
                                                                                            tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTg       >  2:2562420/1‑100 (MQ=255)
                                                                                            tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTg       <  1:3594125/100‑1 (MQ=255)
                                                                                             gtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa      >  2:3907140/1‑100 (MQ=255)
                                                                                             gtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa      <  1:2773521/100‑1 (MQ=255)
                                                                                               ncATTTGTACACAATTGCAGAGAGAGGAAATGACATGAGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaaa    >  1:2365017/2‑100 (MQ=255)
                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc   <  2:3702719/100‑1 (MQ=255)
                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc   <  1:2751731/100‑1 (MQ=255)
                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc   >  1:4182057/1‑100 (MQ=255)
                                                                                                cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc   <  2:3323917/100‑1 (MQ=255)
                                                                                                  tttGTACACAATTGCAGAGAGAGCAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGCAAAAATTATCCAGACCACGCTGGAGGCGGTGAAACt  >  2:3735870/1‑99 (MQ=255)
                                                                                                   |                                                                                                 
CCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACT  >  NC_000913/887294‑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: