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,4260GA96.0% 70.2 / NA 26intergenic (‑97/‑155)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (13/11);  minor base C (0/1);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATG  >  NC_000913/2876291‑2876563
                                                                                                                                       |                                                                                                                                          
cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATGAGAATGAGTAATTCCTTa                                                                                                                                     >  1:453353/1‑143 (MQ=255)
  ccAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTCGAATGCGCAGTTCCTTaag                                                                                                                                   >  1:153040/1‑142 (MQ=255)
         ccGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATaa                                                                                                                            >  2:662401/1‑143 (MQ=255)
          cGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAGTGAGTAATTCCTTAACGGAATAAc                                                                                                                           >  1:222592/1‑143 (MQ=255)
               gTGAGAAACTCGTATATGATCCATACCCGTCCCCTTGCAATCCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGAttt                                                                                                                      <  1:172579/143‑1 (MQ=255)
                      aGTCGTATTTGATCCATAACCGTTCCTGTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATCGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAAGTCCTTAACGGAATAACGATTTGGCAAAg                                                                                                               <  1:612563/143‑1 (MQ=255)
                       gTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGc                                                                                                             <  1:425363/143‑1 (MQ=255)
                          gTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAACGTAAGAGAACAAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                                                                          >  1:305613/1‑143 (MQ=255)
                                  tCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTCTAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATATGGCAAAGCTAATATCaaaa                                                                                                  >  2:98901/1‑143 (MQ=255)
                                           gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTaa                                                                                         <  2:626920/143‑1 (MQ=255)
                                             tCCTTTGCAATACCGCTGTTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAgg                                                                                       <  1:153497/143‑1 (MQ=255)
                                                tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTACTATCAAACAGTGCTTAAGGCAc                                                                                    >  1:381844/1‑143 (MQ=255)
                                                          cGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGACCGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGCTTTCgg                                                                          >  2:75479/1‑143 (MQ=255)
                                                          cGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGTCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCgg                                                                          <  2:611204/143‑1 (MQ=255)
                                                           gCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCggg                                                                         >  1:394530/1‑143 (MQ=255)
                                                               ttttCTTGCCCTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGtt                                                                     <  2:464804/143‑1 (MQ=255)
                                                                       ccATCAGATCTTTCGACTATAGGGAGCGTAAGCGAACGAATGACATTACCAATTAGAATGAGTAATTCCTTAACTGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGa                                                             <  2:259997/143‑1 (MQ=255)
                                                                        cATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAt                                                            <  1:343224/143‑1 (MQ=255)
                                                                         aTCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGAtt                                                           >  1:445517/1‑143 (MQ=255)
                                                                               tGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt                                                     <  1:248744/143‑1 (MQ=255)
                                                                                        tataGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAGGGGACCGGAACACCAGCGTTTAGGAAGATTTGATATTGTTTTAgc                                            >  2:491997/1‑143 (MQ=255)
                                                                                                 cGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGTCGAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCa                                   >  1:627546/1‑143 (MQ=255)
                                                                                                     agagAACGAATGAAATTACCAATTAGAATGAGTAATGCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGTGCAGCGGCAGttc                               >  1:452042/1‑142 (MQ=255)
                                                                                                         aaCGAATGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGCAAAGTACAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATa                           <  1:695077/143‑1 (MQ=255)
                                                                                                             aaTGAAATTACCAATTAGAATGAGTAATTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACtct                       >  1:411383/1‑141 (MQ=255)
                                                                                                                                  gAGTACTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGACTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATg  <  1:517253/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
CTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTTTAGCGCAGCGGCAGTTTCATACTATGGCGGTAAAAAAATTTGCATG  >  NC_000913/2876291‑2876563

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