Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 427,941 T→C 11.6% T357T (ACT→ACC secD → SecYEG protein translocase auxillary subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb427,9410TC11.6% 216.9 / 17.2 86T357T (ACT→ACCsecDSecYEG protein translocase auxillary subunit
Reads supporting (aligned to +/- strand):  ref base T (39/37);  new base C (5/5);  total (44/42)
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

TACAACCAACCGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTAAGAAAGATGCA  >  W3110S.gb/427876‑428010
                                                                 |                                                                     
tACAACCAGCCGCAGGTTAACATCTCGCTCGATAGCGCGGGTGGTAACATCATGTCTAACTTCACCAAGGa                                                                  <  1:22868/71‑1 (MQ=255)
tACAACCAGCCGCAGGTTAACATCTCGCTCGATAGCGCGGGTGGTAACATCATGTCTAACTTCACCAAGGa                                                                  <  1:1969223/71‑1 (MQ=255)
tACAACCAGCCGCAGGTTAACATCTCGCTCGATAGCGCGGGTGGTAACATCATGTCTAACTTCACCAAGGa                                                                  <  1:1274351/71‑1 (MQ=255)
tACAACCAGCCGCAGGTTAACATCTCGCTCGATAGCGCGGGTGGTAACATCATGTCTAACTTCACCAAGGa                                                                  <  1:903108/71‑1 (MQ=255)
tACAACCAACCGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAgg                                                                   >  1:3716776/1‑70 (MQ=255)
  caaccaacCGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGac                                                                 >  1:109330/1‑70 (MQ=255)
  caaccaacCGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGac                                                                 <  1:1619645/70‑1 (MQ=255)
    accaaccGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGacaa                                                               >  1:1749439/1‑70 (MQ=255)
       aaccGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACAt                                                            >  1:876990/1‑70 (MQ=255)
       aaccGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACAt                                                            >  1:2763725/1‑70 (MQ=255)
       aaccGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATc                                                           >  1:405578/1‑71 (MQ=255)
       aaccGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATc                                                           <  1:751270/71‑1 (MQ=255)
       aaccGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATc                                                           <  1:307671/71‑1 (MQ=255)
            cAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTaa                                                                     >  1:1765957/1‑56 (MQ=255)
            cAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCa                                                       >  1:2610595/1‑70 (MQ=255)
               gTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTa                                                                      >  1:1707752/1‑52 (MQ=255)
               gTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAAcc                                                   <  1:1081301/71‑1 (MQ=255)
               gTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAAcc                                                   <  1:566724/71‑1 (MQ=255)
               gTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAAcc                                                   <  1:3444172/71‑1 (MQ=255)
                ttAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAAcc                                                   >  1:3179847/1‑70 (MQ=255)
                 tAACATCTCGCTCGATAGCGCTTGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCg                                                  >  1:3700472/1‑70 (MQ=255)
                 tAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGa                                                 >  1:2422771/1‑71 (MQ=255)
                  aaCATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTa                                                                      <  1:1658315/49‑1 (MQ=255)
                  aaCATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGAt                                                >  1:1775835/1‑71 (MQ=255)
                    cATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATgg                                              >  1:503675/1‑71 (MQ=255)
                     aTCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCaa                                                      <  1:162846/62‑1 (MQ=255)
                     aTCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAAcc                                                   <  1:2713999/65‑1 (MQ=255)
                     aTCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGa                                                 <  1:2797750/67‑1 (MQ=255)
                     aTCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGa                                                 <  1:3856153/67‑1 (MQ=255)
                     aTCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATTGc                                             <  1:1414001/71‑1 (MQ=255)
                     aTCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGc                                             <  1:3573439/71‑1 (MQ=255)
                      tctcGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCa                                                       <  1:768632/60‑1 (MQ=255)
                        tcgctcgATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAAc                                          >  1:603189/1‑71 (MQ=255)
                          gctcgATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCa                                            <  1:3583466/67‑1 (MQ=255)
                          gctcgATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCa                                            <  1:2281790/67‑1 (MQ=255)
                            tcgATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGac                                                                 <  1:1001550/44‑1 (MQ=255)
                            tcgATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGac                                                                 <  1:2281809/44‑1 (MQ=255)
                            tcgATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGac                                                                 <  1:2559475/44‑1 (MQ=255)
                                tAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGacaaca                                                             >  1:234826/1‑44 (MQ=255)
                                tAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtg                                  <  1:3035357/71‑1 (MQ=255)
                                 aGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtgt                                 >  1:1844324/1‑71 (MQ=255)
                                  gcgcTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAAcc                                                   >  1:3425870/1‑52 (MQ=255)
                                        gtggtAACATCATGTCTAACTTCACCAAGGACAACATCGGCAAACCGATGGCAACTCTGTTTGTGGAGTa                           >  1:2556669/1‑70 (MQ=255)
                                          ggtAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACa                         >  1:3110449/1‑70 (MQ=255)
                                           gtAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACaaa                       <  1:1805252/71‑1 (MQ=255)
                                           gtAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACaaa                       <  1:1823909/71‑1 (MQ=255)
                                            tAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGa                                                 <  1:2179082/44‑1 (MQ=255)
                                            tAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAg                      <  1:3128021/71‑1 (MQ=255)
                                             aaCATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCaa                                           >  1:935256/1‑49 (MQ=255)
                                             aaCATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCaa                                           >  1:1196425/1‑49 (MQ=255)
                                             aaCATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtg                                  >  1:3364491/1‑58 (MQ=255)
                                             aaCATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtg                                  >  1:3843115/1‑58 (MQ=255)
                                             aaCATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGa                     >  1:2571893/1‑71 (MQ=255)
                                              aCATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGtt                                    >  1:3056387/1‑55 (MQ=255)
                                              aCATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGtt                                    >  1:3666219/1‑55 (MQ=255)
                                                 tcatGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtgtg                                <  1:1769076/56‑1 (MQ=255)
                                                 tcatGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtgtg                                <  1:1927360/56‑1 (MQ=255)
                                                 tcatGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtgtg                                <  1:1421878/56‑1 (MQ=255)
                                                 tcatGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtgtg                                <  1:2579287/56‑1 (MQ=255)
                                                   atGTCTAACTTCACTAAGGACAACATCGGCAAACCGATgg                                              >  1:1374879/1‑40 (MQ=255)
                                                   atGTCTAACTTCACTAAGGACAACATCGGCAAACCGATg                                               >  1:1401563/1‑39 (MQ=255)
                                                   atGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGtt                                    >  1:3770642/1‑50 (MQ=255)
                                                   atGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGtt                                    >  1:925265/1‑50 (MQ=255)
                                                    tGTCTAACTTCACTAAGGACAACATCGGCAAACCGa                                                 >  1:3812028/1‑36 (MQ=255)
                                                    tGTCTAACTTCACTAAGGACAACATCGGCAAACCGa                                                 >  1:3078780/1‑36 (MQ=255)
                                                    tGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGt              <  1:723080/71‑1 (MQ=255)
                                                    tGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGt              <  1:2873348/71‑1 (MQ=255)
                                                     gTCTAACTTCACCAAGGACAACATCGGCAAACCGATGGCAACTCTGTTTGTGGAGTACAAAGAc                    <  1:2765883/64‑1 (MQ=255)
                                                        taacttaactAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTaag           >  1:1511823/1‑70 (MQ=255)
                                                        tAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGtt                                    <  1:3154222/45‑1 (MQ=255)
                                                        tAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGtt                                    <  1:2939873/45‑1 (MQ=255)
                                                        tAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGtt                                    <  1:2885395/45‑1 (MQ=255)
                                                        tAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTtgtg                                >  1:2011155/1‑49 (MQ=255)
                                                        tAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTaag           >  1:1181367/1‑70 (MQ=255)
                                                         aaCTTCACCAAGGACAACATCGGCAAACCGATGGCAACTCTGTTTGTGGAGTACAAAGACAGCGGTaaga          >  1:2567949/1‑70 (MQ=255)
                                                         aaCTTCACCAAGGACAACATCGGCAAACCGATGGCAACTCTGTTTGTGGAGTACAAAGACAGCGGTaaga          >  1:2453267/1‑70 (MQ=255)
                                                         aaCTTCACCAAGGACAACATCGGCAAACCGATGGCAACTCTGTTTGTGGAGTACAAAGACAGCGGTaaga          >  1:1563661/1‑70 (MQ=255)
                                                         aaCTTCACCAAGGACAACATCGGCAAACCGATGGCAACTCTGTTTGTGGAGTACAAAGACAGCGGTaag           >  1:2709605/1‑69 (MQ=255)
                                                          aCTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTa                           >  1:2528803/1‑52 (MQ=255)
                                                          aCTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTa                           >  1:1056706/1‑52 (MQ=255)
                                                          aCTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACa                         >  1:3383671/1‑54 (MQ=255)
                                                          aCTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACa                         >  1:794609/1‑54 (MQ=255)
                                                          aCTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACa                         >  1:1503514/1‑54 (MQ=255)
                                                               aCTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTAAGAAAGATGc   <  1:932340/71‑1 (MQ=255)
                                                                cTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTAAGAAAGATg    <  1:2926413/69‑1 (MQ=255)
                                                                cTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTAAGAAAGATGCa  <  1:748544/71‑1 (MQ=255)
                                                                cTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTAAGAAAGATGCa  <  1:1699013/71‑1 (MQ=255)
                                                                 |                                                                     
TACAACCAACCGCAGGTTAACATCTCGCTCGATAGCGCTGGTGGTAACATCATGTCTAACTTCACTAAGGACAACATCGGCAAACCGATGGCAACCCTGTTTGTGGAGTACAAAGACAGCGGTAAGAAAGATGCA  >  W3110S.gb/427876‑428010

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: