Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,332,615 T→A 10.2% intergenic (‑25/+102) yhbE ← / ← rpmA conserved inner membrane protein/50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,332,6150TA10.2% 155.8 / 10.2 59intergenic (‑25/+102)yhbE/rpmAconserved inner membrane protein/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base A (2/4);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCG  >  W3110S.gb/3332548‑3332684
                                                                   |                                                                     
tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGTATTTg                                                                      <  1:1256846/69‑1 (MQ=255)
tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    >  1:3113756/1‑71 (MQ=255)
tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    <  1:2345907/71‑1 (MQ=255)
tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTAGCg                                                                    <  1:3096505/71‑1 (MQ=255)
     aGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGAtt                                                                <  1:486623/70‑1 (MQ=255)
     aGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTa                                                               <  1:2940664/71‑1 (MQ=255)
     aGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTa                                                               <  1:2377492/71‑1 (MQ=255)
     aGCGCCAAAAGAATGCCAATGCCTGCCTGCAGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGAtt                                                                >  1:181026/1‑70 (MQ=255)
      gcgcCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    >  1:1096682/1‑65 (MQ=255)
       cgcCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAAc                                                             >  1:4239466/1‑71 (MQ=255)
          cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    <  1:535730/61‑1 (MQ=255)
          cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    <  1:2993129/61‑1 (MQ=255)
          cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTc                                                          >  1:909873/1‑71 (MQ=255)
          cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTc                                                          >  1:368075/1‑71 (MQ=255)
          cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTc                                                          <  1:2280380/71‑1 (MQ=255)
          cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTAGCGGATTAACGt                                                           <  1:3935389/70‑1 (MQ=255)
           aaaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    >  1:3440685/1‑60 (MQ=255)
            aaaGAATGCCAATGCCTGCCTTCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    <  1:610092/59‑1 (MQ=255)
            aaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGa                                                                  <  1:2895123/61‑1 (MQ=255)
            aaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCg                                                         >  1:3661858/1‑70 (MQ=255)
            aaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCg                                                         >  1:1685255/1‑70 (MQ=255)
             aaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGc                                                        >  1:3337297/1‑70 (MQ=255)
                aatgccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCa                                                    >  1:1712114/1‑71 (MQ=255)
                      aatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCa                                              >  1:2660691/1‑71 (MQ=255)
                      aatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCa                                              <  1:1849985/71‑1 (MQ=255)
                       atgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAt                                             >  1:2031929/1‑71 (MQ=255)
                        tgcctgccTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTg                                                                      <  1:3211375/45‑1 (MQ=255)
                         gcctgcctGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGc                                                                     >  1:578543/1‑45 (MQ=255)
                          cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGc                                                                     >  1:644329/1‑44 (MQ=255)
                          cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGc                                                                     >  1:1209174/1‑44 (MQ=255)
                          cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCa                                                       >  1:554700/1‑58 (MQ=255)
                          cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtttt                                          <  1:3273870/71‑1 (MQ=255)
                          cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtttt                                          >  1:3019187/1‑71 (MQ=255)
                              cctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCa                                                       <  1:3805439/54‑1 (MQ=255)
                              cctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTa                                        <  1:639239/69‑1 (MQ=255)
                               ctgctgCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCg                                     <  1:4094786/71‑1 (MQ=255)
                                tgctgcTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCAt                                                   <  1:3828803/56‑1 (MQ=255)
                                   tgctTCATACCGTATTCCCTGCAAAGCGCATTTGCg                                                                    <  1:2784629/36‑1 (MQ=255)
                                   tgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATgcgc                                               <  1:335849/57‑1 (MQ=255)
                                    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAttttt                                         >  1:3686004/1‑62 (MQ=255)
                                    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAttttt                                         >  1:3817654/1‑62 (MQ=255)
                                    gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTa                                        <  1:644553/63‑1 (MQ=255)
                                     ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCAcc                                                     <  1:1987655/49‑1 (MQ=255)
                                        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCa                                                    >  1:1232480/1‑47 (MQ=255)
                                        cATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGaa                            >  1:162433/1‑71 (MQ=255)
                                         aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGaa                            <  1:1208901/70‑1 (MQ=255)
                                          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAttttt                                         <  1:1465627/56‑1 (MQ=255)
                                          tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTAc                                       >  1:3491998/1‑58 (MQ=255)
                                             cGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTaa                       <  1:2383746/71‑1 (MQ=255)
                                              gTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTaa                       >  1:1486006/1‑70 (MQ=255)
                                                  tCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGt                                   >  1:2229597/1‑54 (MQ=255)
                                                     cTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTaa                       <  1:3121940/63‑1 (MQ=255)
                                                       gCAAAGCGCATTTGCGGATTAACGTCGCACCATgcgc                                               >  1:872783/1‑37 (MQ=255)
                                                        cAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCATTTTACCGGTTATCGAATGTAAAAAGCCCCGc             <  1:1325304/69‑1 (MQ=255)
                                                          aaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAAc          <  1:1929903/71‑1 (MQ=255)
                                                           aGCGCATTTGCGGATTAACGTCGCACCATTCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACg         <  1:4216414/71‑1 (MQ=255)
                                                           aGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGt                         >  1:829581/1‑55 (MQ=255)
                                                             cgcATTAGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGcc                 >  1:3669599/1‑61 (MQ=255)
                                                              gcATTTGCGGATTAACGTCGCACCATGCGCATTATTACCGGTTATCGAATGTAAAAAGCCCCGCAACgtgt      >  1:3678928/1‑71 (MQ=255)
                                                                   tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCg  <  1:3197283/70‑1 (MQ=255)
                                                                   tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCg  <  1:2782429/70‑1 (MQ=255)
                                                                   tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCg  <  1:1910826/70‑1 (MQ=255)
                                                                   |                                                                     
TGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCG  >  W3110S.gb/3332548‑3332684

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