Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 929,528 G→C 14.3% A59G (GCG→GGG)  U712_04525 ← putative lipoprotein ygaO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881929,5280GC14.3% 98.4 / 2.7 52A59G (GCG→GGG) U712_04525putative lipoprotein ygaO
Reads supporting (aligned to +/- strand):  ref base G (12/32);  new base C (6/2);  total (18/34)
Fisher's exact test for biased strand distribution p-value = 1.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATATATGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAGGCCAAA  >  CP006881/929400‑929664
                                                                                                                                |                                                                                                                                        
cGATATATGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTATTTAAGTAAAAAAAAACCGCACAcgccgc                                                                                                                                >  6:121737/1‑139 (MQ=255)
 gATATATGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCATTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGca                                                                                                                               <  3:362943/139‑1 (MQ=255)
      atGCCGCGCGACAGGCTTTTGCCGTCATCATCTGCTCTCTTTCACCCTGCTCTTCAAACACATCCGGCAGAAGTCGCGGCTTTCGGGCAGGGCGGCTTTTTATCTGTAAGTACAAATAACCCCCAGACGCCGCaacaaa                                                                                                                          <  6:120729/139‑1 (MQ=255)
       tGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCCGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTACGTACAAATAAACCCCCGACGCCGCaacaaa                                                                                                                          <  1:279883/138‑1 (MQ=255)
       tGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaa                                                                                                                          <  5:121737/138‑1 (MQ=255)
       tGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTATGTAAGAACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                         <  1:215550/139‑1 (MQ=255)
        gACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGGGGGTTTTTTTCTGTAAGTAAAAATAAACCGCAGACGCCGCaacaaac                                                                                                                         <  1:303182/138‑1 (MQ=255)
        gACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGAAAGTACAAATAAACCGCCGACGCCGCaacaaac                                                                                                                         <  3:281483/138‑1 (MQ=255)
        gACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTAAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                         <  3:192908/138‑1 (MQ=255)
          cgcgcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTGTTTTTTTGTAATATAAAAAAAACCCCACACACCGCcaaaaaaac                                                                                                                       >  8:346492/1‑135 (MQ=255)
          cgcgcgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGTTTTTTTTGTTAGTTCAAATAAACCCCCGCCGCCGCCACAAACAg                                                                                                                       >  7:321248/1‑138 (MQ=255)
              cgACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAAGTCCTCTGGCAGAAGTAGCTGCTTTCTGGCTGGGCGGATTTTTTTCTGTAAGTACAAATAACCCGCAGACGCCGCAACAAACAGAATaa                                                                                                                  <  7:346492/139‑1 (MQ=255)
                     cTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGt                                                                                                           <  5:88198/139‑1 (MQ=255)
                        ttGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTg                                                                                                          <  7:121169/137‑1 (MQ=255)
                            cggcttcatcTGCTCTGTTTCATCCTGCTCTTCAAACTCATCAGGCATAAGTCGCTGCTTTCTGGCAGGGCGGATTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGttt                                                                                                    <  4:202491/134‑1 (MQ=255)
                            cGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGttt                                                                                                    <  8:60665/139‑1 (MQ=255)
                            cGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGAACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGttt                                                                                                    <  8:139689/139‑1 (MQ=255)
                            cGTCATCATCTGCTCTCTTTAATCCTGCTCTTCAAACTCATAAGGCAGAATTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAAAAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGttt                                                                                                    <  3:143543/139‑1 (MQ=255)
                               catcatCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGTTTTTTGTGTAATTACAAATAAACCGCAGCCGCCGCAACAAACAGAATAAAGGGAGTGCGGGTTtaa                                                                                                 >  5:287100/1‑139 (MQ=255)
                                atcatcTGCTCTCTTTCATCCTGATCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTTATTAATTACAAAAAAACCGCACACCCCGCAAACAACACAATATAGGGACTTAGGGTTaaa                                                                                                 >  8:89889/1‑138 (MQ=255)
                                 tcatcTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATc                                                                                               <  7:362027/139‑1 (MQ=255)
                                     cTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGGTTTTTATGTTAGTACAAATAAACCCCAGACGCCGCAACAAACAGAATAAAAGGAGTGAGGGTTAAATCaaa                                                                                            >  3:285491/1‑138 (MQ=255)
                                        ctctctTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTTAGTACAAATAAACCGCAGCCGCCGCAACAAACAGCATAAAGGGAGTGAGGTGTTAATCAAATAc                                                                                         >  8:50000/1‑138 (MQ=255)
                                        ctctctTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGGTTTTTTTTTTAAAAAAAATAAAAACCAGCCGCCGCCACAAACACAAAAAAAGAGTTGAGGTTTTAATAAAAAAc                                                                                         >  5:210095/1‑138 (MQ=255)
                                           tctTTCATCCTGCTCTTCAAACTCATCAGGCCGAAGTCGCTGCTTTCGGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCa                                                                                       <  8:321248/137‑1 (MQ=255)
                                             tttCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGca                                                                                                                               >  8:13093/1‑95 (MQ=255)
                                             tttCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGca                                                                                                                               >  8:18414/1‑95 (MQ=255)
                                             tttCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGca                                                                                                                               <  7:18414/95‑1 (MQ=255)
                                             tttCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGca                                                                                                                               <  7:13093/95‑1 (MQ=255)
                                             tttCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTATGTAAGTTCAAATAAACCCCAGACGCCGCCACAAACACAATATAGGGAGTTGGGGTTTAATCAAAAACCATccc                                                                                   >  6:89692/1‑139 (MQ=255)
                                             tttCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCGATTTTTAAGTTAATTAAAATAAAACCCACACCCCCCCACCACAAAAATAAAGGGATTGACGGTTAAAAAAAATAACAAccc                                                                                   >  6:121548/1‑139 (MQ=255)
                                                      gCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCtt                                                                          >  3:347766/1‑139 (MQ=255)
                                                       ctctTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTGTGTAAGTTCAAATAAACCGCAGACGCCGCAACAAACAGAATAAATGGAGTGAGGTTTTAATAAAATACCATCCCGCCGCGCttt                                                                         >  4:174756/1‑139 (MQ=255)
                                                       ctctTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTATTTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCGCttt                                                                         >  5:287952/1‑139 (MQ=255)
                                                                  aatcaGGCATAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGaaaaa                                                              <  3:261653/138‑1 (MQ=255)
                                                                   atcaGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGaaaaa                                                              <  4:347766/138‑1 (MQ=255)
                                                                    tcaGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                         <  8:242105/78‑1 (MQ=255)
                                                                    tcaGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCaacaaac                                                                                                                         >  7:242105/1‑78 (MQ=255)
                                                                             aGTAGCTGCTTTCTGGCAGGGCGGGTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATg                                                    <  6:287952/138‑1 (MQ=255)
                                                                                            gctgggCGGATTTTTTTCTGTCAGTAAAAATAACCCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCa                                     <  5:260921/135‑1 (MQ=255)
                                                                                            gccggggggCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCGGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCa                                     <  7:89889/135‑1 (MQ=255)
                                                                                            gCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCaa                                    <  1:19384/139‑1 (MQ=255)
                                                                                                   ggcTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGa                             <  7:50000/139‑1 (MQ=255)
                                                                                                                tAAGTACAAATAACCCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCa                 <  4:227521/138‑1 (MQ=255)
                                                                                                                     aCAAAAAAACCGACGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATg            <  6:3132/138‑1 (MQ=255)
                                                                                                                      cAAATAAACCGCAGACGCCGCAACAAACAGAATACAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCa          <  2:230447/139‑1 (MQ=255)
                                                                                                                       aaaTAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCaa         <  5:160186/139‑1 (MQ=255)
                                                                                                                        aaTAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGACAAAAGATGCTATGCCTTCTGAGGCTGAAAAGTATGAACTCAGCGCCCAGCATGCAAg        >  6:202352/1‑139 (MQ=255)
                                                                                                                         aTAACCCGCCGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAgg       <  3:66094/139‑1 (MQ=255)
                                                                                                                         aTAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAgg       <  3:174756/139‑1 (MQ=255)
                                                                                                                         aTAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAgg       <  3:174859/139‑1 (MQ=255)
                                                                                                                              ccGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAGGCCaaa  >  5:118511/1‑139 (MQ=255)
                                                                                                                                |                                                                                                                                        
CGATATATGACGCGCGACATGCTTTTGCCGTCATCATCTGCTCTCTTTCATCCTGCTCTTCAAACTCATCAGGCAGAAGTAGCTGCTTTCTGGCAGGGCGGCTTTTTTTCTGTAAGTACAAATAAACCGCAGACGCCGCAACAAACAGAATAAAGGGAGTGAGGTTTAAATCAAATACCATCCCGCCGCTCTTTGTGCCGAAAAAAGATGCTATGCCTTCTGAGGCTTCAAAGTATGAAATCAGCGCCCAGCATGCAAGGCCAAA  >  CP006881/929400‑929664

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