Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,876,388 C→T intergenic (‑59/‑193) 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,3880CT100.0% 105.1 / NA 32intergenic (‑59/‑193)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

TCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAT  >  NC_000913/2876253‑2876512
                                                                                                                                       |                                                                                                                            
tCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTa                                                                                                                           >  2:1027174/1‑139 (MQ=255)
  tGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAag                                                                                                                         >  1:581423/1‑139 (MQ=255)
          ccTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGaa                                                                                                                 <  2:581423/139‑1 (MQ=255)
                      tgtgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCaa                                                                                                     <  1:1027174/139‑1 (MQ=255)
                       gtgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAAt                                                                                                    <  2:243175/139‑1 (MQ=255)
                           aTGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGa                                                                                                <  2:1239081/139‑1 (MQ=255)
                              cTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATg                                                                                             >  2:1165102/1‑139 (MQ=255)
                                 tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGagt                                                                                          >  2:707285/1‑139 (MQ=255)
                                 tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGagt                                                                                          <  1:446944/139‑1 (MQ=255)
                                    gCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGagtagt                                                                                       >  2:843165/1‑139 (MQ=255)
                                        ccAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCt                                                                                   >  1:1223074/1‑139 (MQ=255)
                                           gTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTaa                                                                                <  1:578023/139‑1 (MQ=255)
                                                             gTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTATAATGAGTAGTTCCTTAACGGAATAACGATTTGGCa                                                              >  1:746913/1‑139 (MQ=255)
                                                               cGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaaa                                                            >  1:528848/1‑139 (MQ=255)
                                                                                  ttCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCt                                         >  1:910188/1‑139 (MQ=255)
                                                                                       ttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAgg                                    >  2:1177778/1‑139 (MQ=255)
                                                                                        tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCa                                                              >  2:1126011/1‑112 (MQ=255)
                                                                                        tGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCa                                                              <  1:1126011/112‑1 (MQ=255)
                                                                                             tACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTaa                                      <  2:500651/131‑1 (MQ=255)
                                                                                             tACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATAGCAAAAAGTGCTTaa                                      >  1:500651/1‑131 (MQ=255)
                                                                                                     ttttCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCggg                      <  1:707285/139‑1 (MQ=255)
                                                                                                     ttttCTTGCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCggg                      <  2:528848/139‑1 (MQ=255)
                                                                                                            gCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGag                                                                                           <  1:360176/63‑1 (MQ=255)
                                                                                                            gCCATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGag                                                                                           >  2:360176/1‑63 (MQ=255)
                                                                                                              cATCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGa             <  1:1165102/139‑1 (MQ=255)
                                                                                                               aTCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGag                                                                                           <  1:354868/60‑1 (MQ=255)
                                                                                                               aTCAGATGTTTCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGag                                                                                           >  2:354868/1‑60 (MQ=255)
                                                                                                                         tCGACTATAGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt  <  1:843165/139‑1 (MQ=255)
                                                                                                                               ataGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                        <  1:641218/111‑1 (MQ=255)
                                                                                                                               ataGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCg                        >  2:641218/1‑111 (MQ=255)
                                                                                                                                 aGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                                                                                        >  2:1021690/1‑45 (MQ=255)
                                                                                                                                 aGGGAGTGTAAGAGAACGAATGAAATTACCAATTAGAATGagtag                                                                                        <  1:1021690/45‑1 (MQ=255)
                                                                                                                                       |                                                                                                                            
TCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAT  >  NC_000913/2876253‑2876512

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