Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 887,392 T→A 17.6% 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.6% 391.8 / 75.1 165intergenic (‑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 (70/66);  new base A (15/14);  total (85/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

Only 100 of 166 total aligned reads displayed.

CAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACT  >  NC_000913/887296‑887490
                                                                                                |                                                                                                  
cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 <  1:4374037/100‑1 (MQ=255)
cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 <  1:994404/100‑1 (MQ=255)
cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 >  1:3782302/1‑100 (MQ=255)
cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                 <  1:2056416/100‑1 (MQ=255)
   gTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTac                                                                                              <  2:4107888/100‑1 (MQ=255)
       ggCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACaa                                                                                          >  2:2607626/1‑100 (MQ=255)
       ggCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACaa                                                                                          >  1:1712401/1‑100 (MQ=255)
        ggttttCAATGCATTACGTGAAGAAGTTATGGTCATAGATCACGTCAAAATAAGAAGATGGAAGACAACATTAGCAAGGTTGTGTACAAACGTACACaac                                                                                         <  1:2788194/97‑2 (MQ=255)
        gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAAt                                                                                         >  2:1636050/1‑100 (MQ=255)
        gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  1:3384938/100‑2 (MQ=255)
        gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaac                                                                                         <  2:2392877/100‑2 (MQ=255)
         cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                        <  2:3164584/100‑3 (MQ=255)
         cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                        <  2:1214885/100‑3 (MQ=255)
         cGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaacc                                                                                        <  1:712102/100‑3 (MQ=255)
          gggtCAAGGCATGACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTg                                                                                       <  1:809738/97‑1 (MQ=255)
           tttCAATTCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGc                                                                                      >  1:2451818/1‑100 (MQ=255)
           tttCAATGCATTACGTGAAGAATTTACGGTCATAGCTCACGTCAAAATAAGAAGAGGGAAGACAACATTAGTAAGGTTGTGTACATTTGTACACAATTg                                                                                       <  1:868927/99‑1 (MQ=255)
            ttCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaaccttg                                                                                     <  2:687411/100‑6 (MQ=255)
             tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGATGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCag                                                                                    <  1:1008009/100‑1 (MQ=255)
                  gCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                               <  2:2519309/100‑1 (MQ=255)
                  gCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                               >  1:4174751/1‑100 (MQ=255)
                   cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                              >  1:3937469/1‑100 (MQ=255)
                   cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                              >  1:3626166/1‑100 (MQ=255)
                   cATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagagag                                                                              <  1:3233664/100‑1 (MQ=255)
                     ttACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGa                                                                            <  1:257722/100‑1 (MQ=255)
                     ttACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGa                                                                            >  2:1119694/1‑100 (MQ=255)
                     ttACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGa                                                                            <  1:1132579/100‑1 (MQ=255)
                      tatgTGAAGAATTTACGGTCATAGATCACGTAAAAATACTAAGATTTAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGa                                                                            <  1:1808490/96‑1 (MQ=255)
                      tACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGAAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGaa                                                                           >  1:2133747/1‑100 (MQ=255)
                       aCGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGaaa                                                                          >  1:2964948/1‑100 (MQ=255)
                         gTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATg                                                                        <  1:2128965/100‑1 (MQ=255)
                            aagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACa                                                                     <  2:1933523/100‑1 (MQ=255)
                            aagaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACa                                                                     <  1:572019/100‑1 (MQ=255)
                             agaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACAt                                                                    >  1:1426464/1‑100 (MQ=255)
                             agaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACAt                                                                    >  1:2873414/1‑100 (MQ=255)
                              gaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCATGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATg                                                                   <  2:2535083/100‑1 (MQ=255)
                              gaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATg                                                                   >  2:190526/1‑100 (MQ=255)
                              gaaTTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATg                                                                   >  1:4260620/1‑100 (MQ=255)
                                  ttACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTc                                                               <  1:272470/100‑1 (MQ=255)
                                      ggTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCgcgc                                                           >  2:2468727/1‑100 (MQ=255)
                                      ggTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCgcgc                                                           <  1:3461157/100‑1 (MQ=255)
                                       gTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCt                                                          >  2:213567/1‑100 (MQ=255)
                                         cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                        >  1:4029191/1‑100 (MQ=255)
                                         cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAAACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTa                                                         >  1:2685655/1‑99 (MQ=255)
                                          aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                       <  2:586411/100‑1 (MQ=255)
                                          aTAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAAc                                                       <  2:924388/100‑1 (MQ=255)
                                           tAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACg                                                      >  2:2990083/1‑100 (MQ=255)
                                             gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGa                                                     >  1:1380260/1‑99 (MQ=255)
                                                cACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                 >  2:3120368/1‑100 (MQ=255)
                                                   gTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCaa                                              <  1:117075/100‑1 (MQ=255)
                                                     cAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACg                                            >  1:4049142/1‑100 (MQ=255)
                                                     cAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACg                                            >  1:3393314/1‑100 (MQ=255)
                                                     cAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACg                                            >  2:1410172/1‑100 (MQ=255)
                                                      aaaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACgg                                           <  1:162044/100‑1 (MQ=255)
                                                       aaaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                          >  1:1006776/1‑100 (MQ=255)
                                                        aaTAAGAAGAGGGAAGACAACATTAGCAGGGTTGTGTACATTTGTACACAATTTCAGAGAGAGGAAATGACATGCTTCGCGCTAACGATCCGCAACGGCg                                         <  1:2438317/100‑1 (MQ=255)
                                                        aaTAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                         <  2:3636975/100‑1 (MQ=255)
                                                            agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgaga                                      >  2:2387404/1‑99 (MQ=255)
                                                            agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaa                                     >  1:4329039/1‑100 (MQ=255)
                                                            agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAACTGCAGAGAGAGGAAATGACATGCGTCACGCTAACGACCCTCAACGGCgaga                                      <  1:2680932/99‑1 (MQ=255)
                                                               agagGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaa                                  <  2:3782326/100‑1 (MQ=255)
                                                                gagGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaaa                                 <  2:3247524/100‑1 (MQ=255)
                                                                gagGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaaa                                 >  2:2054979/1‑100 (MQ=255)
                                                                  gggAAGACAACATTAGCAAGGTTGTGTACATTTGTACACGATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAAtt                               >  1:4012074/1‑100 (MQ=255)
                                                                    gAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTAt                             >  2:2861574/1‑100 (MQ=255)
                                                                      aGAAAACATTAGCAAGGTTGTGTACATTTGTAAACAATTGCAGAGAGAGGAAATGACATGCGTGTAGATAACAATCCGCAGGGGCGAGAAAAAATTATc                            >  1:3938406/1‑99 (MQ=255)
                                                                       gACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCa                          >  2:1945888/1‑100 (MQ=255)
                                                                        acaacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                         >  2:1815453/1‑100 (MQ=255)
                                                                           acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAgg                        <  1:1066124/98‑1 (MQ=255)
                                                                           acaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGcc                      >  2:4403343/1‑100 (MQ=255)
                                                                            caTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCa                     >  2:2956279/1‑100 (MQ=255)
                                                                               tAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGc                  >  1:3307162/1‑100 (MQ=255)
                                                                                tctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                 >  1:1988260/8‑100 (MQ=38)
                                                                                aGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                 <  1:4343885/100‑1 (MQ=255)
                                                                                 ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                >  2:440118/7‑100 (MQ=39)
                                                                                 ctgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTg                >  2:2330931/7‑100 (MQ=39)
                                                                                  tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCGACGCTgg               >  2:725700/6‑100 (MQ=255)
                                                                                  tgcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg               >  2:3370618/6‑100 (MQ=255)
                                                                                  cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg               >  2:3525457/1‑100 (MQ=255)
                                                                                   gcaattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGa              >  2:129578/5‑100 (MQ=255)
                                                                                   aaGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGa              >  1:3432397/1‑100 (MQ=255)
                                                                                    caattGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAg             >  1:3724661/4‑100 (MQ=255)
                                                                                     ggTTGTGTACATTTGTACACAATTGCAGAGTGAGGAAATGACATGCGTCACGCTAACGATCCGGAACGGAGAAAAAATATTATCCAGGCCACGCTGGAgg            >  1:2236335/1‑100 (MQ=255)
                                                                                      gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc           >  1:1923444/1‑100 (MQ=255)
                                                                                      attGTGTACAAATGTACACAATTGCAGAGCGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc           >  1:4174151/2‑100 (MQ=255)
                                                                                      attGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGc           >  2:379542/2‑100 (MQ=255)
                                                                                       ttGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCg          <  2:2886221/100‑1 (MQ=255)
                                                                                       ttGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCg          <  2:1787659/100‑1 (MQ=255)
                                                                                        tgtgTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg         <  2:2936312/100‑1 (MQ=255)
                                                                                         gtgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt        >  2:2788169/1‑100 (MQ=255)
                                                                                         gtgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt        <  2:1712407/100‑1 (MQ=255)
                                                                                         gtgtACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGt        >  2:4371082/1‑100 (MQ=255)
                                                                                           gtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa      <  1:799429/100‑1 (MQ=255)
                                                                                           gtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa      >  2:4050743/1‑100 (MQ=255)
                                                                                            tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa     <  1:4059473/100‑1 (MQ=255)
                                                                                            tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaa     <  2:1850376/100‑1 (MQ=255)
                                                                                              cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc   <  2:3903000/100‑1 (MQ=255)
                                                                                              cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc   <  1:1847663/100‑1 (MQ=255)
                                                                                              cATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAAc   <  2:1358739/100‑1 (MQ=255)
                                                                                                tttGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCAAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACt  >  1:611044/1‑99 (MQ=255)
                                                                                                |                                                                                                  
CAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACT  >  NC_000913/887296‑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: