Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,044,421 A→C 31.7% C122W (TGT→TGG U712_15635 ← Putative aminotransferase YugH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,044,4210AC31.7% 30.9 / 5.3 34C122W (TGT→TGGU712_15635Putative aminotransferase YugH
Reads supporting (aligned to +/- strand):  ref base A (15/8);  new base C (2/9);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 2.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAATGGCCTTTGTTTTTTCAGTGACTGCCGCTTCAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTGGCGAGTAGCTAAGGTCAAACCTGTTGCTCAAAT  >  CP006881/3044286‑3044557
                                                                                                                                        |                                                                                                                                        
gagctgggCCTTGGTTTTTTCTGCGACGGCCGCTTCAAAATCCGCAGCCGTTGCTTTAAATCCTTTGTCTGCGGTGTTATGACCGGCACCCGGCAACCCGCCTGCCAAACAACCACGCGCATGGACCCCAACAAAcca                                                                                                                                        <  4:99429/132‑1 (MQ=255)
   aaTGGCCTGTGTTTTTTCAGTGACTGCCGCTTCAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAGCCGGGATTCCGCCTGCCAAAGAAACAAGCGCATCGTACTCAACAAAACACGGc                                                                                                                                    >  2:9177/1‑139 (MQ=255)
       gCCTTTGTTTTTTCAGTGACTGCCGCTTCAAAATCCGCAGCCGTTGCTTTAAATCCTTTGGCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCa                                                                                                                                 <  7:227837/138‑1 (MQ=255)
        ccTTTGTTTTTTCAGTGACTGCCGCTTCAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTGTGCCGTGGTATGACCGTGAACCGGGATCCCGCCTGCCAAAGAACCACGCGCAGCGTACGCACCAAACCACGGCTCAg                                                                                                                                >  3:150261/1‑138 (MQ=255)
        ccTTTGTTTTTTCAGTGACTGCCGCTTCAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACCCGGCTCAg                                                                                                                                >  2:138611/1‑138 (MQ=255)
          tttGTTTTTTCAGTGACTGCCGCTTCAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGaa                                                                                                                             >  7:126984/1‑139 (MQ=255)
                  tacAGTGACTGCCGCTTCAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAACCACGGCTCAGGAATGATGACt                                                                                                                     >  6:117850/3‑139 (MQ=255)
                     aGTGACTGCCGCTTAAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTGTGCCGTTGTATGAACGGGACCCGGGATACCGCCTGCCAAACAAACAAGCGCATCGTACGCAAAAAAACACGGCTCAGGAATGATGACTTc                                                                                                                   <  8:161067/138‑1 (MQ=255)
                           tGCCGATTCAAAATCAGCAGGCGTTGATTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAATAAACAACCGGATAGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCt                                                                                                            <  8:243894/139‑1 (MQ=255)
                                cTTCAAAATCTGCAGCCGTTGCTTTAATTCCTTGGTCTGCCTTTGTATCAAGTGTACCCGGGATACCGCCGGCAAAAGAAACAAGCGCATCGTCCGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGAtt                                                                                                       <  8:210722/139‑1 (MQ=255)
                                          agccgcGGTTGATTTAAAGAATTTGTCTGCCGTTGTATGAAGGTGCACCGGCAAACCGCCTCCCAAAGAAACACGCGCATCGTCCGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGgcg                                                                                             <  2:216205/139‑1 (MQ=255)
                                              gcgtttGGTTTAAATCGTTTGTGTGCGTGTGTATGATGGTGACCCGGGATCCCGCCTGCAAACTAAACACGCGAATCGAACTCAAAAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGa                                                                                         <  8:83543/135‑1 (MQ=255)
                                                   tGCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATTCCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGAAGACTTCCTCCCCTGGAGTCACAATGGCGCGGTTCGc                                                                                     >  8:307794/1‑138 (MQ=255)
                                                    gCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCa                                                                                    <  7:239837/138‑1 (MQ=255)
                                                    gCTTTAAATCCTTTGGCTGCCGTTGTATGAACGTGAACCGGTATACCGCCTCCCAAAGAACCAAGCGCATCGTACGCAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCa                                                                                    <  3:274272/138‑1 (MQ=255)
                                                    gCTTAAAATCTTTTGTGTGCGGTTTAATGACCGTGACCCGGGAACCCGCCTCCAACAGAAACAACCGCATCGAACGAAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCaa                                                                                   <  5:311744/139‑1 (MQ=255)
                                                         aaaTCCTTTGTCTGCCGTTTTATCAACGTGACCCGGGATCCCGCCTGCAAAACAACCACGCGCAGCGACCCCAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATcc                                                                              <  5:338290/139‑1 (MQ=255)
                                                               tttGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGCCTCAGGAATGATGACTTCCTAACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGcc                                                                        >  5:317455/1‑139 (MQ=255)
                                                                   gcTGCCGGTGTATGACCGGCACCCGGGATCCCGCCTGCCAAACAAACACGCGCACCGCCCGCAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGa                                                                     <  4:150261/137‑1 (MQ=255)
                                                                    gggCCGTTGTAGGCACGTGAACCGGGATCCCGCCCGCCAAAGAAACCAGCGACTCGTCCGCAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACt                                                                   <  3:271933/137‑1 (MQ=255)
                                                                          ttGTCTGACCGTGTACCGCGATGCCGCCTCCCAAAGAACCACGCTCCTCGAACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGcccc                                                             <  8:294893/139‑1 (MQ=255)
                                                                             tATCAACGTGACCCGGGAACCCGCCTGCCAAACAAACACGCGCATCGTACCAAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCa                                                           <  7:172134/138‑1 (MQ=255)
                                                                                     tGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTcaca                                                   >  1:32213/1‑138 (MQ=255)
                                                                                        aaCCGGGATACCGCCTGCCAAAGAGACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGACTGAATTGCCCCCACTGTcacatt                                                 >  2:153682/1‑136 (MQ=255)
                                                                                        aaCCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAAt                                                 >  2:13169/1‑138 (MQ=255)
                                                                                           cGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACCCAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGc                                             <  1:114410/139‑1 (MQ=255)
                                                                                                     ccTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGCTGGCTATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTg                                   >  3:114960/1‑139 (MQ=255)
                                                                                                     ccTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTg                                   >  6:286139/1‑139 (MQ=255)
                                                                                                                        cATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCACCACTGTCACAATCAGCTCATTATCTGGCGAGTAGCTAAggctgtc                >  5:194194/1‑134 (MQ=255)
                                                                                                                         aTCGTACGCAACAAAACACGGCTGAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTGGCGAGTAGCTAAGGTCAAAc               >  2:186969/1‑139 (MQ=255)
                                                                                                                         aTCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTGGCGAGTAGCTAAGGTCAAAc               >  1:13274/1‑139 (MQ=255)
                                                                                                                            gaaCGCAACAACACACCGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTGGCGAGTAGCTAAGGTCAAACCTg            <  4:123676/137‑1 (MQ=255)
                                                                                                                               cGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTGGCGAGTAGCTAAGGTCAAACCTGTTg         >  8:28333/1‑139 (MQ=255)
                                                                                                                                gCAACAAACCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCaa                                                                                   <  1:89659/64‑1 (MQ=255)
                                                                                                                                gCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCaa                                                                                   >  2:89659/1‑64 (MQ=255)
                                                                                                                                       aaCACGGCTCAGGAATGATGACTTCCTCCCCTGGATTCACAATGGCGCGGAGCTCAAGATCCAAAGCCTGACATGCCCCCACTTTCACAATCAGCTCATTATCTGGCGAGTAGCTAAGGTCAACCCTGTTGATCaaag  >  2:26385/1‑137 (MQ=255)
                                                                                                                                       aaCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAAt                                                 <  7:207980/91‑1 (MQ=255)
                                                                                                                                       aaCACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAAt                                                 >  8:207980/1‑91 (MQ=255)
                                                                                                                                        |                                                                                                                                        
GAGAATGGCCTTTGTTTTTTCAGTGACTGCCGCTTCAAAATCAGCAGCCGTTGCTTTAAATCCTTTGTCTGCCGTTGTATGAACGTGAACCGGGATACCGCCTGCCAAAGAAACAAGCGCATCGTACGCAACAAAACACGGCTCAGGAATGATGACTTCCTCACCTGGATTCACAATGGCGCGGATCGCAATATCCAAAGCCTGACTTGCCCCCACTGTCACAATCAGCTCATTATCTGGCGAGTAGCTAAGGTCAAACCTGTTGCTCAAAT  >  CP006881/3044286‑3044557

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