Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,594,973 A→G 11.1% pseudogene (807/1401 nt) ydeU ← ECK1502:JW1502:b1509; conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,594,9730AG11.1% 136.8 / 9.7 54pseudogene (807/1401 nt)ydeUECK1502:JW1502:b1509; conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/20);  new base G (2/4);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GTTATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCACACGATTAC  >  W3110S.gb/1594907‑1595043
                                                                  |                                                                      
gTTATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATAc                                                                     >  1:292930/1‑70 (MQ=255)
  tATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAg                                                                  >  1:2445010/1‑71 (MQ=255)
   aTCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATAcc                                                                    <  1:1466823/68‑1 (MQ=255)
    tCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGt                                                                >  1:468623/1‑71 (MQ=255)
      aaaCCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATAcc                                                                    <  1:865676/65‑1 (MQ=255)
          cAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTa                                                          <  1:3177782/71‑1 (MQ=255)
          cAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCCTa                                                          <  1:535572/71‑1 (MQ=255)
            gTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTa                                                          >  1:621926/1‑69 (MQ=255)
               ataCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGa                                                                        >  1:3081558/1‑52 (MQ=255)
               ataCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                     >  1:1307695/1‑71 (MQ=255)
                taCAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                     <  1:3973215/70‑1 (MQ=255)
                  cAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGa                                                                        >  1:3561712/1‑49 (MQ=255)
                  cAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACa                                                        >  1:1236069/1‑65 (MQ=255)
                     cgcCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCAGGTcgc                                                             >  1:584062/1‑57 (MQ=255)
                       cccAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCt                                              <  1:2179963/70‑1 (MQ=255)
                        ccAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGt                                            <  1:4619734/71‑1 (MQ=255)
                         cAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGt                                            >  1:3741057/1‑70 (MQ=255)
                         cAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTa                                           <  1:1112379/71‑1 (MQ=255)
                          aGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTaa                                          <  1:3886773/71‑1 (MQ=255)
                            cTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGa                                                                        >  1:2487338/1‑39 (MQ=255)
                            cTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCAGGTcgcg                                                            <  1:3677348/51‑1 (MQ=255)
                            cTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAAcc                                        <  1:3233290/71‑1 (MQ=255)
                               tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                     <  1:376576/71‑1 (MQ=255)
                               tCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                     <  1:3516886/71‑1 (MQ=255)
                                cAACATAAGCGCCGGTCTTATTCGCATCGTTCTGa                                                                        >  1:2764939/1‑35 (MQ=255)
                                 aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGGTACCa                                                                   >  1:2885422/1‑39 (MQ=255)
                                 aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGc                                      >  1:1002958/1‑68 (MQ=255)
                                 aaCATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                     <  1:928388/69‑1 (MQ=255)
                                    aTAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATgc                                >  1:3648827/1‑71 (MQ=255)
                                     tAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATgcg                               >  1:3527397/1‑71 (MQ=255)
                                       aGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGcc                             >  1:4378512/1‑71 (MQ=255)
                                       aGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGcc                             >  1:1723045/1‑71 (MQ=255)
                                       aGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCAGCGCTGTAACCGCTGATGCGcc                             >  1:2992317/1‑70 (MQ=255)
                                        gcgcCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                     >  1:469277/1‑62 (MQ=255)
                                         cgcCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGa                                   >  1:4637621/1‑63 (MQ=255)
                                          gcCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                     >  1:2847759/1‑44 (MQ=255)
                                          gcCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGcc                                                     >  1:3683579/1‑44 (MQ=255)
                                          gcCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGt                                            <  1:4460369/53‑1 (MQ=255)
                                          gcCGGTCTTATTCGCATCGTTCTGATACCACGTcg                                                              >  1:4465699/1‑35 (MQ=255)
                                            cGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATcc                        <  1:2327962/71‑1 (MQ=255)
                                              gTCTTATTCGCATCGTTCTGGTACCAGGTCGCGTACAGCCCAgcg                                                <  1:2500458/45‑1 (MQ=255)
                                               tCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                     <  1:2376435/55‑1 (MQ=255)
                                                cTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGAtt                    >  1:3220947/1‑71 (MQ=255)
                                                 ttATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAgcg                                                <  1:3612605/42‑1 (MQ=255)
                                                  tATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTta                  >  1:932661/1‑71 (MQ=255)
                                                  tATTCGCATCGTTCTGATACCAGGTCGCGTACAGCACAGCGCTGTAAACGCTGATGCGCCCATCCGATTta                  >  1:3115163/1‑71 (MQ=255)
                                                    ttCGCATCGTTCTGGTACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTtata                <  1:3476873/71‑1 (MQ=255)
                                                    ttCGCATCGTTCTGGTACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTtata                <  1:469957/71‑1 (MQ=255)
                                                    ttCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAg                                                  >  1:3672014/1‑37 (MQ=255)
                                                    ttCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAg                                                  >  1:3204217/1‑37 (MQ=255)
                                                    ttCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTtata                <  1:3239930/71‑1 (MQ=255)
                                                    ttCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTtat                 >  1:2611081/1‑70 (MQ=255)
                                                      cGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTa                                           >  1:2997370/1‑42 (MQ=255)
                                                          tCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCca           >  1:1785296/1‑70 (MQ=255)
                                                            gTTCTGATACCAGGTCGCGTATAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCacac        <  1:720592/71‑1 (MQ=255)
                                                            gTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCacac        <  1:1609266/71‑1 (MQ=255)
                                                               cTGATACCAGGTCGCGTACAGGCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCACACGAt     >  1:1766845/1‑71 (MQ=255)
                                                               cTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCt                                     >  1:1549409/1‑39 (MQ=255)
                                                                  aTACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCACACGATTAc  <  1:4434336/71‑1 (MQ=255)
                                                                  |                                                                      
GTTATCAAACCAGTTATACAGCGCCCAGCTGTCAACATAAGCGCCGGTCTTATTCGCATCGTTCTGATACCAGGTCGCGTACAGCCCAGCGCTGTAACCGCTGATGCGCCCATCCGATTTATAACCCACACGATTAC  >  W3110S.gb/1594907‑1595043

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: