Predicted mutation
evidence seq id position mutation annotation gene description
RA USA300TCH1516_ALE 1,902,155 C→T A105A (GCC→GCT arsB → Arsenical pump membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,902,1550CT92.7% 103.5 / 0.6 41A105A (GCC→GCTarsBArsenical pump membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (21/17);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 2.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

TTTTGTCTCAGTCATTCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCCATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGCTGA  >  USA300TCH1516_ALE/1902023‑1902289
                                                                                                                                    |                                                                                                                                      
ttttGTCTCAGTCATTCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTtgc                                                                                                                                >  2:193925/1‑141 (MQ=255)
    gTCTCAGTCATTCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCa                                                                                                                            <  2:133276/141‑1 (MQ=255)
    gTCTCAGTCATTCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCa                                                                                                                            <  2:403406/141‑1 (MQ=255)
    gTCTCAGTCATTCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCa                                                                                                                            <  2:401397/141‑1 (MQ=255)
           tCATTCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCg                                                                                                                     >  1:9830/1‑141 (MQ=255)
              ttCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCaa                                                                                                                  >  2:355318/1‑141 (MQ=255)
                 ttATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAatt                                                                                                               >  2:283300/1‑140 (MQ=255)
                      tCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGgc                                                                                                          >  2:76493/1‑141 (MQ=255)
                        aTTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGgcgc                                                                                                        <  1:195750/141‑1 (MQ=255)
                                 aGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAAt                                                                                               <  1:335035/141‑1 (MQ=255)
                                  gATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATc                                                                                              >  1:259179/1‑141 (MQ=255)
                                  gATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATc                                                                                              >  1:139901/1‑141 (MQ=255)
                                  gATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATc                                                                                              >  1:133700/1‑141 (MQ=255)
                                           ggATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCt                                                                                     >  1:193417/1‑141 (MQ=255)
                                                  tCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTAt                                                                              >  2:258138/1‑141 (MQ=255)
                                                  tCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTAt                                                                              >  2:19923/1‑141 (MQ=255)
                                                  tCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTAt                                                                              >  2:150219/1‑141 (MQ=255)
                                                  tCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGATGCTATTGTTGCTGCATTTTTCCCTTTTGGTTGCGCAGCGTTACTCTTAACGGCTGTAGTAt                                                                              >  2:219614/1‑141 (MQ=255)
                                                                  tCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTaaaaa                                                              <  1:224331/141‑1 (MQ=255)
                                                                    acaTGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAata                                                            <  2:9830/141‑1 (MQ=255)
                                                                     catgcttcatgcttCAAAAGGAAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCCATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAatat                                                           >  2:265582/1‑141 (MQ=255)
                                                                             atgcttcaAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTa                                                   >  1:324501/1‑141 (MQ=255)
                                                                             atgcttcaAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTa                                                   >  1:204356/1‑141 (MQ=255)
                                                                                    aaaaGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTTTGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCg                                            <  1:316220/141‑1 (MQ=255)
                                                                                    aaaaGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCg                                            <  1:332697/141‑1 (MQ=255)
                                                                                       aGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGTGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGc                                         >  1:183050/1‑141 (MQ=255)
                                                                                       aGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGc                                         <  2:270175/141‑1 (MQ=255)
                                                                                       aGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGc                                         <  2:341720/141‑1 (MQ=255)
                                                                                           aaTGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCATATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCata                                     >  2:288476/1‑141 (MQ=255)
                                                                                                    aaaaTGTTCGTTTATATCATATTATTGGGTGCCATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCtttt                            >  1:341729/1‑141 (MQ=255)
                                                                                                    aaaaTGTTCGTTTATATCATATTATTGGGTGCCATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCtttt                            >  2:389232/1‑141 (MQ=255)
                                                                                                        tGTTCGTTTATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTAGCACCAGATGGCGCAGCGTTAATATTTACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTttatta                        <  1:143798/141‑1 (MQ=255)
                                                                                                              tttATATCATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGa                  >  1:224878/1‑141 (MQ=255)
                                                                                                                     cATATTATTGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGttt           <  2:316766/141‑1 (MQ=255)
                                                                                                                           attGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGc     >  1:42629/1‑141 (MQ=255)
                                                                                                                           attGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGc     >  1:129081/1‑141 (MQ=255)
                                                                                                                           attGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGc     >  2:197761/1‑141 (MQ=255)
                                                                                                                            ttGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGCt    <  2:279513/141‑1 (MQ=255)
                                                                                                                            ttGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGCt    <  2:275942/141‑1 (MQ=255)
                                                                                                                            ttGGGTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGCt    <  2:259647/141‑1 (MQ=255)
                                                                                                                              gggTGCTATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGCTGa  <  1:122871/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TTTTGTCTCAGTCATTCTTATTTCATTAATATTAGATAAAGTTGGATTATTCGAATGGTCAGCTATTCACATGCTTCATGCTTCAAAAGGCAATGGTTTAAAAATGTTCGTTTATATCATATTATTGGGTGCCATTGTTGCTGCATTTTTCGCAAATGATGGCGCAGCGTTAATCTTAACGCCTATTGTATTAGCGATGGTTAAAAATATAGGTTTTAGTAAGCGGGCCATATTCCCTTTTATTATTGCGAGTGGTTTTATAGCTGA  >  USA300TCH1516_ALE/1902023‑1902289

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