Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 515,934 C→T 100% A43V (GCG→GTG)  ybbM → predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb515,9340CT91.2% 89.7 / 1.9 34A43V (GCG→GTG) ybbMpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (19/12);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCC  >  W3110S.gb/515872‑515995
                                                              |                                                             
cTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATa                                                           >  1:3309189/1‑67 (MQ=255)
 ttAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGgcgagtga                                                             >  1:1151065/1‑64 (MQ=255)
 ttAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATaa                                                          >  1:3611290/1‑67 (MQ=255)
          cATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGat                                                <  1:1052351/68‑1 (MQ=255)
           aTAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTAGGGCGAGTGATAATTCAGCTGatta                                              <  1:2322460/69‑1 (MQ=255)
             aaaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGattatt                                            <  1:1302554/69‑1 (MQ=255)
              aaGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTg                                           >  1:3407280/1‑69 (MQ=255)
                     aCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGc                                      <  1:1300147/67‑1 (MQ=255)
                     aCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGc                                      <  1:1642565/67‑1 (MQ=255)
                     aCTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGc                                      <  1:3184727/67‑1 (MQ=255)
                       tGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTAtg                                  <  1:3368142/69‑1 (MQ=255)
                                 gaAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatta                                              >  1:1068924/1‑47 (MQ=255)
                                 gaAAGATATTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGatta                                              >  1:540545/1‑47 (MQ=255)
                                  aaaGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGtat                       <  1:2393452/69‑1 (MQ=255)
                                      atatTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTatat                     >  1:1845716/1‑67 (MQ=255)
                                      atatTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAtt                    >  1:1619503/1‑68 (MQ=255)
                                      atatTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATAtt                    >  1:2819141/1‑68 (MQ=255)
                                        atTCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGc                                      >  1:1858696/1‑48 (MQ=255)
                                        atTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTc                 >  1:2313128/1‑69 (MQ=255)
                                         ttCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGc                                      <  1:3259678/47‑1 (MQ=255)
                                         ttCTCTGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCa                >  1:1992930/1‑69 (MQ=255)
                                           ctctGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTg                             >  1:3138437/1‑54 (MQ=255)
                                           ctctGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGta                        >  1:1521616/1‑59 (MQ=255)
                                           ctctGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAg               <  1:2401454/68‑1 (MQ=255)
                                              tGGAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg           >  1:866907/1‑69 (MQ=255)
                                                gAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTgg          >  1:19692/1‑68 (MQ=255)
                                                gAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTgg          >  1:3710978/1‑68 (MQ=255)
                                                gAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTgg          >  1:1109258/1‑68 (MQ=255)
                                                gAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg           <  1:3413547/67‑1 (MQ=255)
                                                gAGCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTg           <  1:3729313/67‑1 (MQ=255)
                                                  gCGTCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat        <  1:222020/68‑1 (MQ=255)
                                                     tCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCATCGTGgatgat     >  1:3683398/1‑68 (MQ=255)
                                                     tCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgatgat     >  1:969051/1‑68 (MQ=255)
                                                     tCGGGCGAGTGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGgat        >  1:1225040/1‑65 (MQ=255)
                                                       gggCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGcc  >  1:3714181/1‑69 (MQ=255)
                                                              |                                                             
CTTAATTAGCCATAAAGAAAAACTGGCGCTGGAGAAAGATATTCTCTGGAGCGTCGGGCGAGCGATAATTCAGCTGATTATTGTCGGCTATGTGCTGAAGTATATTTTCAGCGTGGATGATGCC  >  W3110S.gb/515872‑515995

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: