Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,527,924 G→A 7.3% R306C (CGC→TGC)  mrpD_2 ← Na( )/H( ) antiporter subunit D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,527,9240GA 7.3% 166.6 / 3.5 55R306C (CGC→TGC) mrpD_2Na( )/H( ) antiporter subunit D
Reads supporting (aligned to +/- strand):  ref base G (36/15);  new base A (2/2);  total (38/17)
Fisher's exact test for biased strand distribution p-value = 5.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

CACAGTGCAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGC  >  ACB122/2527794‑2528056
                                                                                                                                  |                                                                                                                                    
cacaGTGCAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCCGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAg                                                                                                                                 >  6:27355/1‑136 (MQ=255)
    gTGCAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGc                                                                                                                             >  1:5981/1‑136 (MQ=255)
    gTGCAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGc                                                                                                                             >  4:68944/1‑136 (MQ=255)
      gCAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGc                                                                                                                             >  7:47565/1‑134 (MQ=255)
       cAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAg                                                                                                                          >  1:126044/1‑136 (MQ=255)
       cAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGAGCCGGCCGAACCCAGGATACTGAGTGCCGCCATAAGCTCCAGGCGCTGAGCTGCCAg                                                                                                                          >  1:134992/1‑136 (MQ=255)
         aGGGTCGAGTTGACCAGGTAGAACAGAGCAGCGCCGGTCAGGGCCGCCAGGGCAAAACCCACCGCACCCAGAAGCGTGCCGCCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCAGCCAg                                                                                                                          >  7:85427/1‑134 (MQ=255)
               gAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATa                                                                                                                  <  2:40800/136‑1 (MQ=255)
                aGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATAc                                                                                                                 >  8:65117/1‑136 (MQ=255)
                 gTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGCTAcc                                                                                                                >  5:51316/1‑136 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGCAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCg                                                                                                               >  2:78282/1‑136 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCt                                                                                                                                    >  7:5931/1‑115 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCt                                                                                                                                    <  8:5931/115‑1 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCg                                                                                                               >  2:72572/1‑136 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCg                                                                                                               >  7:190219/1‑136 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCg                                                                                                               >  7:156399/1‑136 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCg                                                                                                               >  7:148381/1‑136 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCg                                                                                                               >  5:112746/1‑136 (MQ=255)
                  ttGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCg                                                                                                               >  5:124467/1‑136 (MQ=255)
                        aGGTATAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCa                                                                                                         <  6:25295/136‑1 (MQ=255)
                        aGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCa                                                                                                         >  7:206213/1‑136 (MQ=255)
                                  gcgcagcgcCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGc                                                                                                              <  6:127725/121‑1 (MQ=255)
                                  gcgcagcgcCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGc                                                                                                              >  5:127725/1‑121 (MQ=255)
                                  gcgcagcgcCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCACTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGt                                                                                               >  7:161378/1‑136 (MQ=255)
                                  gcgcagcgcCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCACTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGt                                                                                               >  3:95955/1‑136 (MQ=255)
                                      agcgcCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGAcc                                                                                           >  5:194669/1‑136 (MQ=255)
                                      agcgcCGGTCAGGGCCGCCTGCGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGCGCTGCCAGGAAAGATACCGCTGCCACGGCCAGGGTGAc                                                                                            >  5:77550/1‑135 (MQ=255)
                                                aGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGccgccg                                                                                  >  7:86704/1‑135 (MQ=255)
                                                aGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGccgccg                                                                                  >  1:52242/1‑135 (MQ=255)
                                                          gggCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCa                                                                       >  8:158740/1‑136 (MQ=255)
                                                           ggCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCAt                                                                      <  7:158740/136‑1 (MQ=255)
                                                               aaaCCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGCCCAGGCCGCCGTAAACCAGCCATTcgc                                                                  >  7:61783/1‑136 (MQ=255)
                                                                        gCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAAt                                                         >  8:10984/1‑136 (MQ=255)
                                                                        gCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAAt                                                         >  3:127286/1‑136 (MQ=255)
                                                                         cACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATg                                                        >  5:202507/1‑136 (MQ=255)
                                                                           cccAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATgcg                                                      <  8:156399/136‑1 (MQ=255)
                                                                                     tGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCACTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGcccc                                            <  8:161378/136‑1 (MQ=255)
                                                                                     tGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCACTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGcccc                                            <  4:95955/136‑1 (MQ=255)
                                                                                        cggccgAAACCAGGATACTGAGTGCCGCCATACGATCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAACCCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGc                                         >  4:113452/1‑136 (MQ=255)
                                                                                          gccgAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCt                                       <  8:190219/136‑1 (MQ=255)
                                                                                            cgAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGc                                     <  4:127286/136‑1 (MQ=255)
                                                                                            cgAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGc                                     <  6:202507/136‑1 (MQ=255)
                                                                                                aCCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCg                                  >  1:128191/1‑135 (MQ=255)
                                                                                                             gTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGt                    >  7:106782/1‑136 (MQ=255)
                                                                                                              tGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTc                   <  1:72572/136‑1 (MQ=255)
                                                                                                              tGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGACCGGAGAACAGCAGCGTc                   <  6:175517/136‑1 (MQ=255)
                                                                                                               gccgccATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCAcc                                                                                                       <  2:37785/51‑1 (MQ=255)
                                                                                                               gccgccATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCAcc                                                                                                       >  1:37785/1‑51 (MQ=255)
                                                                                                                 cgacATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCca                 >  6:2966/4‑135 (MQ=255)
                                                                                                                   ccATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAg              >  6:87797/1‑136 (MQ=255)
                                                                                                                    cATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAgg             >  6:148754/1‑136 (MQ=255)
                                                                                                                         gCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCTGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCgcag        >  3:161302/1‑136 (MQ=255)
                                                                                                                          cTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCgcagc       <  2:128191/136‑1 (MQ=255)
                                                                                                                          cTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCgcagc       <  5:2966/136‑1 (MQ=255)
                                                                                                                               ggCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGc  <  7:26608/136‑1 (MQ=255)
                                                                                                                                gcgcTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGCCCGCCTTAAACCAGCCATACGCCGCCTAAATGCGCCGAGGCCCCGGCCTGCCCGTATAACAGCAGCGTCCACAGGCGCAGCACGGc  <  6:112746/135‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                    
CACAGTGCAAGGGTCGAGTTGACCAGGTAGAACAGCGCAGCGCCGGTCAGGGCCGCCTGGGCAAAACCCACCGCACCCAGCAGCGTGCCGGCCGAAACCAGGATACTGAGTGCCGCCATACGCTCCAGGCGCTGAGCTGCCAGGATCGATACCGCTGCCACCGCCAGGGTGACCAGGCCGCCGTAAACCAGCCATTCGCCGCCGAAATGCGCCGAGGCCCCGGCCTGCCCGGAGAACAGCAGCGTCCACAGGCGCAGCACGGC  >  ACB122/2527794‑2528056

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