Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,511,637 A→G intergenic (+101/‑438) acrD → / → yffB aminoglycoside/multidrug efflux system/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,511,6370AG95.2% 239.4 / 1.0 62intergenic (+101/‑438)acrD/yffBaminoglycoside/multidrug efflux system/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (37/22);  total (40/22)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTAT  >  REL606/2511564‑2511709
                                                                         |                                                                        
gtCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttgct                                                                        >  1:473771/1‑76 (MQ=255)
 tCGCGTTGTTCAACTTTCCGATTAACAACCTGCTCAGCGGGTTCTTGCTGTTTGTCCTTTGTCTCAGGAAttgctt                                                                       <  2:956992/76‑1 (MQ=255)
 tCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttgctt                                                                       >  1:1626004/1‑76 (MQ=255)
  cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttgctta                                                                      <  1:1185922/76‑1 (MQ=255)
  cGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAAttactta                                                                      >  1:621064/1‑76 (MQ=255)
    ccTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACg                                                                    >  2:830815/1‑76 (MQ=255)
      tttttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGaa                                                                  <  2:1344522/76‑1 (MQ=255)
       ttttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAg                                                                 >  2:1184946/1‑76 (MQ=255)
       ttttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAg                                                                 >  2:1213641/1‑76 (MQ=255)
        tttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGc                                                                <  1:656185/76‑1 (MQ=255)
        tttCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGc                                                                <  2:1611277/76‑1 (MQ=255)
          tCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCAt                                                              <  1:314335/76‑1 (MQ=255)
          tCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCAt                                                              >  2:998414/1‑76 (MQ=255)
            aaCTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATaa                                                            >  1:1210638/1‑76 (MQ=255)
                  ccGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGa                                                      <  1:1543212/76‑1 (MQ=255)
                   cGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGAt                                                     >  2:1114969/1‑76 (MQ=255)
                      ttAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATaaa                                                  <  1:203657/76‑1 (MQ=255)
                      ttAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATaaa                                                  >  1:206934/1‑76 (MQ=255)
                        aaGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGt                                                >  2:609854/1‑76 (MQ=255)
                         aGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTc                                               >  1:1024770/1‑76 (MQ=255)
                            aCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCttt                                            <  1:211126/76‑1 (MQ=255)
                            aCCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCttt                                            >  1:258287/1‑76 (MQ=255)
                             ccTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTc                                           >  1:1060707/1‑76 (MQ=255)
                              cTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTcc                                          <  1:1506881/76‑1 (MQ=255)
                               tGCTCAGCGGGTTCTTGCTGTTTGTACTATGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCa                                         <  2:66716/76‑1 (MQ=255)
                                 cTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGt                                       >  2:379132/1‑76 (MQ=255)
                                  tCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGtt                                      >  2:808705/1‑76 (MQ=255)
                                    aGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcc                                    >  2:482397/1‑76 (MQ=255)
                                    aGCGGGTTCTTGCTGTTTGTACTTGGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcc                                    <  2:246727/76‑1 (MQ=255)
                                      cGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTcccc                                  >  2:1615579/1‑76 (MQ=255)
                                       gggTTCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCa                                 >  1:739179/1‑76 (MQ=255)
                                          ttCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGtt                              >  1:137482/1‑76 (MQ=255)
                                          ttCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGtt                              >  1:1394073/1‑76 (MQ=255)
                                           tCTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTc                             >  2:885562/1‑76 (MQ=255)
                                            cTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCa                            >  2:809556/1‑76 (MQ=255)
                                            cTTGCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCa                            >  2:809571/1‑76 (MQ=255)
                                               gCTGTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGt                         >  1:1141134/1‑76 (MQ=255)
                                                 tGTTTGTACTTGGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTc                       >  1:289867/1‑76 (MQ=255)
                                                  gTTTGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCa                      >  2:1205589/1‑76 (MQ=255)
                                                   tttGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCaa                     >  1:1624840/1‑76 (MQ=255)
                                                    ttGTACTTTGTCTCAGGAATTGCTTCCGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt                    >  2:1642994/1‑76 (MQ=255)
                                                    ttGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt                    >  1:932220/1‑76 (MQ=255)
                                                    ttGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAAt                    >  2:106349/1‑76 (MQ=255)
                                                     tGTACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATc                   >  1:306572/1‑76 (MQ=255)
                                                       tACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCAt                 >  1:638504/1‑76 (MQ=255)
                                                       tACTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCAt                 >  1:1350716/1‑76 (MQ=255)
                                                        aCTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATa                >  1:1602057/1‑76 (MQ=255)
                                                        aCTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATa                >  2:1301236/1‑76 (MQ=255)
                                                         cTTTGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaa               <  1:1575759/76‑1 (MQ=255)
                                                          tttGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaac              <  1:83426/76‑1 (MQ=255)
                                                          tttGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaac              <  1:758409/76‑1 (MQ=255)
                                                          tttGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaac              >  1:760795/1‑76 (MQ=255)
                                                           ttGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaaca             <  1:1616942/76‑1 (MQ=255)
                                                            tGTCTCAGGAATTGCTTGCGAAGCATGACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaacaa            <  2:1183228/76‑1 (MQ=255)
                                                            tGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaacaa            >  1:998392/1‑76 (MQ=255)
                                                            tGTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaacaa            >  2:423809/1‑76 (MQ=255)
                                                             gTCTCAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATaacaac           >  2:626065/1‑76 (MQ=255)
                                                               ctcAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACct         <  1:1092191/76‑1 (MQ=255)
                                                               ctcAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACct         <  2:515489/76‑1 (MQ=255)
                                                               ctcAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACct         <  1:896745/76‑1 (MQ=255)
                                                                 cAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctct       >  1:488001/1‑76 (MQ=255)
                                                                 cAGGAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctct       >  1:487995/1‑76 (MQ=255)
                                                                   ggAATTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACctctct     <  1:487286/76‑1 (MQ=255)
                                                                      aTTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCtt    >  2:294258/1‑74 (MQ=255)
                                                                      aTTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCtt    <  1:294258/74‑1 (MQ=255)
                                                                      aTTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTtat  <  2:698721/76‑1 (MQ=255)
                                                                      aTTGCTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTtat  <  2:874275/76‑1 (MQ=255)
                                                                         |                                                                        
GTCGCCTTTTTCAACTTTCCGATTAAGAACCTGCTCAGCGGGTTCTTGCTGTTTGTACTTTGTCTCAGGAATTACTTACGAAGCATAACTTCGATAAAGTCTTTCCAGTTCCCCAGTTCACGTTCAATCATAACAACCTCTCTTAT  >  REL606/2511564‑2511709

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: