Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 790,778 G→A intergenic (+277/‑197) kipA → / → ltaS possible allophanate hydrolase subunit 2/possible phosphatidylglycerol‑‑membrane‑oligosaccharide glycerophosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730790,7780GA96.5% 87.3 / ‑3.6 29intergenic (+277/‑197)kipA/ltaSpossible allophanate hydrolase subunit 2/possible phosphatidylglycerol‑‑membrane‑oligosaccharide glycerophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (15/13);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

TTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTAC  >  CP000730/790645‑790904
                                                                                                                                     |                                                                                                                              
ttGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCa                                                                                                                         <  2:74037/141‑1 (MQ=255)
  gAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCaaa                                                                                                                       <  2:210095/141‑1 (MQ=255)
       ttGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACt                                                                                                                  >  2:357441/1‑141 (MQ=255)
        tGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACtt                                                                                                                 <  1:239816/141‑1 (MQ=255)
                aTTACTCTATATTTATTAATTGATTTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATc                                                                                                         <  2:370476/141‑1 (MQ=255)
                aTTACTCTATATTTAGTAATTGATCTTTGTTATTTATTACAATTTGATTACAAAATATAAATGTGTTAATTGGAAGATAATATTAAATAAAGAAACTTACACAAGGAAATATGAGTTATAGCCCAAATTACTTGTTCAATc                                                                                                         <  2:75923/141‑1 (MQ=255)
                    ctctATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAg                                                                                                     <  1:288486/141‑1 (MQ=255)
                                taATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACa                                                                                         >  1:158002/1‑141 (MQ=255)
                                   ttGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAAt                                                                                      <  2:264991/141‑1 (MQ=255)
                                      aGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAAc                                                                                   >  1:59523/1‑141 (MQ=255)
                                        cTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGa                                                                                 >  2:61812/1‑141 (MQ=255)
                                        cTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGa                                                                                 >  2:282789/1‑141 (MQ=255)
                                         tATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGaa                                                                                >  1:36335/1‑141 (MQ=255)
                                                    ttaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACt                                                                     >  2:344815/1‑141 (MQ=255)
                                                       cAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGa                                                                  >  2:215558/1‑141 (MQ=255)
                                                         aTTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAAt                                                                >  1:76137/1‑141 (MQ=255)
                                                                          aaaTTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACtat                                               >  1:135700/1‑141 (MQ=255)
                                                                            aTTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATaa                                             >  2:315830/1‑141 (MQ=255)
                                                                             tttGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAAt                                            <  2:339983/141‑1 (MQ=255)
                                                                                    aTTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAACTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACa                                     <  2:280286/141‑1 (MQ=255)
                                                                                        aaTGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGt                                 <  1:171922/141‑1 (MQ=255)
                                                                                          tGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTa                               <  2:46185/141‑1 (MQ=255)
                                                                                            ataataTTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTGCTATAATTATGACAATGTTAAc                             >  1:309430/1‑141 (MQ=255)
                                                                                               atatTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACttt                          <  2:66975/141‑1 (MQ=255)
                                                                                                    aaataaaGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAAc                     <  2:181080/141‑1 (MQ=255)
                                                                                                         aaGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACt                >  2:223378/1‑141 (MQ=255)
                                                                                                                    acaAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAAc     >  1:222722/1‑141 (MQ=255)
                                                                                                                     caAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACt    <  2:103942/141‑1 (MQ=255)
                                                                                                                       aGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAAAGCACTTATTAATTAACTaa  >  2:263576/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                              
TTGAATTTTGACTAAAATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTAC  >  CP000730/790645‑790904

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: