Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 699,559 A→G 21.6% intergenic (‑382/+15) 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,5590AG21.6% 74.4 / 18.0 37intergenic (‑382/+15)asnB/umpHasparagine synthetase B/UMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (15/14);  new base G (4/4);  total (19/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.25e-01

ATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCA  >  NC_000913/699448‑699704
                                                                                                               |                                                                                                                                                 
aTCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                                   <  1:152191/128‑1 (MQ=255)
aTCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTc                                                                                                                                   >  2:152191/1‑128 (MQ=255)
                            atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                               >  4:260576/1‑88 (MQ=255)
                            atatatTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACg                                                                                                                                               <  3:260576/88‑1 (MQ=255)
                                         gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                             >  4:747424/1‑141 (MQ=255)
                                         gCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                             <  3:747424/141‑1 (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                                                                                                                        <  2:622875/84‑1 (MQ=255)
                                                       aTAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCg                                                                                                                        >  1:622875/1‑84 (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                                                                             >  4:905953/1‑108 (MQ=255)
                                                                          cAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCa                                                                             <  3:905953/108‑1 (MQ=255)
                                                                                                    aCCACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGa           >  3:709500/1‑148 (MQ=255)
                                                                                                     ccACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                                                                 >  3:1080641/1‑93 (MQ=255)
                                                                                                     ccACGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGatat                                                                 <  4:1080641/93‑1 (MQ=255)
                                                                                                     ccACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg                                                            >  1:1105047/1‑98 (MQ=255)
                                                                                                     ccACGTCCGCAGACGTGGTTGCTAATCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATAtcgtcg                                                            <  2:1105047/98‑1 (MQ=255)
                                                                                                      cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                                 <  1:951294/44‑1 (MQ=255)
                                                                                                      cACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAg                                                                                                                 >  2:951294/1‑44 (MQ=255)
                                                                                                       aCGTCTGCGGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGc        <  4:709500/148‑1 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                                                                   <  2:156263/38‑1 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTc                                                                                                                   >  1:156263/1‑38 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                                 <  3:776809/104‑1 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAac                                                 >  4:776809/1‑104 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                            >  4:447288/1‑125 (MQ=255)
                                                                                                          tCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGt                            <  3:447288/125‑1 (MQ=255)
                                                                                                             gCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCa  >  2:421100/1‑148 (MQ=255)
                                                                                                               |                                                                                                                                                 
ATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGCTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTGGCGTTTGTCATCAGAGCCAACCACGTCCGCAGACGTGGTTGCTATTCAGATAACGTCGATTTCAGCGACTGACGGGTAAATCCAGCTGGGGCGGAAAGGCATACTGTCGATATCGTCGAGCGACGAAACACCAGAAAGCACCAGAATCGTCTCCAGACCTGCCTGGAAGCCGGCCA  >  NC_000913/699448‑699704

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: