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,7780GA97.0% 103.5 / ‑3.6 33intergenic (+277/‑197)kipA/ltaSpossible allophanate hydrolase subunit 2/possible phosphatidylglycerol‑‑membrane‑oligosaccharide glycerophosphotransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (17/15);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

ATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATAATGTTAA  >  CP000730/790661‑790914
                                                                                                                     |                                                                                                                                        
attattTTATATTTATTAATTGATCTATGCTTATTATTACAATTTGATTACAAATTGTAAATTGGTTAACTGAATGATAATATTAAATAAAGAAACTCTCACAATCAAATATGAGATATAGCCCAAAATACTTGTTAAATc                                                                                                                   <  2:38987/141‑1 (MQ=255)
 ttACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCa                                                                                                                  >  2:194310/1‑141 (MQ=255)
              attaATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTa                                                                                                     >  2:361972/1‑141 (MQ=255)
               ttaATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTAc                                                                                                    <  1:167387/141‑1 (MQ=255)
                          aTGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGaaa                                                                                         <  2:112748/141‑1 (MQ=255)
                          aTGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGaaa                                                                                         >  1:101262/1‑141 (MQ=255)
                              ttattattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATta                                                                                     >  2:84291/1‑141 (MQ=255)
                                attattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATtata                                                                                   <  1:282409/141‑1 (MQ=255)
                                attattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATtata                                                                                   <  1:51155/141‑1 (MQ=255)
                                  tattaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATaaa                                                                                 >  2:203765/1‑141 (MQ=255)
                                   attaCAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTATACAAGCAAATATGAGTTATAGCCAAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATaaaa                                                                                >  2:356301/1‑140 (MQ=255)
                                                  cAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTa                                                                 >  2:144587/1‑141 (MQ=255)
                                                   aaaTTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTaa                                                                >  2:290723/1‑141 (MQ=255)
                                                       tttAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTAc                                                            >  2:142657/1‑141 (MQ=255)
                                                           aaTTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACtata                                                        >  1:88758/1‑141 (MQ=255)
                                                           aaTTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACtata                                                        >  1:124337/1‑141 (MQ=255)
                                                              ttGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAAtt                                                     <  2:157599/141‑1 (MQ=255)
                                                                     ttGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACaa                                              <  2:172233/141‑1 (MQ=255)
                                                                       gAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATg                                            >  1:72341/1‑141 (MQ=255)
                                                                       gAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATg                                            <  1:287771/141‑1 (MQ=255)
                                                                         aTGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGtt                                          >  1:44785/1‑141 (MQ=255)
                                                                             taataTTAAATAAAGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACt                                      >  1:123135/1‑141 (MQ=255)
                                                                                        aaaGAAACTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAAGCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCAc                           <  1:144587/141‑1 (MQ=255)
                                                                                            aaaCTTACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACttat                       <  2:111213/141‑1 (MQ=255)
                                                                                                 tACACAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTattaatt                  <  1:189127/141‑1 (MQ=255)
                                                                                                  acacAAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTattaatta                 >  2:298510/1‑141 (MQ=255)
                                                                                                     caAGCAAATATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACt              <  2:316353/141‑1 (MQ=255)
                                                                                                           aaTATGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATaa        >  1:179869/1‑141 (MQ=255)
                                                                                                               tGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATAATGtt    >  2:312933/1‑141 (MQ=255)
                                                                                                               tGAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATAATGtt    <  1:361972/141‑1 (MQ=255)
                                                                                                                gAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATAATGTTa   <  2:15619/141‑1 (MQ=255)
                                                                                                                gAGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATAATGTTa   <  2:124337/141‑1 (MQ=255)
                                                                                                                 aGTTATAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATAATGTTaa  >  1:188760/1‑141 (MQ=255)
                                                                                                                     |                                                                                                                                        
ATTACTCTATATTTATTAATTGAGCTATGCTTATTATTACAATTTGATTACAAATTTTAAATTTGTTAATTGAATGATAATATTAAATAAAGAAACTTACACAAGCAAATATGAGTTGTAGCCCAAAATACTTGTTAAATCAAAGTTGAAAGCTACAAATAACGAAAATTATAAACTTGAATCTGAAAGTAATTACTATAATTATGACAATGTTAACTTTTAAACGCACTTATTAATTAACTACATAATGTTAA  >  CP000730/790661‑790914

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