Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 553,029 T→G 17.3% A45A (GCA→GCC lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb553,0290TG17.3% 43.4 / 9.9 29A45A (GCA→GCClpxHUDP‑2,3‑diacylglucosamine pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (13/10);  new base G (3/2);  total (17/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.40e-01

ATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAA  >  W3110S.gb/552960‑553089
                                                                     |                                                            
aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg                                                             >  1:743152/1‑71 (MQ=255)
aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg                                                             >  1:2250739/1‑71 (MQ=255)
     acacCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa                                                        <  1:3441602/71‑1 (MQ=255)
      tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                                                        >  1:2131296/2‑70 (MQ=255)
      tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                                                        >  1:380728/2‑70 (MQ=255)
      tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                                                        >  1:700069/2‑70 (MQ=255)
      cacCGCTTTGAGCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa                                                        <  1:1702986/70‑1 (MQ=255)
        ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc                                                     <  1:141476/71‑1 (MQ=255)
        ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc                                                     <  1:1991868/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)
                                                                     |                                                            
ATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATATACAGCGCGTCGGCCTTGCGGGCTTCCCCCGCTAAA  >  W3110S.gb/552960‑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: