Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,810,304 G→A *160* (TAG→TAA coaD → pantetheine‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,810,3040GA96.5% 83.4 / ‑6.4 29*160* (TAG→TAAcoaDpantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (13/15);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTT  >  NC_000913/3810176‑3810436
                                                                                                                                |                                                                                                                                    
cTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGcc                                                                                                                            <  1:247926/139‑1 (MQ=255)
   aTGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGa                                                                                                                         <  1:695617/139‑1 (MQ=255)
    tGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGAt                                                                                                                        <  2:747644/139‑1 (MQ=255)
     gCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATg                                                                                                                       >  1:329136/1‑139 (MQ=255)
     gCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATg                                                                                                                       <  2:600715/139‑1 (MQ=255)
               gagTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATc                                                                                                             <  1:554546/139‑1 (MQ=255)
                       gTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGgcgcgc                                                                                                     >  1:807314/1‑139 (MQ=255)
                        tttATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCa                                                                                                    <  2:743157/139‑1 (MQ=255)
                          tATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCg                                                                                                                           <  1:493692/114‑1 (MQ=255)
                          tATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCg                                                                                                                           >  2:493692/1‑114 (MQ=255)
                            tctcTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGaa                                                                                                <  2:807314/139‑1 (MQ=255)
                            tctcTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGaa                                                                                                <  2:326365/139‑1 (MQ=255)
                             ctctTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAt                                                                                               <  2:329136/139‑1 (MQ=255)
                                          gTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGt                                                                                                                     <  1:421474/104‑1 (MQ=255)
                                          gTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGt                                                                                                                     >  2:421474/1‑104 (MQ=255)
                                          gTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACt                                                                                  >  1:670405/1‑139 (MQ=255)
                                              aaGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGcc                                                                               <  2:868968/138‑1 (MQ=255)
                                              aaGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGcc                                                                               >  1:868968/1‑138 (MQ=255)
                                                                      aTGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                                              >  2:993847/1‑99 (MQ=255)
                                                                      aTGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAAtt                                                                                              <  1:993847/99‑1 (MQ=255)
                                                                        gTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                                       <  1:639294/104‑1 (MQ=255)
                                                                        gTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTg                                                                                       >  2:639294/1‑104 (MQ=255)
                                                                            cccATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGttt                                                <  2:670405/139‑1 (MQ=255)
                                                                                       ccGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACaa                                     >  1:561962/1‑139 (MQ=255)
                                                                                                 tCCATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTAcc                           >  2:969739/1‑139 (MQ=255)
                                                                                                   cATCAGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGc                         >  1:112393/1‑139 (MQ=255)
                                                                                                       aGGCGCTGATGGCGAAGTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCacac                     >  1:259796/1‑139 (MQ=255)
                                                                                                               aTGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGACGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACg             >  1:577358/1‑139 (MQ=255)
                                                                                                                       gTTAGCGTAACGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCt     >  1:483712/1‑139 (MQ=255)
                                                                                                                         tagcgtaacgTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCttt   >  1:809317/1‑139 (MQ=255)
                                                                                                                          agcgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCtttt  <  2:259796/139‑1 (MQ=255)
                                                                                                                           gcgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGaa                                                                                                <  2:446469/44‑1 (MQ=255)
                                                                                                                           gcgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGaa                                                                                                >  1:446469/1‑44 (MQ=255)
                                                                                                                            cgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGt                                             >  1:884854/1‑94 (MQ=255)
                                                                                                                            cgtaacgtTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGt                                             <  2:884854/94‑1 (MQ=255)
                                                                                                                                |                                                                                                                                    
CTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTTGGTGAAAGAGGTGGCGCGCCATCAGGGCGATGTCACCCATTTCCTGCCGGAGAATGTCCATCAGGCGCTGATGGCGAAGTTAGCGTAGCGTTTATGCCGGATGGTATGCCATCCGGCGCGCATGAATTACTTCTGGCACTGCCGACAATAAAACGTTGCCCGCTGCGCATGTTTAGTCGCCACAATCGGCGTACCGCACACCCGACACGGCTCACCTTTT  >  NC_000913/3810176‑3810436

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