Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 553,023 A→G 16.7% I47I (ATT→ATC lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb553,0230AG16.7% 62.1 / 6.7 30I47I (ATT→ATClpxHUDP‑2,3‑diacylglucosamine pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (15/10);  new base G (3/2);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAA  >  W3110S.gb/552955‑553089
                                                                    |                                                                  
ccGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAAt                                                                  >  1:2104609/1‑71 (MQ=255)
     aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg                                                             >  1:743152/1‑71 (MQ=255)
     aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg                                                             >  1:2250739/1‑71 (MQ=255)
          acacCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa                                                        <  1:3441602/71‑1 (MQ=255)
           tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                                                        >  1:700069/2‑70 (MQ=255)
           tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                                                        >  1:380728/2‑70 (MQ=255)
           tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                                                        >  1:2131296/2‑70 (MQ=255)
           cacCGCTTTGAGCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa                                                        <  1:1702986/70‑1 (MQ=255)
             ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc                                                     <  1:1991868/71‑1 (MQ=255)
             ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc                                                     <  1:141476/71‑1 (MQ=255)
             ccGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCaaa                                                      >  1:1629352/1‑70 (MQ=255)
              cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACa                                                    <  1:1679859/71‑1 (MQ=255)
              cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACa                                                    <  1:235084/71‑1 (MQ=255)
               gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCaaa                                                      >  1:1808479/1‑68 (MQ=255)
                   tGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCg                                               >  1:16782/1‑71 (MQ=255)
                    gATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCAATGCTACAAACAGATCGc                                              >  1:3305373/1‑71 (MQ=255)
                            ccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAata                                      <  1:2444161/71‑1 (MQ=255)
                                  tCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACa                                  <  1:1541132/69‑1 (MQ=255)
                                    ttGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAata                                      >  1:3773586/1‑63 (MQ=255)
                                    ttGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAgcgc                              >  1:168845/1‑71 (MQ=255)
                                      gCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGt                            >  1:2878433/1‑71 (MQ=255)
                                          tGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAAtata                                    >  1:2194525/1‑59 (MQ=255)
                                               tGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTg                    <  1:1197391/70‑1 (MQ=255)
                                                    tGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGcgcg                             >  1:1967935/1‑56 (MQ=255)
                                                           gtcgtcgCCAATCCAAGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTg                    >  1:1021664/1‑58 (MQ=255)
                                                            tcgtcgCCAATCCATGCTTCAAACAGATCGCCAAGaa                                        <  1:1850897/37‑1 (MQ=255)
                                                            tcgtcgCCAATCCATGCTTCAAACAGAGCGCCAATAATATACAGCGCGTCGGCCTTg                    <  1:3431283/57‑1 (MQ=255)
                                                              gtcgCCAATCCATGCTTCAAACAGATCGCCAAGAATATACa                                  >  1:871631/1‑41 (MQ=255)
                                                                 gCCAATCCATGCTTCAAACAGATCGCCAAGAATATAc                                   >  1:362575/1‑37 (MQ=255)
                                                                 gCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTaaa  <  1:2204911/70‑1 (MQ=255)
                                                                    |                                                                  
CCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAA  >  W3110S.gb/552955‑553089

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: