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
                                               |                                                            
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTc                                       <  1:3590215/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTc                                       <  1:1035379/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTc                                       <  1:3211373/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTc                                       <  1:17725/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTc                                       <  1:1880914/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTc                                       <  1:1386592/70‑1 (MQ=255)
acaTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTc                                       <  1:4056874/70‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAAtt                                          <  1:306341/66‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAAtt                                          <  1:1876109/66‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAAtt                                          <  1:273943/66‑1 (MQ=255)
 caTTCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAAtt                                          <  1:3226897/66‑1 (MQ=255)
   ttCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAAtttt                                        >  1:45249/1‑66 (MQ=255)
   ttCCCTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAAtttt                                        >  1:1649387/1‑66 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:1444205/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:3414035/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:2498872/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:690506/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATc                                 >  1:3662852/1‑71 (MQ=255)
     cccTGTTATTACATACTTTAGAGGCGATAATTTTATCAATGC‑AAAAAAATATCCAACGAATTTTCTTGATc                                 >  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)
               aCATACTTTAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGatat                       >  1:3047274/1‑71 (MQ=255)
                       tAGAGGCGATAATTTTCTCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:55857/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:377783/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3897391/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:2281505/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3457740/70‑1 (MQ=255)
                       tAGAGGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3535844/70‑1 (MQ=255)
                         gagGCTATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:2856854/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:3921245/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:500849/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGc                              <  1:2221321/54‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:992796/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:3089790/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:140205/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1828827/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1417644/68‑1 (MQ=255)
                         gagGCGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1310178/68‑1 (MQ=255)
                             cGATAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGaa                <  1:1237583/64‑1 (MQ=255)
                               aTAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:2726719/71‑1 (MQ=255)
                               aTAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:4119888/71‑1 (MQ=255)
                               aTAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:759893/71‑1 (MQ=255)
                               aTAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:1102902/71‑1 (MQ=255)
                               aTAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:3925603/71‑1 (MQ=255)
                                tAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:2756054/70‑1 (MQ=255)
                                tAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:19840/70‑1 (MQ=255)
                                tAATTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCg       <  1:4050674/70‑1 (MQ=255)
                                 aaTTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATAtt                      >  1:3783812/1‑54 (MQ=255)
                                 aaTTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATAtt                      >  1:1716630/1‑54 (MQ=255)
                                 aaTTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGaa     >  1:1650414/1‑71 (MQ=255)
                                 aaTTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGaa     >  1:1090391/1‑71 (MQ=255)
                                 aaTTTTATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGaa     >  1:1085698/1‑71 (MQ=255)
                                    tttATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:3069518/1‑71 (MQ=255)
                                    tttATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:1681554/1‑71 (MQ=255)
                                    tttATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:1775592/1‑71 (MQ=255)
                                    tttATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:1844518/1‑71 (MQ=255)
                                    tttATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  1:2086411/1‑71 (MQ=255)
                                    tttATCCATGC‑AAAAAAATATCCAACGAATTTTCTTGATCCGCTGGATATTTCAGAATATGACTCGAATAg  >  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: ‑