Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 344,381 G→A 72.4% intergenic (+166/‑247) yahL → / → yahM hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb344,3810GA72.4% 88.8 / 41.9 58intergenic (+166/‑247)yahL/yahMhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base A (22/20);  total (30/28)
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

CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGAC  >  W3110S.gb/344314‑344446
                                                                    |                                                                 
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGt                                                                  >  1:263721/1‑69 (MQ=255)
ccATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGt                                                                  >  1:178921/1‑69 (MQ=255)
   tctcTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                <  1:587919/68‑1 (MQ=255)
   tctcTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                >  1:3575790/1‑68 (MQ=255)
   tctcTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGAttt                                                                >  1:3539795/1‑68 (MQ=255)
    ctctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGttt                                                                >  1:225516/1‑67 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGc                                                             >  1:820921/1‑69 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGc                                                             >  1:1947419/1‑69 (MQ=255)
     tctTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGc                                                             >  1:1676675/1‑69 (MQ=255)
         ttttCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGcc                                                         >  1:1826268/1‑69 (MQ=255)
           ttCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGc                                                          <  1:2988766/66‑1 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaa                                            >  1:431810/1‑68 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaa                                            >  1:1051135/1‑68 (MQ=255)
                       ttgttgTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaa                                            >  1:2845936/1‑68 (MQ=255)
                          ttgtGATCGGCGACACTTCGCTTAAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATaaaa                                          <  1:2220578/68‑1 (MQ=255)
                            gtgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATg                                      <  1:1720384/69‑1 (MQ=255)
                            gtgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATg                                      <  1:1688298/69‑1 (MQ=255)
                            gtgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATg                                      <  1:1969782/69‑1 (MQ=255)
                            gtgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATg                                      <  1:1166779/69‑1 (MQ=255)
                            gtgATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATg                                      <  1:202530/69‑1 (MQ=255)
                                 cGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                  <  1:306709/68‑1 (MQ=255)
                                 cGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                  >  1:2625854/1‑68 (MQ=255)
                                 cGGCGAAACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                  <  1:1776075/68‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                  <  1:1332750/67‑1 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCg                                                    >  1:2959875/1‑49 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCg                                                    >  1:18165/1‑49 (MQ=255)
                                  ggCGACACTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCg                                                    >  1:752479/1‑49 (MQ=255)
                                      acacTTCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCa                            <  1:2192334/69‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                       <  1:1981529/69‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGGTTTCGCACCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                       <  1:820227/69‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCa                                                        >  1:2143689/1‑36 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCa                                                        >  1:2891085/1‑36 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCa                                                        >  1:1893463/1‑36 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCa                                                        <  1:2199425/36‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATg                        <  1:122464/68‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATg                        <  1:882794/68‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATg                        <  1:1770929/68‑1 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                       >  1:675624/1‑69 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                       >  1:2141245/1‑69 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                       >  1:981506/1‑69 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                       >  1:2935419/1‑69 (MQ=255)
                                           tCGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGt                       <  1:2690112/69‑1 (MQ=255)
                                            cGCTTAAAAAAGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGtt                      >  1:2764219/1‑69 (MQ=255)
                                                  aaaaaaGCACCAGTAGTGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 >  1:781875/1‑69 (MQ=255)
                                                    aaaaGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 <  1:95379/67‑1 (MQ=255)
                                                     aaaGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 >  1:3549780/1‑66 (MQ=255)
                                                        gCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAAtgtg           <  1:3402237/69‑1 (MQ=255)
                                                        gCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAAtgtg           <  1:2259747/69‑1 (MQ=255)
                                                         cACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAAtg             >  1:1887495/1‑66 (MQ=255)
                                                               tagtGGTTTCGCAGCCATGCGGTGTATAAAAAATGAtct                                  <  1:88739/39‑1 (MQ=255)
                                                               tagtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 <  1:1240590/56‑1 (MQ=255)
                                                               tagtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 <  1:2678680/56‑1 (MQ=255)
                                                                 gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 <  1:3523810/54‑1 (MQ=255)
                                                                 gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 <  1:1289276/54‑1 (MQ=255)
                                                                 gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGAc  >  1:2999137/1‑69 (MQ=255)
                                                                 gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGAc  >  1:1350734/1‑69 (MQ=255)
                                                                 gtGATTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGAc  >  1:3024556/1‑69 (MQ=255)
                                                                 gtGATTTAGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTtgtg                 <  1:2047237/54‑1 (MQ=255)
                                                                    |                                                                 
CCATCTCTTTTTTCGCCGATATTTTGTTGTGATCGGCGACACTTCGCTTAAAAAAGCACCAGTAGTGGTTTCGCAGCCATGCGGTGTATAAAAAATGATCTCATGCAGATGTTTTGTGAATGTGTTGGTTGAC  >  W3110S.gb/344314‑344446

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: