Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,556 C→T 22.2% intergenic (‑379/+18) asnB ← / ← umpH asparagine synthetase B/UMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913699,5560CT22.2% 72.2 / 19.0 36intergenic (‑379/+18)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base C (14/14);  new base T (4/4);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

ATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGC  >  NC_000913/699448‑699698
                                                                                                            |                                                                                                                                              
aTCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                             <  1:152191/128‑1 (MQ=255)
aTCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                             >  2:152191/1‑128 (MQ=255)
                            atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                         <  3:260576/88‑1 (MQ=255)
                            atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                         >  4:260576/1‑88 (MQ=255)
                                         gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                       <  3:747424/141‑1 (MQ=255)
                                         gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                       >  4:747424/1‑141 (MQ=255)
                                                     tGATAAGCGTAGCCCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATa                                                                                                                        >  1:993164/1‑80 (MQ=255)
                                                     tGATAAGCGTAGCCCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATa                                                                                                                        <  2:993164/80‑1 (MQ=255)
                                                      gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTgg                                                                                                                                      <  1:201512/65‑1 (MQ=21)
                                                      gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTgg                                                                                                                                      >  2:201512/1‑65 (MQ=21)
                                                      gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGt                                                                                                                                     <  3:219070/66‑1 (MQ=255)
                                                      gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCTGCGGACGTGGt                                                                                                                                     >  4:219070/1‑66 (MQ=21)
                                                      gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGa                                                                                                        >  1:230286/1‑95 (MQ=255)
                                                      gATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGa                                                                                                        <  2:230286/95‑1 (MQ=255)
                                                       aTAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCg                                                                                                                  >  1:622875/1‑84 (MQ=255)
                                                       aTAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCg                                                                                                                  <  2:622875/84‑1 (MQ=255)
                                                           gCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAgg                                                                         >  4:777557/1‑121 (MQ=255)
                                                           gCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAgg                                                                         <  3:777557/121‑1 (MQ=255)
                                                                        ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa                                                                                                                            <  4:724921/57‑1 (MQ=255)
                                                                        ggCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCa                                                                                                                            >  3:724921/1‑57 (MQ=255)
                                                                          cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                       <  3:905953/108‑1 (MQ=255)
                                                                          cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                       >  4:905953/1‑108 (MQ=255)
                                                                                                    aCCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGa     >  3:709500/1‑148 (MQ=255)
                                                                                                     ccACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                                                           <  4:1080641/93‑1 (MQ=255)
                                                                                                     ccACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                                                           >  3:1080641/1‑93 (MQ=255)
                                                                                                     ccACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg                                                      >  1:1105047/1‑98 (MQ=255)
                                                                                                     ccACGTCCGCAGACGTGGTTGCTAATCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg                                                      <  2:1105047/98‑1 (MQ=255)
                                                                                                      cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                           >  2:951294/1‑44 (MQ=255)
                                                                                                      cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                           <  1:951294/44‑1 (MQ=255)
                                                                                                       aCGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGc  <  4:709500/148‑1 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                                                             >  1:156263/1‑38 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                                                             <  2:156263/38‑1 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                           <  3:776809/104‑1 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                           >  4:776809/1‑104 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                      <  3:447288/125‑1 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                      >  4:447288/1‑125 (MQ=255)
                                                                                                            |                                                                                                                                              
ATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGC  >  NC_000913/699448‑699698

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