Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 644,304 C→T intergenic (‑78/+36) rna ← / ← citT ribonuclease I/citrate:succinate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb644,3040CT96.7% 80.9 / ‑3.7 30intergenic (‑78/+36)rna/citTribonuclease I/citrate:succinate antiporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (2/27);  total (3/27)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCA  >  W3110S.gb/644237‑644370
                                                                   |                                                                  
gAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCtggtgg                                                                  >  1:2165695/1‑70 (MQ=255)
   cTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGt                                                              <  1:908569/71‑1 (MQ=255)
   cTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGt                                                              <  1:687770/71‑1 (MQ=255)
   cTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGt                                                              <  1:2210274/71‑1 (MQ=255)
   cTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGt                                                              <  1:1767286/71‑1 (MQ=255)
           tacaCTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGAttt                                                      <  1:1633080/71‑1 (MQ=255)
              acTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGt                                                              <  1:1360746/60‑1 (MQ=255)
                   aTACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTa                                                     <  1:1071282/64‑1 (MQ=255)
                     aCTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTAtt                                                   >  1:774497/1‑64 (MQ=255)
                      ctatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGttt                                           >  1:4018346/1‑71 (MQ=255)
                       tatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGtg                                                             <  1:699961/52‑1 (MQ=255)
                       tatctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGtg                                                             <  1:1867443/52‑1 (MQ=255)
                         tctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:1246465/69‑1 (MQ=255)
                         tctatTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:139158/69‑1 (MQ=255)
                          ctatAAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTa                                                     <  1:1569631/57‑1 (MQ=39)
                               aaaTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGtga                                                            >  1:1035599/1‑45 (MQ=255)
                                 aTCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:699867/61‑1 (MQ=255)
                                  tCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:1628246/60‑1 (MQ=255)
                                    atatAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:938550/58‑1 (MQ=39)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:171352/58‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:2826077/58‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:3416193/58‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:1529411/58‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:1100932/58‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTg                                          <  1:1035594/58‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTa                              <  1:2519715/70‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTa                              <  1:2982625/70‑1 (MQ=255)
                                    aTAAAAAAGCCCGCCAGGTGGGCTTACCTGGTGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTa                              <  1:1344005/70‑1 (MQ=255)
                                                              ccTGGTGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCca                        <  1:187536/50‑1 (MQ=255)
                                                               cTGGTGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGTGAGAATCGgccagcca  <  1:621051/71‑1 (MQ=255)
                                                               cTGGTGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGgccagcca  <  1:1290745/71‑1 (MQ=255)
                                                                   |                                                                  
GAACTCATACATACACTGAATACTATCTATTAAATCATAAAAAAGCCCGCCAGGTGGGCTTACCTGGCGGGCGTGATGATTTATTCAGCGTTTGGCGAACGTATTAGTTCCACATGGCGAGAATCGGCCAGCCA  >  W3110S.gb/644237‑644370

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: