Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 1,846,641 (A)8→7 intergenic (‑56/+72) ynjI ← / ← topB predicted inner membrane protein/DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,846,6340A.91.8% 225.6 / 16.0 61intergenic (‑49/+79)ynjI/topBpredicted inner membrane protein/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base A (5/0);  new base . (20/36);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 8.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

ACATTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAG  >  W3110S.gb/1846587‑1846694
                                               |                                                            
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTc                                       <  1:3590215/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTc                                       <  1:1035379/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTc                                       <  1:3211373/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTc                                       <  1:17725/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTc                                       <  1:1880914/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTc                                       <  1:1386592/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTc                                       <  1:4056874/70‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAAtt                                          <  1:306341/66‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAAtt                                          <  1:1876109/66‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAAtt                                          <  1:273943/66‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAAtt                                          <  1:3226897/66‑1 (MQ=255)
   ttCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAAtttt                                        >  1:45249/1‑66 (MQ=255)
   ttCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAAtttt                                        >  1:1649387/1‑66 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:1444205/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:3414035/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:2498872/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:690506/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:3662852/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCAATGCAAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:1058048/1‑71 (MQ=255)
            attaCATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAAATATCCAACGAATTTTCTTGATCCGCTg                            >  1:1815997/1‑70 (MQ=255)
            attaCATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAAATATCCAACGAATTTTCTTGATCCGCTg                            >  1:3853145/1‑70 (MQ=255)
            attaCATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAAATATCCAACGAATTTTCTTGATCCGCTg                            >  1:1910043/1‑70 (MQ=255)
            attaCATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAAATATCCAACGAATTTTCTTGATCCGCTg                            >  1:308919/1‑70 (MQ=255)
            attaCATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAAATATCCAACGAATTTTCTTGATCCGCTg                            >  1:781747/1‑70 (MQ=255)
               aCATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGatat                       >  1:3047274/1‑71 (MQ=255)
                       tAGAGGCGATAATTTTCTCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:55857/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:377783/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3897391/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:2281505/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3457740/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3535844/70‑1 (MQ=255)
                         gagGCTATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:2856854/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:3921245/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:500849/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:2221321/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:992796/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3089790/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:140205/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1828827/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1417644/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1310178/68‑1 (MQ=255)
                             cGATAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1237583/64‑1 (MQ=255)
                               aTAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:2726719/71‑1 (MQ=255)
                               aTAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:4119888/71‑1 (MQ=255)
                               aTAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:759893/71‑1 (MQ=255)
                               aTAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:1102902/71‑1 (MQ=255)
                               aTAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:3925603/71‑1 (MQ=255)
                                tAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:2756054/70‑1 (MQ=255)
                                tAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:19840/70‑1 (MQ=255)
                                tAATTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:4050674/70‑1 (MQ=255)
                                 aaTTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATAtt                      >  1:3783812/1‑54 (MQ=255)
                                 aaTTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATAtt                      >  1:1716630/1‑54 (MQ=255)
                                 aaTTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGaa     >  1:1650414/1‑71 (MQ=255)
                                 aaTTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGaa     >  1:1090391/1‑71 (MQ=255)
                                 aaTTTTATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGaa     >  1:1085698/1‑71 (MQ=255)
                                    tttATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:3069518/1‑71 (MQ=255)
                                    tttATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:1681554/1‑71 (MQ=255)
                                    tttATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:1775592/1‑71 (MQ=255)
                                    tttATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:1844518/1‑71 (MQ=255)
                                    tttATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:2086411/1‑71 (MQ=255)
                                    tttATCCATGCAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:2914518/1‑71 (MQ=255)
                                               |                                                            
ACATTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGCAAAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAG  >  W3110S.gb/1846587‑1846694

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: