Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,426 G→A intergenic (‑97/‑155) cysD ← / → iap sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,4260GA97.3% 105.0 / NA 37intergenic (‑97/‑155)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (23/13);  minor base C (1/0);  total (24/13)
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

CCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATT  >  NC_000913/2876287‑2876557
                                                                                                                                           |                                                                                                                                   
tcGCCTCCTGTTGTCTCAGGTGTGTTGGTCGTTTTTGATCCTTTACCGTTCCTTTGCTATTTCGCTATTTTGTTGCCATCAGATGTTTCGACTATAGGGAGCGTAATAGAACGAATGTATTTTCCAATTAGAATGAGTAATTc                                                                                                                                  <  2:574915/142‑1 (MQ=255)
 cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGCTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGATTACTTcc                                                                                                                                 >  1:644915/1‑143 (MQ=255)
   ccTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGCTCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTAGAGGGGGCGTGAGAGAACGAATGCAATCACCAATTAGAACGAGCGATTCCtt                                                                                                                               >  2:352638/1‑143 (MQ=255)
       cAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACg                                                                                                                           <  1:447588/143‑1 (MQ=255)
             ccGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATCACCAATTAGAATGAGTAATTCCTTAACGGAATaa                                                                                                                     >  1:431266/1‑143 (MQ=255)
                        tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCaa                                                                                                          >  1:416689/1‑143 (MQ=255)
                           gtaggaTTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGc                                                                                                       <  2:32114/138‑1 (MQ=255)
                            tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAATGAAATAACGATTTGGCAAAGCt                                                                                                      >  2:436148/1‑143 (MQ=255)
                            tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACAAATTAGAATGAGTAATTCCTTAACGGAATAACGCTTTGGCGAACCt                                                                                                      >  2:24292/1‑143 (MQ=255)
                                       ccATATCCGTTCCTTTGCAATACCGCTATTTTCTCGCCATCAGACGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                           <  1:1751/143‑1 (MQ=255)
                                        cATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAg                                                                                          >  1:348821/1‑143 (MQ=255)
                                                ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                                                                  >  1:502836/1‑143 (MQ=255)
                                                    tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACAGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                                                              >  2:41732/1‑143 (MQ=255)
                                                       gCAATTCCGCTATTTTCTTGCCATCAGATGTTTCGACCATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCgg                                                                           <  1:264855/143‑1 (MQ=255)
                                                        cAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCCATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCCAAAAGTGCTTAAGGCACCGGa                                                                          >  2:630616/1‑143 (MQ=255)
                                                            aCCGCTATTGTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTc                                                                      <  1:384393/143‑1 (MQ=255)
                                                             tcgttATTTGTTTGCGCTCAGTCGTTCCGGCTATCGGGAGCGTAAGTGAACGTATGAAATTACCAATTAGAATCCTTCATTCGTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                     <  2:634315/139‑1 (MQ=255)
                                                                cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGCGTGAGTAATTCCTTAACGGAATAACGACTTGGCAAAGCTAATATGAAAAAGTGCTTGAGACACCGGATTTCGGGc                                                                  >  1:202839/1‑143 (MQ=255)
                                                                     ttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAATTCGGGCGttca                                                             >  2:101903/1‑141 (MQ=255)
                                                                             aTCAGATGTTTCGACTATAGGGAGCGGAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAAGGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAtt                                                     <  1:144081/143‑1 (MQ=255)
                                                                                aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGGCTTCGGGCGTTCAGGAAGATTTGa                                                  >  2:53569/1‑143 (MQ=255)
                                                                                aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGa                                                  >  1:485762/1‑143 (MQ=255)
                                                                                       tCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGAGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATGTGAAATTGtt                                           >  1:344035/1‑143 (MQ=255)
                                                                                        cGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGGTTTGACATTGttt                                          >  1:563411/1‑143 (MQ=255)
                                                                                           cTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATATCGGGCGTCTAGGAAGAGTTGAAATCGTTTTAg                                       >  2:123859/1‑143 (MQ=255)
                                                                                                     cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                             >  1:69751/1‑143 (MQ=255)
                                                                                                        aaGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                          >  1:425834/1‑143 (MQ=255)
                                                                                                         agagAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGttt                         >  1:13906/1‑143 (MQ=255)
                                                                                                         agagAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGttt                         >  2:220188/1‑143 (MQ=255)
                                                                                                         agagAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGAGGCCGttt                         >  2:107867/1‑143 (MQ=255)
                                                                                                            gaACGAATGAAATCACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCAt                      <  2:431266/143‑1 (MQ=255)
                                                                                                             aaCGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGACAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                     >  2:455742/1‑143 (MQ=255)
                                                                                                                  aTGAATTTACCAGTTGGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                <  2:516072/143‑1 (MQ=255)
                                                                                                                    gAAATTACCAATGAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGc              <  1:320912/143‑1 (MQ=255)
                                                                                                                         tACTAAGTAGAATGAGTAATTCCTTAGCGGAATAACGATTTGGCAAAGCTAGTATCAAAAAGTGCTTAAGGCTCCGGAGGTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaa         <  2:110948/143‑1 (MQ=255)
                                                                                                                             aaTTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaa     >  2:174114/1‑143 (MQ=255)
                                                                                                                                tAGAATGAGTAATTCCTTAACTGATTAATGATATGGAAAAGCTAACATAAAAAAGTGCTTAAGGCACCGGATGTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAtt  <  1:281636/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
CCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATT  >  NC_000913/2876287‑2876557

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