Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,628,259 T→C T315T (ACT→ACC)  nadR → DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,2590TC96.3% 80.0 / ‑6.1 27T315T (ACT→ACC) nadRDNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/11);  minor base G (0/1);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTT  >  NC_000913/4628019‑4628498
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
gATCCTTGGGGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAGTATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGGCTTTTCACACCTCGGCGGAGATGAGATCGCATTGCAGTGTTCCGATTACGATAAAATCGCGCTGGGCCCCGCACAATACATTGATTTTGCAGTGAAATATGCCCATAGAGTGGCGTTTCTTGACACTGATTTTGTCACTCCCCAGGCGTTCt                                                                                                                                                                                                                                       >  2:9285/1‑251 (MQ=255)
                   aGCGGTAAATCCACCCTGGTACACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCg                                                                                                                                                                                                                    >  1:194884/1‑251 (MQ=255)
                    gCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCgt                                                                                                                                                                                                                   <  1:69449/251‑2 (MQ=255)
                       gTAAATCCACCCTTGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCgtgag                                                                                                                                                                                                                <  1:23220/251‑5 (MQ=255)
                        tAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCCGGCGGTCTGCAAAAAGTACGAAGGGCgtgagg                                                                                                                                                                                                               >  1:47777/1‑246 (MQ=255)
                                          caaacTTGCCAATATCTTCAATACCACCAGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCAACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGgcgt                                                                                                                                                                                             >  2:98432/1‑250 (MQ=255)
                                                                      aGTGCGTGGGAATATGGTCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCAACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGt                                                                                                                                                                 <  1:98432/251‑1 (MQ=255)
                                                                          cGTGGGAATATGGTCGCGATTATGTCTTTTCACTCCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGTGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATc                                                                                                                                                             <  1:131525/251‑1 (MQ=255)
                                                                             gggAATATTGCCGCGATTATGTCTTGTCACACCTCGGCTGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATTCCAATAAAGTGGCGTTTATTGACGCTTATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCTTCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCctg                                                                                                                                                          <  1:154004/251‑1 (MQ=255)
                                                                                               aTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGGGAACAACACGCCg                                                                                                                                        >  1:114163/1‑251 (MQ=255)
                                                                                                    ttttCACACCTCGGCGGTGATGAGATCGCCTTGCAGTATTTCGATTACGATAACATCGCGCTGGGCCCCGCACCATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCCCTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGGTTCGATCTGGGGATCCTGCTGGAGAACCACACGCCgtgggg                                                                                                                                   >  2:241303/1‑250 (MQ=255)
                                                                                                        cacaCCGCGGCGGTGATGAGATCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGc                                                                                                                                <  2:77642/250‑1 (MQ=255)
                                                                                                               cggcggTGATGAGGTCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGGGGCACGCCCCATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTTCAAAAAGTAAGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGCTGGtt                                                                                                                        >  2:98672/1‑251 (MQ=255)
                                                                                                                             tCGCATTGCAGTATTCCGATTACGATAAAATCGCGCTGTGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGATGGAGAACAACACGCCGTGGGTGGCGGTTGGTTTAAGCAGCCTctgc                                                                                                          >  2:5586/1‑248 (MQ=255)
                                                                                                                                         attccgattACGATAAAATCGCGCTGGGCCACGCACAATACAGTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGAt                                                                                              <  1:9285/243‑1 (MQ=255)
                                                                                                                                               attACGATAAAATCGCGCTGGGCCACGCACAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTAACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCaa                                                                                         >  1:138130/3‑250 (MQ=255)
                                                                                                                                                 gactaTAAAATCGCGCTGGGCCACGCACAATACATGGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACGGAGTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCGGCTGGAGAACAACACGCCGTGGGTGGCGGAAGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAaga                                                                                      <  2:245919/247‑1 (MQ=255)
                                                                                                                                                       aaaaTCGCGCTGGGCCACGAACAATACATTGATTTTGCATTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACGCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTTATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGTCTCTGCAGTTCGGTGGATCGCAAACAGTTcc                                                                                 <  2:166949/250‑1 (MQ=255)
                                                                                                                                                                            cAATACATTGATTTTGCAGTGACATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACCCGCCGTGGGGGGCGGATGGTTTACGCAGCCTCGGAAGTTCGGTGGATCGCAACGAGTTCCAGAACTTGCTGGGGGAGAGGCt                                                           >  2:192368/1‑251 (MQ=255)
                                                                                                                                                                            cAATACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCt                                                           >  1:92297/1‑251 (MQ=255)
                                                                                                                                                                              aTACATTGATTTTGCAGTGAAATATGCCAATAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCtgg                                                         >  1:121243/1‑249 (MQ=255)
                                                                                                                                                                                                       gCCAATAAAGTGGCGTTTATTGACACTGATTTTTTCACTACCCAGGCGTTCTGCAAAAAGGACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGc                                 <  2:47777/250‑1 (MQ=255)
                                                                                                                                                                                                          aaTAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGt                             >  1:210024/1‑251 (MQ=255)
                                                                                                                                                                                                           aTAAAGTGGCGTTTATTGACACTGATTTTGTCACTACCCAGGCGGTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGtt                            <  2:114163/251‑1 (MQ=255)
                                                                                                                                                                                                                         attgacactgATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTGCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGAtt               <  1:98672/241‑1 (MQ=255)
                                                                                                                                                                                                                             acactgATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTAATGAATACCGTTTCGATCTGGTGATCCTGCTGGTGAACAACGCGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAagcgaacgattatgac          >  1:259727/6‑237 (MQ=255)
                                                                                                                                                                                                                                     tttGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGttt  >  2:70569/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                               
GATCCTTGGCGGCGAGTCGAGCGGTAAATCCACCCTGGTAAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTT  >  NC_000913/4628019‑4628498

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