Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 430,031 G→A 14.3% C68Y (TGC→TAC)  yajD → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb430,0310GA14.3% 186.9 / 20.2 70C68Y (TGC→TAC) yajDconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (6/4);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGAAGAT  >  W3110S.gb/429964‑430101
                                                                   |                                                                      
cacaTTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCa                                                                     >  1:1476713/1‑71 (MQ=255)
    ttGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGAt                                                                 >  1:2434965/1‑71 (MQ=255)
    ttGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGAt                                                                 >  1:2657150/1‑71 (MQ=255)
    ttGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGAt                                                                 >  1:1494803/1‑71 (MQ=255)
     tGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCa                                                                     >  1:790638/1‑66 (MQ=255)
      gATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCa                                                               >  1:281048/1‑71 (MQ=255)
      gATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTACCACGATCa                                                               >  1:1150748/1‑71 (MQ=255)
      gATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTACCACGATCa                                                               >  1:432024/1‑71 (MQ=255)
      gATCAAGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATc                                                                <  1:1560581/70‑1 (MQ=255)
       aTCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCAt                                                              <  1:2402669/71‑1 (MQ=255)
       aTCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCAt                                                              <  1:2326271/71‑1 (MQ=255)
       aTCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCAt                                                              <  1:3821961/71‑1 (MQ=255)
        tCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATg                                                             <  1:3439888/71‑1 (MQ=255)
        tCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATg                                                             <  1:1295963/71‑1 (MQ=255)
          aCGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAg                                                           >  1:1756243/1‑71 (MQ=255)
                aTACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCg                                                     <  1:1408589/71‑1 (MQ=255)
                 tACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGa                                                    <  1:2679041/71‑1 (MQ=255)
                 tACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGa                                                    <  1:823666/71‑1 (MQ=255)
                 tACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGa                                                    <  1:1131568/71‑1 (MQ=255)
                  aCCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGa                                                    >  1:1891425/1‑70 (MQ=255)
                   ccAATAACCCGGAAGATGGCAGTAGCTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  >  1:38122/1‑71 (MQ=255)
                   ccAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCa                                                               >  1:3745941/1‑58 (MQ=255)
                   ccAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCa                                                               >  1:942269/1‑58 (MQ=255)
                   ccAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  >  1:3035570/1‑71 (MQ=255)
                   ccAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  >  1:3798515/1‑71 (MQ=255)
                   ccAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  >  1:3524221/1‑71 (MQ=255)
                   ccAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTTGTCTCTATTGCCACGATCa                                                               >  1:4063583/1‑57 (MQ=255)
                    cAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  <  1:362838/70‑1 (MQ=255)
                     aaTAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAat                                                 >  1:734290/1‑70 (MQ=255)
                       tAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAAtata                                              >  1:2866988/1‑71 (MQ=255)
                       tAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTACCACGATCATGAGCATTCGAAAtata                                              >  1:2860479/1‑71 (MQ=255)
                       tAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTACCACGATCATGAGCATTCGAAAtata                                              >  1:2519533/1‑71 (MQ=255)
                        aaCCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAc                                             >  1:1286856/1‑71 (MQ=255)
                         aCCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAcc                                            >  1:3327924/1‑71 (MQ=255)
                         aCCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAcc                                            >  1:1804118/1‑71 (MQ=255)
                         aCCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTACCACGATCATGAGCATTCGAAATATAcc                                            <  1:3670252/71‑1 (MQ=255)
                                 gATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGa                                    >  1:427783/1‑71 (MQ=255)
                                   tGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGAt                                   >  1:2953966/1‑70 (MQ=255)
                                     gCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAg                                <  1:1982798/71‑1 (MQ=255)
                                      cAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGt                               <  1:244936/71‑1 (MQ=255)
                                       aGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTa                              >  1:2056673/1‑71 (MQ=255)
                                         tAACTGGGAATTGTTGTGTCTCTATTACCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAg                                <  1:2781636/67‑1 (MQ=255)
                                            cTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat                                               <  1:3991740/49‑1 (MQ=255)
                                            cTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat                                               <  1:446286/49‑1 (MQ=255)
                                            cTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAata                                                <  1:1187826/48‑1 (MQ=255)
                                             tGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  <  1:2287924/45‑1 (MQ=255)
                                              gggAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAAtata                                              >  1:2763832/1‑48 (MQ=255)
                                              gggAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAAtata                                              >  1:2577712/1‑48 (MQ=255)
                                              gggAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGt                               >  1:3645140/1‑63 (MQ=255)
                                              gggAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTAc                        >  1:1716646/1‑70 (MQ=255)
                                              gggAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTAc                        >  1:2667459/1‑70 (MQ=255)
                                               ggAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAc                                             >  1:2948467/1‑48 (MQ=255)
                                                gAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat                                               <  1:574045/45‑1 (MQ=255)
                                                   ttgttgTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  <  1:274940/39‑1 (MQ=255)
                                                   ttgttgTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  <  1:3420166/39‑1 (MQ=255)
                                                   ttgttgTGTCTCTATTGCCACGATCATGAGCATTCGaaa                                                  <  1:2474022/39‑1 (MQ=255)
                                                   ttgttgTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGt                  <  1:2506485/71‑1 (MQ=255)
                                                     gttgtGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGTGGATCAGTATGGTACGACCGTTa                >  1:2373184/1‑71 (MQ=255)
                                                     gttgtGTCTCTATTACCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTa                >  1:488053/1‑71 (MQ=255)
                                                     gttgtGTCTCTATTACCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTa                >  1:255347/1‑71 (MQ=255)
                                                         tgtCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCg             <  1:117140/70‑1 (MQ=255)
                                                             tctATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATgg                           >  1:928580/1‑52 (MQ=255)
                                                               tATTACCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAg                                <  1:1641024/45‑1 (MQ=255)
                                                               tATTACCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAg                                <  1:771165/45‑1 (MQ=255)
                                                                 ttGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGt                  <  1:2352918/57‑1 (MQ=255)
                                                                 ttGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAgggg       <  1:1890639/68‑1 (MQ=255)
                                                                 ttGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAgggg       <  1:3162990/68‑1 (MQ=255)
                                                                 ttGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGaa     <  1:3861755/70‑1 (MQ=255)
                                                                 ttGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGaa     <  1:452813/70‑1 (MQ=255)
                                                                 ttGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGaa     <  1:3402899/70‑1 (MQ=255)
                                                                   gCCACGATCGTGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGAAGa   >  1:267906/1‑70 (MQ=255)
                                                                   gCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTa                >  1:3896453/1‑57 (MQ=255)
                                                                   gCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTa                >  1:3480619/1‑57 (MQ=255)
                                                                   gCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTa                >  1:2713012/1‑57 (MQ=255)
                                                                   gCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGAAGa   >  1:495491/1‑70 (MQ=255)
                                                                   gCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGAAGAt  >  1:28632/1‑71 (MQ=255)
                                                                   gCCAAGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGAAGAt  >  1:2612141/1‑71 (MQ=255)
                                                                   |                                                                      
CACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCAGTATGGTACGACCGTTATCGCAGGGGAAGAT  >  W3110S.gb/429964‑430101

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: