Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,029,729 (C)6→5 coding (191/921 nt) yedA → amino acid exporter for phenylalanine, threonine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,029,7240C.96.5% 110.2 / ‑2.8 29coding (186/921 nt)yedAamino acid exporter for phenylalanine, threonine
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (16/12);  minor base A (1/0);  total (17/12)
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.00e+00

ATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029588‑2029858
                                                                                                                                        |                                                                                                                                      
atatCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccc                                                                                                                                      >  1:64772/1‑139 (MQ=255)
 tatCATTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTccccc                                                                                                                                    <  1:529960/139‑1 (MQ=255)
     aTTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTcccccgct                                                                                                                                 <  1:395446/139‑4 (MQ=255)
         ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  2:646934/1‑139 (MQ=255)
         ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:552972/1‑139 (MQ=255)
         ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:798074/1‑139 (MQ=255)
         ggggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGt                                                                                                                            >  1:140700/1‑139 (MQ=255)
           ggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCTGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTcc                                                                                                                          >  1:363412/1‑139 (MQ=255)
           ggCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTcc                                                                                                                          >  1:657821/1‑139 (MQ=255)
            gCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCg                                                                                                                         <  1:57716/139‑1 (MQ=255)
              tCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCgct                                                                                                                       >  2:692519/1‑139 (MQ=255)
                aaCCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCgctgc                                                                                                                     >  1:210658/1‑139 (MQ=255)
                          gTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCcgc                                                                                                           <  1:692519/139‑1 (MQ=255)
                                cGGATTGGCGTGGAAAGCTGGCCTCCGGTAATGTTGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACACACCCGCTACGTCCGCTGCTCAATGCCGCGCTGAt                                                                                                     >  1:802996/1‑139 (MQ=255)
                                             aaaGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTg                                                                                        <  2:798074/139‑1 (MQ=255)
                                                                    gcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGGTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGt                                                                 >  2:529197/1‑139 (MQ=255)
                                                                    gcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGt                                                                 >  2:284450/1‑139 (MQ=255)
                                                                    gcgggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGt                                                                 >  2:684607/1‑139 (MQ=255)
                                                                      gggcgTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGa                                                               <  2:657821/139‑1 (MQ=255)
                                                                                   cTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACa                                                  <  2:128073/139‑1 (MQ=255)
                                                                                          ccGGTATATTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAAt                                           >  1:524405/1‑139 (MQ=255)
                                                                                                        cTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCAt                             >  2:595314/1‑139 (MQ=255)
                                                                                                          ggCATTTTTGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATcg                           >  1:646200/1‑139 (MQ=255)
                                                                                                                tttGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGccgccg                     <  1:595314/139‑1 (MQ=255)
                                                                                                                  tGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTa                   <  1:87163/139‑1 (MQ=255)
                                                                                                                  tGCTACTGCGCGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTa                   <  1:684607/139‑1 (MQ=255)
                                                                                                                         gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  <  1:554380/53‑1 (MQ=255)
                                                                                                                         gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTgg                                                                                                  >  2:554380/1‑53 (MQ=255)
                                                                                                                         gcgcGGACACAAACTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTg            >  2:320940/1‑139 (MQ=255)
                                                                                                                                   aaaCTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  2:524405/139‑1 (MQ=255)
                                                                                                                                   aaaCTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  2:646200/139‑1 (MQ=255)
                                                                                                                                   aaaCTCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTCGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:529197/139‑1 (MQ=255)
                                                                                                                                   aaaCTCCCCCGCTACGACCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGATAATGGCATGGTGACGGTTGACGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGcc  <  1:620450/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
ATATCATTTGGGGCTCAACCTATTTTGTCATTCGGATTGGCGTGGAAAGCTGGCCTCCGTTAATGATGGCGGGCGTTCGATTCCTGGCAGCCGGTATTTTATTGCTGGCATTTTTGCTACTGCGCGGACACAAACTCCCCCCGCTACGTCCGCTGCTCAATGCCGCGCTGATTGGCCTGTTATTGCTGGCTGTCGGTAATGGCATGGTGACGGTTGCCGAACATCAAAATGTTCCTTCCGGCATCGCCGCCGTAGTGGTTGCAACCGTGCC  >  NC_000913/2029588‑2029858

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