Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,693,744 T→C 8.9% E274E (GAA→GAG uidC ← predicted outer membrane porin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,693,7440TC 8.9% 114.4 / 3.3 45E274E (GAA→GAGuidCpredicted outer membrane porin protein
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base C (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAACCACCACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGCAACTCCCCTTTC  >  W3110S.gb/1693675‑1693813
                                                                     |                                                                     
aTAACCACCACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGttt                                                                      >  1:3155430/1‑71 (MQ=255)
 tAACCACCACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGttt                                                                      >  1:308379/1‑70 (MQ=255)
 tAACCACCACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGttt                                                                      >  1:995153/1‑70 (MQ=255)
 tAACCACCACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGttt                                                                      >  1:4112004/1‑70 (MQ=255)
 tAACCACCACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTc                                                                     <  1:4181300/71‑1 (MQ=255)
   accaccACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCAt                                                                   <  1:3344954/71‑1 (MQ=255)
            gCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGggctggc                                                          >  1:1044932/1‑71 (MQ=255)
               aaaTACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCAt                                                                   >  1:2777570/1‑59 (MQ=255)
               aaaTACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGctggctgg                                                       <  1:1679062/71‑1 (MQ=255)
                 aTACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTCTCATTGGGCTGGCTGGTa                                                     >  1:2318218/1‑71 (MQ=255)
                  tACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTc                                                                     <  1:2969860/54‑1 (MQ=255)
                  tACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTCTCATTGGGCTGGCTGGTAt                                                    <  1:700484/71‑1 (MQ=255)
                   acacTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTg                                                                 >  1:1147287/1‑57 (MQ=255)
                    cacTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCAt                                                                   >  1:4188249/1‑54 (MQ=255)
                     acTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTAc                                                 >  1:434109/1‑71 (MQ=255)
                        tCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGAc                                              >  1:409361/1‑71 (MQ=255)
                         ccccACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACt                                             >  1:2726008/1‑71 (MQ=255)
                         ccccACGGCGCATTCCAGGTAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACt                                             <  1:2027555/70‑1 (MQ=255)
                           ccACGGCGCATTCCAGGTCTATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTa                                           <  1:2634445/71‑1 (MQ=255)
                             aCGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTAt                                                    <  1:120863/60‑1 (MQ=255)
                                 cgcATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCt                                      <  1:4032653/70‑1 (MQ=255)
                                     ttCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGt                                                      <  1:1912574/50‑1 (MQ=255)
                                      tCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTa                                                     >  1:1552322/1‑50 (MQ=255)
                                       ccAGGTCAATTGTCCACTAACCAACGCCGTTTCAt                                                                   >  1:2668436/1‑35 (MQ=255)
                                          ggTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTgcg                            >  1:4010282/1‑71 (MQ=255)
                                             cAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTa                                                     >  1:2008267/1‑43 (MQ=255)
                                             cAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTc                                    >  1:3611514/1‑60 (MQ=255)
                                             cAATTGTCCACTAACCAACGCCGTTTAATTGGGctggctgg                                                       >  1:3060766/1‑41 (MQ=255)
                                             cAATTGTCCACTAACCAACGCCGTCTCATTGGGCTGGCTGGTATTACGACTTAAACCtt                                     >  1:3242969/1‑59 (MQ=255)
                                                ttGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGt                                                      <  1:45168/39‑1 (MQ=255)
                                                ttGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGt                                                      <  1:947561/39‑1 (MQ=255)
                                                ttGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTAtt                                                   <  1:510280/42‑1 (MQ=255)
                                                ttGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTaaa                      <  1:4625850/71‑1 (MQ=255)
                                                 tGTCCACTAACCAACGCCGTTTCATTGGGctggctgg                                                       <  1:88591/37‑1 (MQ=255)
                                                  gTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTaaaa                     <  1:921913/70‑1 (MQ=255)
                                                  gTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAg                    >  1:2052437/1‑71 (MQ=255)
                                                     cACTAACCAACGCCGTTTCATTGGGCTGGCTGTATTACGACTTaa                                          <  1:2174991/45‑1 (MQ=38)
                                                       cTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTaa                                          <  1:1179134/44‑1 (MQ=255)
                                                       cTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTaa                                          <  1:4396884/44‑1 (MQ=255)
                                                        tAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTaa                                          <  1:2609429/43‑1 (MQ=255)
                                                        tAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGc                                <  1:3796619/53‑1 (MQ=255)
                                                          aCCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTaa                       >  1:2707609/1‑60 (MQ=255)
                                                          aCCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGc              <  1:3175342/69‑1 (MQ=255)
                                                          aCCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGc              <  1:118215/69‑1 (MQ=255)
                                                          aCCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGc              <  1:1434296/69‑1 (MQ=255)
                                                          aCCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGCaa            >  1:1349211/1‑71 (MQ=255)
                                                             aaCGCCGTCTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGCAACTc         <  1:3043167/71‑1 (MQ=255)
                                                             aaCGCCCTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCg                 >  1:593640/1‑63 (MQ=255)
                                                              aCGCCGTTTCATTGGGCTGGCTGGTATTACGACTTaa                                          >  1:3570223/1‑37 (MQ=255)
                                                                     ttCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGCAACTCCCCTTTc  >  1:1487997/1‑70 (MQ=255)
                                                                     |                                                                     
ATAACCACCACTGCCAAATACACTTCCCCACGGCGCATTCCAGGTCAATTGTCCACTAACCAACGCCGTTTCATTGGGCTGGCTGGTATTACGACTTAAACCTTCGAGCTGCGCGTAAAAGCCGAGCAACTCCCCTTTC  >  W3110S.gb/1693675‑1693813

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: