Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,620,907 (A)8→9 intergenic (‑35/‑290) upp ← / → purM uracil phosphoribosyltransferase/phosphoribosylaminoimidazole synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,620,8991.A84.9% 97.7 / 10.7 33intergenic (‑27/‑298)upp/purMuracil phosphoribosyltransferase/phosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base A (13/15);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

TACCCACTTCGGAAGCGAGTTCGCGAAAGCGCTTGGTGCTGATATCTTGCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑AAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATGGTTATGAAAGAAGATTATTC  >  NC_000913/2620764‑2621028
                                                                                                                                        |                                                                                                                                   
tACCCACTTCGGAAGCGAGTTCGCGAAAGCGCTTGGTGCTGATATCTTGCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑aaa                                                                                                                                 >  1:258306/1‑139 (MQ=255)
tACCCACTTCGGAAGCGAGTTCGCGAAAGCGCTTGGTGCTGATATCTTGCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑aaa                                                                                                                                 >  1:665168/1‑139 (MQ=255)
        tCGGAAGCGAGTTCGCGAAAGCGCTTGGTGCTGATATCTTGCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑AAAAAAAATCg                                                                                                                        >  1:43037/1‑139 (MQ=255)
            aaGCGAGTTCGCGAAAGCGCTTGGTGCTGATATCTTGCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACg                                                                                                                     <  2:346424/139‑1 (MQ=255)
                                  ggTGCTGATATCTTGCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑AAAAAAAATCGACGGATTATACCTCCTTTCTTCAAgg                                                                                              <  2:43037/139‑1 (MQ=255)
                                                gCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑AAAGAAAAATCGAc                                                                                                                      >  2:417354/1‑102 (MQ=255)
                                                gCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTGGAGGGGCAGCCACAAAAAAAAATCGAc                                                                                                                      <  1:417354/102‑1 (MQ=255)
                                                 cTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCttt                                                                                >  1:265154/1‑139 (MQ=255)
                                                   cACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTTAAGGCGGCAATATTCtt                                                                                 >  2:74901/1‑136 (MQ=255)
                                                   cACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCtt                                                                                 <  1:74901/136‑1 (MQ=255)
                                                                       tgtTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATTCTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGt                                                                  <  1:23575/131‑1 (MQ=255)
                                                                       tgtTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGt                                                                  >  2:23575/1‑131 (MQ=255)
                                                                       tgtTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑aaaaaaaaatcga                                                                                                                      <  2:480672/78‑6 (MQ=255)
                                                                       tgtTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑aaaaaaaaatcga                                                                                                                      >  1:480672/1‑73 (MQ=255)
                                                                         tttGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATa                                                        <  2:665168/139‑1 (MQ=255)
                                                                          ttGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATaa                                                       <  2:265154/139‑1 (MQ=255)
                                                                                 gTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAgg                                                                                              >  1:594976/1‑93 (MQ=255)
                                                                                 gTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAgg                                                                                              <  2:594976/93‑1 (MQ=255)
                                                                                                         tCATACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGa                        <  2:258306/139‑1 (MQ=255)
                                                                                                           aTACTCTTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATg                      >  1:384619/1‑139 (MQ=255)
                                                                                                               tctTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATGGTTa                  <  2:384619/139‑1 (MQ=255)
                                                                                                                ctTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTc                                                                                   <  2:31065/73‑1 (MQ=255)
                                                                                                                ctTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTc                                                                                   >  1:31065/1‑73 (MQ=255)
                                                                                                                ctTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCaa                                                               <  1:554084/93‑1 (MQ=255)
                                                                                                                ctTTCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCaa                                                               >  2:554084/1‑93 (MQ=255)
                                                                                                                 tttCTCCTTTGAGGGGCAGCCACAAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAg                         <  2:180625/131‑1 (MQ=255)
                                                                                                                 tttCTCCTTTGAGGGGCAGCCACAAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAg                         >  1:180625/1‑131 (MQ=255)
                                                                                                                 tttCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATGGTTATg                >  2:659528/1‑139 (MQ=255)
                                                                                                                 tttCTCCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATGGTTATg                >  2:396664/1‑139 (MQ=255)
                                                                                                                    ctcCTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATGGTTATGaaa             >  1:551911/1‑139 (MQ=255)
                                                                                                                       cTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCa                                                                                         >  2:660758/1‑60 (MQ=255)
                                                                                                                       cTTTGAGGGGCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCa                                                                                         <  1:660758/60‑1 (MQ=255)
                                                                                                                              gggCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATa                                           <  2:398260/99‑1 (MQ=255)
                                                                                                                              gggCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATa                                           >  1:398260/1‑99 (MQ=255)
                                                                                                                               ggCAGCCACAAAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATGGTTATGAAAGAAGATTATTc  <  2:551911/139‑1 (MQ=255)
                                                                                                                                   gCCAC‑‑AAAAAAAATCGACGGATTATACCTCCTTTCTTCAAgg                                                                                              >  2:564164/1‑42 (MQ=255)
                                                                                                                                   gCCAC‑‑AAAAAAAATCGACGGATTATACCTCCTTTCTTCAAgg                                                                                              <  1:564164/42‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                   
TACCCACTTCGGAAGCGAGTTCGCGAAAGCGCTTGGTGCTGATATCTTGCTCACGCATCAGTCCCAGCTTGTGTTTGACGAGTGGGTGTTTGACTTCCACGATCTTCATACTCTTTCTCCTTTGAGGGGCAGCCAC‑‑AAAAAAAATCGACGGATTATACCTCCTTTCTTCAAGGCGGCAATATTCTTTTCGTTGACTTTAGTCAAAATGATAACGGTTTGAGATAAAGTTATTTTATATTCAGATGGTTATGAAAGAAGATTATTC  >  NC_000913/2620764‑2621028

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: