Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,026,635 T→G 20.7% I41S (ATC→AGC)  U712_05035 → Uncharacterized protein yhaP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,026,6350TG20.7% 50.6 / 2.1 41I41S (ATC→AGC) U712_05035Uncharacterized protein yhaP
Reads supporting (aligned to +/- strand):  ref base T (18/14);  new base G (7/2);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGGGGCCGTCTATGAATAAATTTTGGATCATGCTTTCTCATACATACAAAAATAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCTTGATGTGAAGCTGTCGAAG  >  CP006881/1026501‑1026771
                                                                                                                                      |                                                                                                                                        
gTGGGGGCCGTCTATGAATAAATTTTGGATCATGCTTTCTCATACATACAAAAATTAGATTATGGCAAAATCCTTTATTATTTCGACTGTGATCACGGTGCGGCTTGTTTTAGTTGTAACGAATCTGGCGTCGATccct                                                                                                                                      >  8:257757/1‑136 (MQ=255)
  gggggCCGTCTATGAATAAATTTTGGATCATGCTTTCTCATACATACAAAAATAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCAttta                                                                                                                                    >  7:352156/1‑138 (MQ=255)
        cGTCTATGAATAAATTTTGGATCATGCGTTCACATACATACAAAAATAAGATTATGGCAAAATCCTTTTTTATTTCGCCTGTAATCACGGTGCTGCGTGTTTTAGTTGTAACGAATCTGGAGTCGATGCTTTCtttttt                                                                                                                              >  8:83477/1‑139 (MQ=255)
                                          tacatacaAAAATAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGAGGTTTTCTTTGTTTGTAGGTGGCGATGCGAAATAAAAAATAGCCgtt                                                                                            >  2:152010/1‑139 (MQ=255)
                                            catacaAAAATAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGAGGAAATATTTGTTTGAAGTGGACCATAGGAAAGAAAAAATAGCCggtgt                                                                                          >  1:90122/1‑139 (MQ=255)
                                               acaAAAATAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTTGTTCAAGGTGACGATGAGAAATAAAAAATAGCCGGTGTtga                                                                                       >  5:343597/1‑139 (MQ=255)
                                                 aaaaaTAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATGAATTATTTTTTTGAAAGGGGAGAAGCGAAAGAAAAAAAAGCCGGTTTtgatg                                                                                     >  6:208179/1‑139 (MQ=255)
                                                 aaaaaTAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGGGAAAGAAAAAAACGCCGTTGTtgat                                                                                      >  1:259565/1‑138 (MQ=255)
                                                 aaaaaTAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGAGCATTTATTTTTTTCAAGGTGACGATGCGAAAGAAAAAAACGCCGTTGTtgatt                                                                                     >  1:251119/1‑138 (MQ=255)
                                                   aaaTAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTGTTTGTTTGAAGGTGACGATGCGAAATAAAAAATCGCCgttgtttgtga                                                                                    >  2:66175/1‑133 (MQ=255)
                                                   aaaTAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGAGGATTTATTTTTTTTAATGTGAAGATGCGAAAGAAAAAATCGCCGTTGGtgttga                                                                                    >  2:346748/1‑138 (MQ=255)
                                                         gATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGaaacaaac                                                                             <  2:208120/139‑1 (MQ=255)
                                                              tGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAAtt                                                                        <  1:346748/139‑1 (MQ=255)
                                                              tGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAAtt                                                                        <  1:66175/139‑1 (MQ=255)
                                                                    aaTCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTGAAGGTGACGATGGGAAAGAAAAAATAGCCCTTGTTGATGAAACAAACGAATTatatac                                                                  >  3:63284/1‑137 (MQ=255)
                                                                                tttCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATGATTTGTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATAGCCGTTGTTGATGAAACAAACGAATTATATGCGGTTTTTTCt                                                        >  2:25599/1‑137 (MQ=255)
                                                                                  tCGTCTGGCATCACGGTGCTGCCTGTTTTAGTTGTAACGAACCTGGAGTCGAGCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAACTCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAag                                                     <  7:316269/138‑1 (MQ=255)
                                                                                  tCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAag                                                     <  4:354246/138‑1 (MQ=255)
                                                                                  tCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAag                                                     <  7:83477/138‑1 (MQ=255)
                                                                                  tCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGAGCAATTCTATTTTTCAAGGGGACGCTGCGAAAGAAAAAATAGCCCTTGTTGATGCAAAAAACCAAATATATACCTTTTTTtata                                                       >  5:77039/1‑136 (MQ=255)
                                                                                  tCGACTGTCATCACGGGGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAag                                                     <  5:208179/138‑1 (MQ=255)
                                                                                        gTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTATTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTATAATCAGCTg                                               >  8:255662/1‑138 (MQ=255)
                                                                                        gTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGAGGAATTGTTTGTTTACAGGTTACCCTGCGAAAGAAAAAATAGCCCGTTTTGATGTAACAAAACAATTATTTACGGTTTTTTATAagccgcag                                               >  2:302382/1‑136 (MQ=255)
                                                                                           atcaTGGTGTTGTTTGTTTTTTTTTTAACGATTGTGTTTTGTATATTTTATTTTTTTAAATGTGACGAGGGGAATTAAAAAATCGCCTTTGTTGATGAAACAAACTAATTATATCCGGTTTTTTCTAAGCAGCTGaaa                                            <  2:213914/138‑1 (MQ=255)
                                                                                           atcaGGGTGCTGCTTGTTTTTGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATTCGAAATAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGaaa                                            <  6:118203/138‑1 (MQ=255)
                                                                                           atcaCGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGaaa                                            <  7:155024/138‑1 (MQ=255)
                                                                                            gcaCGGGGGTGCTTGGTTTCGTTGTCACGACTCTGGCGGCGAGCATTTCTTTGTTTCAAGGTGACGAGGCGAAACAAAAACTCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGc                                          <  6:23478/138‑1 (MQ=255)
                                                                                             caCGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCg                                         <  2:338927/139‑1 (MQ=255)
                                                                                                   gctgctTGTTTTAGTTGTAACGAATCTGGAGTCGATGATTTATTTGTTTAAATGTGACGATGCGAAATAAAAAATCCCCCGTGTTGATGAAAAAAAACAAATATATACGGTTTTTTTTAAGCAGCTGAAAACGGTGGa                                    >  5:118203/1‑138 (MQ=255)
                                                                                                     tgctTGTTTTAGTTGTAACGAATCTGGAGTCGAGGATTTTGTTTTTTAAAGGTGAGGCTGGGAAAAAAAAAAAAGCCCGTGTTGTTTTTAAAAAACATATATATTCGGTTTTTTGTTAGACACAGAAAGAGGTGGATa                                  >  2:345172/1‑138 (MQ=255)
                                                                                                            tttAGTTGTAACGAATCTGGAGTCGATCATTTATTTGTTTCAAGGTGACGATGCGAAAGAAAAAATAGCCGTTGTTGATGAAACAAAAGAATTATATACGGTTTTTTCTAAGCAGCAGAAAGCGGTGGATACAGACGGc                          >  6:352442/1‑139 (MQ=255)
                                                                                                            tttAGTTGAAACGAATCTGGAGTCGATAATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACgg                           <  4:63284/138‑1 (MQ=255)
                                                                                                                 ttGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAATAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCt                     >  3:93399/1‑139 (MQ=255)
                                                                                                                 ttGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCt                     <  6:91739/139‑1 (MQ=255)
                                                                                                                    tAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGGCGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAAAAAACGAAATTTATACGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATATTGa                  >  1:51734/1‑139 (MQ=255)
                                                                                                                      aCGAATCTGGAGTCGATCATTTATTTGTTTGAAGGTGCCGATGCGAAAGAAAAAATACCCGTTGTTGATGAAACAAACCAATTATTTACCGTTTTTTCTAAGCACCTGAAAGCCGTGGATACAGACGGCGATATTGAtt                >  6:224483/1‑138 (MQ=255)
                                                                                                                                aGTCTATCATTTATTTGTTTCATGGTGAGCATGCTAACCAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCTTGATGTGAAGCTGTc      <  6:30196/139‑1 (MQ=255)
                                                                                                                                    gATCATTTCTTTGTTTCAAGGTGACGATGCGAACTAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCTTGATGTGAAGCTGTCGAAg  <  4:258180/139‑1 (MQ=255)
                                                                                                                                    gATCATTTCTTTGTTTCAAGGTGACGATGCGAAATAAAAAATAGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCTTGATGTGAAGCTGTCGAAg  >  7:154087/1‑139 (MQ=255)
                                                                                                                                    gATCATTTATTTGTTTCAAGGTGACGATTCGAAACAAAAAATCGCCGTTGTTGATGAAACAAAAGAATTATATACGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCTTGATGTGAAGCTGTCGAAg  >  3:120933/1‑139 (MQ=255)
                                                                                                                                     aTCATTTATTTGTTTCAAAGTGACCATGCGAAATAAAAAATCCCCGTTGTTGATGAAACAAACGAATTATATACGGTTTTTTCTAATCAGCTGAAAGCGGTTGATAAAGACGGCGAACTTGATGTGAAGATGTCGaat  >  8:258925/1‑137 (MQ=255)
                                                                                                                                      |                                                                                                                                        
GTGGGGGCCGTCTATGAATAAATTTTGGATCATGCTTTCTCATACATACAAAAATAAGATTATGGCAAAATCCTTTATTATTTCGACTGTCATCACGGTGCTGCTTGTTTTAGTTGTAACGAATCTGGAGTCGATCATTTCTTTGTTTCAAGGTGACGATGCGAAAGAAAAAATCGCCGTTGTTGATGAAACAAACGAATTATATCCGGTTTTTTCTAAGCAGCTGAAAGCGGTGGATACAGACGGCGATCTTGATGTGAAGCTGTCGAAG  >  CP006881/1026501‑1026771

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