Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,487,466 C→T 8.5% D77N (GAC→AAC)  uspE_1 ← Universal stress protein E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,487,4660CT 8.5% 144.8 / 3.8 47D77N (GAC→AAC) uspE_1Universal stress protein E
Reads supporting (aligned to +/- strand):  ref base C (22/21);  new base T (2/2);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

CACCGCACACGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATCAGCAGGTGCAGGCGC  >  ACB122/2487332‑2487596
                                                                                                                                      |                                                                                                                                  
cACCGCACACGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTGAGCGGATTGTCCGGGCGGTGGTCCTTGATCGCCAGTTCGCAGCCGTCGGCCTGCTGTACATGGATGAGGGTGTc                                                                                                                                    >  1:137313/1‑135 (MQ=255)
cACCGCACACGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTc                                                                                                                                    >  7:219890/1‑135 (MQ=255)
cACCGCACACGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTc                                                                                                                                    >  3:17402/1‑135 (MQ=255)
      acacGGACATTGGCGCAGCAGCTTCCAGTCGCTCGGGGTGAGCAGGGCCTTTTTAAGCGGATTTTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCGGCTGTACATGGATGATGGTGTCGTGGAg                                                                                                                              >  1:139517/1‑135 (MQ=255)
        acGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGt                                                                                                                            >  8:13209/1‑135 (MQ=255)
        acGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGt                                                                                                                            >  5:96293/1‑135 (MQ=255)
        acGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGt                                                                                                                            <  4:17402/135‑1 (MQ=255)
             cATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCCGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGGAGGTTGTcg                                                                                                                       >  3:142152/1‑135 (MQ=255)
             cATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTcg                                                                                                                       >  6:154376/1‑135 (MQ=255)
             cATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTcg                                                                                                                       <  5:87724/135‑1 (MQ=255)
                  gcgcAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGGGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCa                                                                                                                  <  8:267107/135‑1 (MQ=255)
                   cgcAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCa                                                                                                                  <  3:90551/134‑1 (MQ=255)
                     cagcagCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGcc                                                                                                               >  5:191514/1‑135 (MQ=255)
                        cagcTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCctc                                                                                                            <  7:67902/135‑1 (MQ=255)
                        cagcTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCctc                                                                                                            <  3:11710/135‑1 (MQ=255)
                        cagcTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCctc                                                                                                            <  4:218536/135‑1 (MQ=255)
                                  tCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGGGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATggggg                                                                                                  <  4:122744/135‑3 (MQ=255)
                                  tCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATGGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTg                                                                                                  <  2:232514/135‑1 (MQ=255)
                                     cTGGGGGTGAGCAGGGCCGTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCCGGTTGTCGCGCCAGCCCTCGTCATGGGTGACa                                                                                               >  5:220/1‑135 (MQ=255)
                                                 aGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCg                                                                                   >  5:110522/1‑135 (MQ=255)
                                                           tCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCGTCATGGGTGACATTGTCGAAGCCGTCATCATGCa                                                                         >  5:262335/1‑135 (MQ=255)
                                                                  aTTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACAGgct                                                                  >  5:249286/1‑135 (MQ=255)
                                                                  aTTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCGGCTGTACATGGATGATGGTGTCGTGCCGGTTTTCGCGCCAGCCCACGTCATGGGTGACATTGTCGCAGCCCGCATCATGCAACtctct                                                                  >  4:180092/1‑131 (MQ=255)
                                                                    tGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGgctgc                                                                <  3:89296/135‑1 (MQ=255)
                                                                       ccGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctc                                                              <  5:221433/134‑1 (MQ=255)
                                                                           gCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctcagcag                                                         <  1:137655/135‑1 (MQ=255)
                                                                              gtggtCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctcagcaggct                                                      <  7:13209/135‑1 (MQ=255)
                                                                                  tCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctcagcaggctcag                                                   <  5:154376/134‑1 (MQ=255)
                                                                                   ccTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctcagcaggctcagca                                                 <  4:130355/135‑1 (MQ=255)
                                                                                   ccTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctcagcaggctcagca                                                 >  7:100800/1‑135 (MQ=255)
                                                                                        aTCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAgcg                                              <  1:128989/133‑1 (MQ=255)
                                                                                        aTCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAgcg                                              <  6:96293/133‑1 (MQ=255)
                                                                                                ttagCAGCCTTCGGCCTGCTGTAAATGGATGATGGTTTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCTTCATCTTGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATc                                    <  7:165957/132‑1 (MQ=255)
                                                                                                            ggCGTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCGAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTg                          <  3:135206/133‑1 (MQ=255)
                                                                                                            ggCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTg                          <  7:52121/133‑1 (MQ=255)
                                                                                                             gCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTc                        >  3:137994/1‑134 (MQ=255)
                                                                                                                               atgGTGTTGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATc                 >  6:186871/1‑123 (MQ=255)
                                                                                                                               atgGTGTTGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATc                 >  6:190518/1‑123 (MQ=255)
                                                                                                                               atgGTGTTGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATc                 <  5:190518/123‑1 (MQ=255)
                                                                                                                               atgGTGTTGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATc                 <  5:186871/123‑1 (MQ=255)
                                                                                                                               atgGTGTCGTGCAGGTTGTCGCGCCAGCCCTGGTGAAGGGTGACATTGTCGAAGCCGTCCTCAAGCAACTGGCTGCTCAGCAGGCTCCGCAGGGCACGGTGCTGAAGCCGGTTGTCCCACCTCAGCAGGGGCCgg     >  4:175450/1‑135 (MQ=255)
                                                                                                                               atgGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATCAGCAGGTGCAgg     >  4:20220/1‑135 (MQ=255)
                                                                                                                               atgGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGGCATCATGCAACTGGCGGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTAGCCCATCCGCAGGTGCCgg     >  6:37640/1‑135 (MQ=255)
                                                                                                                               atgGTGTCGTGCAGGTTGTCGCGCCAGCCCACTTCATGGGTGACATTGTCGAAGCCGGCATCATGCAACTGGCTGCTCAGCAGGCTCCCCAGAGCACTGTGATCATGCCGTTTGTCGCCCATCAGCAGGTGCAgg     >  2:73746/1‑135 (MQ=255)
                                                                                                                                  gtgtCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctc                                                              <  8:5948/75‑1 (MQ=255)
                                                                                                                                  gtgtCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTgctc                                                              >  7:5948/1‑75 (MQ=255)
                                                                                                                                  gtgtCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATCAGCAGGTGCAGgcgc  >  7:136847/1‑135 (MQ=255)
                                                                                                                                  gtgtCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATCAGCAGGTGCAGgcgc  >  1:222642/1‑135 (MQ=255)
                                                                                                                                  gtgtCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATCAGCACGTGCAGgcgc  >  5:286085/1‑135 (MQ=255)
                                                                                                                                  gtgtCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTCGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATCAGCAGGTGCAGgcgc  >  5:74445/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CACCGCACACGGACATTGGCGCAGCAGCTTCCAGTCGCTGGGGGTGAGCAGGGCCTTTTTCAGCGGATTGTCCGGGCGGTGGTCCTTGATCACCAGTTCGCAGCCTTCGGCCTGCTGTACATGGATGATGGTGTCGTGCAGGTTGTCGCGCCAGCCCTCTTCATGGGTGACATTGTCGAAGCCGTCATCATGCAACTGGCTGCTCAGCAGGCTCAGCAGCGCACTGTGATCATGCCGTTTGTCGCACATCAGCAGGTGCAGGCGC  >  ACB122/2487332‑2487596

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