Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,389,320 T→C 16.7% N358S (AAC→AGC)  yfbK ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,389,3200TC16.7% 66.3 / 8.5 30N358S (AAC→AGC) yfbKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (12/13);  new base C (3/2);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

AGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTG  >  W3110S.gb/2389253‑2389382
                                                                   |                                                              
agagagGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGt                                                             <  1:1302804/71‑1 (MQ=255)
       gTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGCTGTAATTGc                                                       >  1:1625547/1‑70 (MQ=255)
          tCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTAtt                                                   <  1:1137817/71‑1 (MQ=255)
          tCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTAtt                                                   <  1:1438493/71‑1 (MQ=255)
               gTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGCTGTAATTGCTATTCCCCa                                              >  1:1735127/1‑71 (MQ=255)
                        gTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACg                                     <  1:533183/71‑1 (MQ=255)
                        gTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACg                                     <  1:1337508/71‑1 (MQ=255)
                              gTTACCAACATCGGCAATTCGCACCATCATTGCCTCGCTGTAATTGCTATTCCCCACGCCAAACGTCGACa                               <  1:584408/71‑1 (MQ=255)
                                 aCCAACATCGGCAATTCGCACCATCATTGCCTCGtt                                                               >  1:1612551/1‑36 (MQ=255)
                                 aCCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACagag                            >  1:1342299/1‑71 (MQ=255)
                                  ccAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACagag                            <  1:1696029/70‑1 (MQ=255)
                                  ccAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACagag                            <  1:1694052/70‑1 (MQ=255)
                                     aCATCGGCAATTCGCAGCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTaa                         >  1:454061/1‑70 (MQ=255)
                                       aTCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTaa                         >  1:1146318/1‑68 (MQ=255)
                                       aTCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTaa                         >  1:329122/1‑68 (MQ=255)
                                       aTCGGCAATTCGCACCATCAATGCCTCGTTGTAATTGCTATTCCCAACGCCAAACGTCGACAGAGTaa                         >  1:551237/1‑68 (MQ=255)
                                        tCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAacaca                     <  1:1718047/71‑2 (MQ=255)
                                         cGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTaa                         >  1:1054080/1‑66 (MQ=255)
                                             aaTTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCCCGCCAAACGTCGACAGAGTAaccc                      >  1:862633/1‑63 (MQ=255)
                                             aaTTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACAcc                     >  1:1132814/1‑66 (MQ=255)
                                                 cGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCg            >  1:1440928/1‑71 (MQ=255)
                                                 cGCACCATCATTGCCTCGTTGTAATTCCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTccc             <  1:262148/70‑1 (MQ=255)
                                                       atcatTGCCTCGCTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGttt      <  1:982641/71‑1 (MQ=255)
                                                          atTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTaa                         >  1:1014602/1‑49 (MQ=255)
                                                          atTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTaa                         >  1:1771379/1‑49 (MQ=255)
                                                          atTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGtttttt   >  1:1442294/1‑71 (MQ=255)
                                                          atTGCCTCGCTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTaa                         >  1:1403179/1‑49 (MQ=255)
                                                           ttGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTc               <  1:554608/58‑1 (MQ=255)
                                                           ttGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTc               <  1:98286/58‑1 (MQ=255)
                                                           ttGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTg  <  1:1003590/71‑1 (MQ=255)
                                                             gCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAac                        <  1:844330/47‑1 (MQ=255)
                                                                   ttGTAATTGCTATTCCACACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGttt      <  1:449973/59‑1 (MQ=255)
                                                                   |                                                              
AGAGAGGGTATCAATGTAGCTGTAGTTGCCGTTACCAACATCGGCAATTCGCACCATCATTGCCTCGTTGTAATTGCTATTCCCCACGCCAAACGTCGACAGAGTAACACCAGACTCCCGCTGTTTTTTG  >  W3110S.gb/2389253‑2389382

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