Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,424 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,4240GA97.4% 107.7 / NA 38intergenic (‑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/14);  minor base C (1/0);  total (24/14)
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/2876285‑2876555
                                                                                                                                           |                                                                                                                                   
tcGCCTCCTGTTGTCTCAGGTGTGTTGGTCGTTTTTGATCCTTTACCGTTCCTTTGCTATTTCGCTATTTTGTTGCCATCAGATGTTTCGACTATAGGGAGCGTAATAGAACGAATGTATTTTCCAATTAGAATGAGTAATTc                                                                                                                                  <  1:574915/142‑1 (MQ=255)
 cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGCTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGATTACTTcc                                                                                                                                 >  2:644915/1‑143 (MQ=255)
   ccTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGCTCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTAGAGGGGGCGTGAGAGAACGAATGCAATCACCAATTAGAACGAGCGATTCCtt                                                                                                                               >  1:352638/1‑143 (MQ=255)
       cAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACg                                                                                                                           <  2:447588/143‑1 (MQ=255)
             ccGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATCACCAATTAGAATGAGTAATTCCTTAACGGAATaa                                                                                                                     >  2:431266/1‑143 (MQ=255)
                        tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCaa                                                                                                          >  2:416689/1‑143 (MQ=255)
                           gtaggaTTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGc                                                                                                       <  1:32114/138‑1 (MQ=255)
                            tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAATGAAATAACGATTTGGCAAAGCt                                                                                                      >  1:436148/1‑143 (MQ=255)
                            tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACAAATTAGAATGAGTAATTCCTTAACGGAATAACGCTTTGGCGAACCt                                                                                                      >  1:24292/1‑143 (MQ=255)
                                       ccATATCCGTTCCTTTGCAATACCGCTATTTTCTCGCCATCAGACGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCaaaaa                                                                                           <  2:1751/143‑1 (MQ=255)
                                        cATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAg                                                                                          >  2:348821/1‑143 (MQ=255)
                                                ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                                                                  >  2:502836/1‑143 (MQ=255)
                                                    tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACAGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                                                              >  1:41732/1‑143 (MQ=255)
                                                       gCAATTCCGCTATTTTCTTGCCATCAGATGTTTCGACCATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCgg                                                                           <  2:264855/143‑1 (MQ=255)
                                                        cAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCCATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCCAAAAGTGCTTAAGGCACCGGa                                                                          >  1:630616/1‑143 (MQ=255)
                                                            aCCGCTATTGTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTc                                                                      <  2:384393/143‑1 (MQ=255)
                                                             tcgttATTTGTTTGCGCTCAGTCGTTCCGGCTATCGGGAGCGTAAGTGAACGTATGAAATTACCAATTAGAATCCTTCATTCGTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                     <  1:634315/139‑1 (MQ=255)
                                                                cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGCGTGAGTAATTCCTTAACGGAATAACGACTTGGCAAAGCTAATATGAAAAAGTGCTTGAGACACCGGATTTCGGGc                                                                  >  2:202839/1‑143 (MQ=255)
                                                                     ttCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAATTCGGGCGttca                                                             >  1:101903/1‑141 (MQ=255)
                                                                             aTCAGATGTTTCGACTATAGGGAGCGGAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAAGGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAtt                                                     <  2:144081/143‑1 (MQ=255)
                                                                                aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGGCTTCGGGCGTTCAGGAAGATTTGa                                                  >  1:53569/1‑143 (MQ=255)
                                                                                aGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGa                                                  >  2:485762/1‑143 (MQ=255)
                                                                                       tCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGAGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATGTGAAATTGtt                                           >  2:344035/1‑143 (MQ=255)
                                                                                        cGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGGTTTGACATTGttt                                          >  2:563411/1‑143 (MQ=255)
                                                                                           cTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATATCGGGCGTCTAGGAAGAGTTGAAATCGTTTTAg                                       >  1:123859/1‑143 (MQ=255)
                                                                                                ggggggggaaGAGAACGAATGAAATTAGCAATGAGAATGAGCAATTCCTTAACGGAGTAACGGTTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAg                                  <  1:623616/135‑1 (MQ=255)
                                                                                                     cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                             >  2:69751/1‑143 (MQ=255)
                                                                                                        aaGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                          >  2:425834/1‑143 (MQ=255)
                                                                                                         agagAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGttt                         >  2:13906/1‑143 (MQ=255)
                                                                                                         agagAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGttt                         >  1:220188/1‑143 (MQ=255)
                                                                                                         agagAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGAGGCCGttt                         >  1:107867/1‑143 (MQ=255)
                                                                                                            gaACGAATGAAATCACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCAt                      <  1:431266/143‑1 (MQ=255)
                                                                                                             aaCGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGACAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                     >  1:455742/1‑143 (MQ=255)
                                                                                                                  aTGAATTTACCAGTTGGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATg                <  1:516072/143‑1 (MQ=255)
                                                                                                                    gAAATTACCAATGAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGc              <  2:320912/143‑1 (MQ=255)
                                                                                                                         tACTAAGTAGAATGAGTAATTCCTTAGCGGAATAACGATTTGGCAAAGCTAGTATCAAAAAGTGCTTAAGGCTCCGGAGGTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaa         <  1:110948/143‑1 (MQ=255)
                                                                                                                             aaTTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaaaaaa     >  1:174114/1‑143 (MQ=255)
                                                                                                                                tAGAATGAGTAATTCCTTAACTGATTAATGATATGGAAAAGCTAACATAAAAAAGTGCTTAAGGCACCGGATGTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAAtt  <  2:281636/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
CCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATT  >  NC_000913/2876285‑2876555

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: ‑