Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,427 T→C intergenic (‑98/‑154) cysD ← / → iap sulfate adenylyltransferase subunit 2/alkaline phosphatase isozyme conversion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,876,4270TC97.5% 132.3 / ‑4.1 40intergenic (‑98/‑154)cysD/iapsulfate adenylyltransferase subunit 2/alkaline phosphatase isozyme conversion protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (23/16);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCAT  >  NC_000913/2876290‑2876562
                                                                                                                                         |                                                                                                                                       
ccTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTc                                                                                                                                       >  1:355127/1‑140 (MQ=255)
    cAGTTGCCGCAGGTGAGTAAGTCGTATTAGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTa                                                                                                                                   >  1:484820/1‑140 (MQ=255)
     aGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTaa                                                                                                                                  <  1:116040/140‑1 (MQ=255)
     aGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTaa                                                                                                                                  <  2:588879/140‑1 (MQ=255)
                   aGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGAttt                                                                                                                    >  2:799947/1‑140 (MQ=255)
                      aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGc                                                                                                                 <  2:355127/140‑1 (MQ=255)
                      aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGc                                                                                                                 <  2:470247/140‑1 (MQ=255)
                                             ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCtt                                                                                          >  1:256764/1‑140 (MQ=255)
                                                    gCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                                                                   <  2:850381/140‑1 (MQ=255)
                                                          ccGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAtt                                                                             >  2:373507/1‑140 (MQ=255)
                                                             cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                          >  1:26267/1‑140 (MQ=255)
                                                             cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                          >  2:782560/1‑140 (MQ=255)
                                                             cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                                                                          >  2:295663/1‑140 (MQ=255)
                                                               aTTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCggg                                                                        <  1:53171/140‑1 (MQ=255)
                                                                       gCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAgg                                                                >  1:578327/1‑140 (MQ=255)
                                                                        ccATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGa                                                               <  1:782560/140‑1 (MQ=255)
                                                                        ccATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGa                                                               <  1:295663/140‑1 (MQ=255)
                                                                            cAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCgg                                                                         <  2:55930/126‑1 (MQ=255)
                                                                            cAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCgg                                                                         >  1:55930/1‑126 (MQ=255)
                                                                                       aCTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGtt                                                >  2:219691/1‑140 (MQ=255)
                                                                                             gggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAgcg                                          >  1:154528/1‑140 (MQ=255)
                                                                                              ggAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAgcgc                                         <  1:799947/140‑1 (MQ=255)
                                                                                                 gCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGc                                      >  1:645428/1‑140 (MQ=255)
                                                                                                  cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                     >  1:280296/1‑140 (MQ=255)
                                                                                                  cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                     >  2:69738/1‑140 (MQ=255)
                                                                                                  cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                     >  2:486441/1‑140 (MQ=255)
                                                                                                  cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCg                                     >  2:427675/1‑140 (MQ=255)
                                                                                                   gTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCgg                                    >  1:160733/1‑140 (MQ=255)
                                                                                                   gTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCgg                                    <  2:645428/140‑1 (MQ=255)
                                                                                                   gTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCgg                                    >  1:571622/1‑140 (MQ=255)
                                                                                                     aaGAGAACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                                  >  1:567034/1‑140 (MQ=255)
                                                                                                        agaACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                               <  2:26267/140‑1 (MQ=255)
                                                                                                        agaACGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGtt                               <  2:154528/140‑1 (MQ=255)
                                                                                                            cGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAg                                 >  1:587417/1‑134 (MQ=255)
                                                                                                            cGAATGAAATTACCAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAg                                 <  2:587417/134‑1 (MQ=255)
                                                                                                                         cAATTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTaa              >  2:167734/1‑140 (MQ=255)
                                                                                                                           aTTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCg                  >  1:601163/1‑134 (MQ=255)
                                                                                                                           aTTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCg                  <  2:601163/134‑1 (MQ=255)
                                                                                                                           aTTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCg                  <  2:18506/134‑1 (MQ=255)
                                                                                                                           aTTAGAATGAGTAGCTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCg                  >  1:18506/1‑134 (MQ=255)
                                                                                                                                    agtagcTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa   >  1:73874/1‑140 (MQ=255)
                                                                                                                                    agtagcTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa   <  2:567034/140‑1 (MQ=255)
                                                                                                                                    agtagcTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa   <  2:208154/140‑1 (MQ=255)
                                                                                                                                    agtagcTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCa   >  2:756507/1‑140 (MQ=255)
                                                                                                                                     gtagcTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCAt  >  1:502085/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCAT  >  NC_000913/2876290‑2876562

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