Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,783,965 C→T 9.2% G13D (GGC→GAC)  yidI ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,783,9650CT 9.2% 237.9 / 10.2 76G13D (GGC→GAC) yidIpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (37/32);  new base T (4/3);  total (41/35)
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

AGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAAAAGCAT  >  W3110S.gb/3783907‑3784033
                                                          |                                                                    
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACt                                                            >  1:2011648/1‑69 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACt                                                            >  1:3478166/1‑69 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACt                                                            >  1:3001710/1‑69 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATg                                                                <  1:2879744/65‑1 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATg                                                                <  1:2230399/65‑1 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACt                                                            >  1:2311487/1‑69 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACt                                                            >  1:2293066/1‑69 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACt                                                            >  1:1841737/1‑69 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACt                                                            >  1:2341720/1‑69 (MQ=255)
aGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACt                                                            >  1:3178181/1‑69 (MQ=255)
      gAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGatttt                                                      >  1:1552717/1‑69 (MQ=255)
        aTGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGccc                                                                    >  1:948446/1‑53 (MQ=255)
           aaaaGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGccc                                                                    >  1:92261/1‑50 (MQ=255)
            aaaGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTg                                                >  1:1628197/1‑69 (MQ=255)
            aaaGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTg                                                >  1:3498456/1‑69 (MQ=255)
             aaGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTg                                                <  1:737551/68‑1 (MQ=255)
             aaGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGTCCAATGAACTGATTTTATTTTg                                                <  1:3180220/68‑1 (MQ=255)
               gTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGc                                             >  1:3615186/1‑69 (MQ=255)
                tGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCt                                            <  1:1123638/69‑1 (MQ=255)
                  gATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          >  1:2979694/1‑69 (MQ=255)
                  gATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:945067/69‑1 (MQ=255)
                    tAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGa                                        >  1:2250794/1‑69 (MQ=255)
                     aaTCCCCATCATTAAAGCGATGGCACCGAACAGCATGccc                                                                    <  1:2619725/40‑1 (MQ=255)
                       tCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCt                                            <  1:1995475/62‑1 (MQ=255)
                        ccccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATAcc                                    <  1:760733/69‑1 (MQ=255)
                         cccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCt                                            <  1:3105886/60‑1 (MQ=255)
                         cccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCt                                            <  1:1287100/60‑1 (MQ=255)
                         cccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCt                                            <  1:2882076/60‑1 (MQ=255)
                         cccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCt                                            <  1:3239927/60‑1 (MQ=255)
                         cccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATAccc                                   <  1:2232037/69‑1 (MQ=255)
                         cccATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATAccc                                   <  1:2829812/69‑1 (MQ=255)
                             tcatTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               >  1:2239441/1‑69 (MQ=255)
                             tcatTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               >  1:3569281/1‑69 (MQ=255)
                             tcatTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               >  1:1852061/1‑69 (MQ=255)
                             tcatAAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATc                                >  1:3657483/1‑68 (MQ=255)
                               atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:972595/56‑1 (MQ=255)
                               atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               <  1:2632697/67‑1 (MQ=255)
                               atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               <  1:143456/67‑1 (MQ=255)
                               atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               >  1:373790/1‑67 (MQ=255)
                               atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               <  1:1603831/67‑1 (MQ=255)
                               atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGa                             >  1:881356/1‑69 (MQ=255)
                               atTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGa                             >  1:316494/1‑69 (MQ=255)
                                 tAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGAcc                           >  1:553431/1‑69 (MQ=255)
                                 tAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGAcc                           >  1:2138343/1‑69 (MQ=255)
                                  aaaGCGATGGCACCGAACAGCATGCCCAATGAACTGatttta                                                     >  1:121694/1‑42 (MQ=255)
                                  aaaGCGATGGCACCGAACAGCATGCCCAATGAACTGatttta                                                     >  1:2975270/1‑42 (MQ=255)
                                  aaaGCGATGGCACCGAACAGCATGCCCAATGAACTGatttta                                                     >  1:3462947/1‑42 (MQ=255)
                                  aaaGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:1144089/53‑1 (MQ=255)
                                  aaaGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCt                               >  1:1181008/1‑64 (MQ=255)
                                  aaaGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTg                              <  1:728637/65‑1 (MQ=255)
                                  aaaGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTg                              >  1:3426301/1‑65 (MQ=255)
                                          ggCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAg                  >  1:3151172/1‑69 (MQ=255)
                                          ggCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAg                  >  1:2730537/1‑69 (MQ=255)
                                          ggCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAg                  >  1:3212836/1‑69 (MQ=255)
                                          ggCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAg                  >  1:2371312/1‑69 (MQ=255)
                                             aCCGAACAGCATGTCCAATGAACTGATTTTATTTTg                                                <  1:1535760/36‑1 (MQ=255)
                                             aCCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCg                <  1:2304548/68‑1 (MQ=255)
                                             aCCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCg                <  1:3465123/68‑1 (MQ=255)
                                             aCCGAACAGCATGCCCAATGAACTGATTATATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCg                <  1:3353999/68‑1 (MQ=255)
                                             aCCGAACAGCATGCCCAATGAACTGATTATATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCg                <  1:3358338/68‑1 (MQ=255)
                                              ccGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCAt                                 >  1:87181/1‑50 (MQ=255)
                                              ccGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCAt                                 >  1:2017371/1‑50 (MQ=255)
                                              ccGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCg                <  1:2582939/67‑1 (MQ=255)
                                                gAACAGTATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGAt             >  1:477934/1‑68 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:2478120/39‑1 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:2638098/39‑1 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:223677/39‑1 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          >  1:41233/1‑39 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          >  1:3680749/1‑39 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          >  1:3149061/1‑39 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:3017908/39‑1 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGAt             >  1:15978/1‑68 (MQ=255)
                                                gAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACATCCTTGTAGCGGAt             >  1:190951/1‑68 (MQ=255)
                                                 aaCAGCATGTCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATg            >  1:3103418/1‑68 (MQ=255)
                                                 aaCAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:1834435/38‑1 (MQ=255)
                                                 aaCAGCATGCCCAATGAACTGATTTTATTTTGAGCTAt                                          <  1:1849883/38‑1 (MQ=255)
                                                 aaCAGCATGCCCAATGAACTGATTTATTTTGAGCTAt                                          <  1:630629/37‑1 (MQ=38)
                                                       aTGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAAAAg     <  1:856530/69‑1 (MQ=255)
                                                         gCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTACCGGATGGCAAAAGCa   <  1:3342196/69‑1 (MQ=255)
                                                          cccAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAAAAGCat  >  1:2101414/1‑69 (MQ=255)
                                                          cccAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAAAAGCat  <  1:2536834/69‑1 (MQ=255)
                                                          cccAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAAAAGCat  <  1:1993649/69‑1 (MQ=255)
                                                          |                                                                    
AGGGCCGAATGAAAAGTGGATAATCCCCATCATTAAAGCGATGGCACCGAACAGCATGCCCAATGAACTGATTTTATTTTGAGCTATGATACCCATCTGACCTCCTTGTAGCGGATGGCAAAAGCAT  >  W3110S.gb/3783907‑3784033

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: