Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,022,325 Δ1 bp 14.3% intergenic (‑20/+648) PP_3547 ← / ← emrB putative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,022,3250A.14.3% 102.9 / 12.7 35intergenic (‑20/+648)PP_3547/emrBputative NAD(P)‑binding oxidoreductase with NAD(P)‑binding Rossmann‑fold domain/multidrug efflux transport system ‑ permease subunit
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base . (3/2);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01

ACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCG  >  AE015451/4022100‑4022518
                                                                                                                                                                                                                                 |                                                                                                                                                                                                 
aCGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGAt                                                                                                                                                                                        >  1:764152/1‑236 (MQ=255)
       ccATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGGTGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                                                                                   >  1:453652/1‑235 (MQ=255)
        cATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                                                                                   <  2:975186/234‑1 (MQ=255)
                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGAt                                                                                                                                                                                        <  2:818595/217‑1 (MQ=255)
                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGAAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGAt                                                                                                                                                                                        >  1:818596/1‑217 (MQ=255)
                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGt                                                                                                                                                                       >  1:962575/1‑235 (MQ=255)
                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTc                                                                                                                                                                      >  1:947271/1‑236 (MQ=255)
                   cAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTc                                                                                                                                                                      >  1:59524/1‑236 (MQ=255)
                                     aaCATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCg                                                                                                                                                    >  1:1015379/1‑236 (MQ=255)
                                      aCATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAgg                                                                                                                                                                                          >  2:216508/1‑197 (MQ=255)
                                      aCATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAAGCGCTGGGCCCCAgg                                                                                                                                                                                          >  2:478554/1‑196 (MQ=255)
                                          cGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGa                                                                                                                                                                                         <  1:365374/194‑1 (MQ=255)
                                               aTAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGGCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                                                                                         <  2:931911/205‑1 (MQ=255)
                                                     gAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGGCCGGCCGGGAGCAAGGCCGGAGTGAggg                                                                                                                                           >  1:217693/1‑229 (MQ=255)
                                                      aaaTCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAgg                                                                                                                                            >  2:758227/1‑227 (MQ=255)
                                                                 cGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGAAcc                                                                                                                                       <  2:764151/220‑1 (MQ=18)
                                                                    tGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgg                                                                                                                      <  1:594748/235‑1 (MQ=255)
                                                                                    accacaGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGACGCCGAACGCAGGCCTTGCGGAGGCAGGTGCCGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGa                                                                                                                                                             <  2:549714/178‑1 (MQ=255)
                                                                                        caGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGCGTGAGGGa                                                                                                                                          >  2:9003/1‑195 (MQ=255)
                                                                                          gCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                                                                                   <  2:487572/152‑1 (MQ=255)
                                                                                                    tGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgg                                                                                                                      <  2:962574/203‑1 (MQ=255)
                                                                                                       ttttgatctggcttttgatctggcTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGa                                                                                                                                  <  1:426367/188‑1 (MQ=255)
                                                                                                                    tttgatctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGc                                                                                                                                                      >  1:1087049/1‑155 (MQ=38)
                                                                                                                    tttgatctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGc                                                                                                                                                      <  2:1087048/155‑1 (MQ=39)
                                                                                                                    tttgatctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGC‑‑‑‑‑‑‑‑‑TCCCGGCAGGA‑‑‑CCGC‑‑GCCGAAGGCCCCATCCTGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcagcg                                                                                                                          >  1:399267/1‑165 (MQ=9)
                                                                                                                           tggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCGAGGATGGGGCCTTCGGCGCGGTCCTGGCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgg                                                                                                                      <  1:170675/180‑1 (MQ=25)
                                                                                                                                 ttgatcttgatcttgatcttgatcTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTAGGTTACTTTTTGTCGCGTTTGACAAAAAGTg                                                          <  2:1015378/234‑1 (MQ=14)
                                                                                                                                          atcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATgg                                                                                                         >  2:547335/1‑178 (MQ=31)
                                                                                                                                              tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgg                                                                                                                      <  2:401824/161‑1 (MQ=31)
                                                                                                                                                  cttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCAAGCGCTGAGGCCCCATCCTGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTgacaaaaagtaaccaaaaaacgctgcgctc                                       >  1:897070‑M2/1‑206 (MQ=255)
                                                                                                                                                                 ttCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGc                                           <  1:540936/217‑1 (MQ=31)
                                                                                                                                                                                                cgcttgacgGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTGTTTGGGTACTTTTTGTCGCGTTTggcaaa                                                               <  1:257993‑M2/162‑7 (MQ=255)
                                                                                                                                                                                                 aGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGc                                                                                                >  2:61523/1‑132 (MQ=32)
                                                                                                                                                                                                               gaggCCCGGCAGGCGTTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCg                                                                                                                 <  2:264319/98‑1 (MQ=11)
                                                                                                                                                                                                                          aGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCg  >  2:265291/1‑201 (MQ=38)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                 
ACGTTCTCCATGGGTTTGCCAATGGTTTGGGGAAGGGAACATCGACAATAGGCGAAATCGCGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGACAAAAAGTGACTCGCCGTAAGGGCGAAAAGGTGAGTCAGCGTCGCCATCGCAAACGGACGGTTCG  >  AE015451/4022100‑4022518

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