Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,832,787 T→C 21.6% M304V (ATG→GTG)  astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,832,7870TC21.6% 70.5 / 17.1 37M304V (ATG→GTG) astCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base C (5/3);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGATTACCGCCA  >  W3110S.gb/1832717‑1832851
                                                                      |                                                                
ccATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCAt                                                                  <  1:578080/71‑1 (MQ=255)
           aTTAATAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGgtgtg                                                        <  1:1676169/70‑1 (MQ=255)
            ttAATAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGtt                                                      <  1:1321037/71‑1 (MQ=255)
             tAATAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGtt                                                      <  1:1828227/70‑1 (MQ=255)
                tAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGgtgt                                                         <  1:969850/64‑1 (MQ=255)
                tAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGgtgt                                                         <  1:640367/64‑1 (MQ=255)
                tAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTtgatg                                                  <  1:2217271/71‑1 (MQ=255)
                tAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTtgatg                                                  <  1:2543044/71‑1 (MQ=255)
                   tATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCACCTCTGGTGTGtt                                                      <  1:1776000/64‑1 (MQ=255)
                   tATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCACCTCTGGTGTGtt                                                      <  1:1779254/64‑1 (MQ=255)
                         gACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCACCTCTGGTGTGTTGATGAGCTCCAGc                                         <  1:1516146/71‑1 (MQ=255)
                                  caaACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTtga                                                    >  1:1921578/1‑51 (MQ=255)
                                  caaACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTtga                                                    >  1:1946787/1‑51 (MQ=255)
                                  caaACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTtga                                                    >  1:430783/1‑51 (MQ=255)
                                   aaaCCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCAc                                       >  1:1811983/1‑63 (MQ=255)
                                   aaaCCAGTCGTGACGCTGTTTAACGCCATTAAGCACCTCTGGTGTGTTGATGAGCTCCAGCACTTtgcctg                               >  1:1261106/1‑71 (MQ=255)
                                         gTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTTGATGAGc                                               >  1:118764/1‑49 (MQ=255)
                                          tCGTGACGCTGTTTAACGCCATTAAGCATCTCTGgtg                                                          <  1:1870593/37‑1 (MQ=255)
                                            gTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCAccgcc                       >  1:972262/1‑70 (MQ=255)
                                            gTGACGCTGTTTAACGCCATTAAGCACCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCAccgcc                       >  1:1286806/1‑70 (MQ=255)
                                                     tttAACGCCATTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCACttt                                    <  1:641608/48‑1 (MQ=255)
                                                     tttAACGCCATTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAg               <  1:1432458/69‑1 (MQ=255)
                                                     tttAACGCCATTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCg             <  1:525194/71‑1 (MQ=255)
                                                       tAACGCCATTAAGCACCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGa           >  1:377953/1‑71 (MQ=255)
                                                       tAACGCCATTAAGCACCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGa           >  1:381201/1‑71 (MQ=255)
                                                         aCGCCATTAAGCACCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGCCCAGCGGAt          >  1:2501418/1‑70 (MQ=255)
                                                              aTTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCAc                                       <  1:1288359/36‑1 (MQ=255)
                                                                tAAGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGa                     >  1:41388/1‑52 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTtgcc                                 >  1:548727/1‑39 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTtgcc                                 >  1:863192/1‑39 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTtgcc                                 >  1:2798934/1‑39 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTtgc                                  >  1:15974/1‑38 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGATTACCGCCa  >  1:579172/1‑70 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGATTACCGCCa  <  1:2050181/70‑1 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGATTACCGCCa  >  1:21680/1‑70 (MQ=255)
                                                                 aaGCATCTCTGGTGTGTTGATGAACTCCAGCACTTtgcc                                 >  1:263560/1‑39 (MQ=255)
                                                                   gCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGATTACCGCCa  <  1:657198/68‑1 (MQ=255)
                                                                   gCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGATTACCGCCa  <  1:866504/68‑1 (MQ=255)
                                                                      |                                                                
CCATAGCGGTGATTAATAGTATTAAGACGCTCAACAAACCAGTCGTGACGCTGTTTAACGCCATTAAGCATCTCTGGTGTGTTGATGAGCTCCAGCACTTTGCCTGCCACCGCCGAGGCCAGCGGATTACCGCCA  >  W3110S.gb/1832717‑1832851

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