Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,474,941 T→C W241R (TGG→CGG)  pdxI → pyridoxine 4‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,474,9410TC96.8% 102.8 / ‑6.1 31W241R (TGG→CGG) pdxIpyridoxine 4‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/18);  minor base G (0/1);  total (12/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCG  >  NC_000913/1474809‑1475072
                                                                                                                                    |                                                                                                                                   
gatGCTTTGGCCCACGATGGCATTGCCTACGCGCCCTTCCTCCCGCTCCGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGtt                                                                                                                              <  2:222368/140‑1 (MQ=255)
 atGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAAAGCCGGTGGCACTGGCGCGGCTGTTa                                                                                                                             >  1:295709/1‑140 (MQ=255)
      ttGGCCCACGATGGCATCGCCTTCGTGCCGTTCTTGCCGTTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCg                                                                                                                        <  2:246460/140‑1 (MQ=255)
            cACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCAcc                                                                                                                  <  1:125243/140‑1 (MQ=255)
                    cATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGACTGTTACAGCGTTCACCGAATAttt                                                                                                          >  1:313536/1‑140 (MQ=255)
                        gCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACATCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCt                                                                                                      <  1:60483/140‑1 (MQ=255)
                        gCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCt                                                                                                      <  2:653468/140‑1 (MQ=255)
                        gCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCt                                                                                                      <  2:62348/140‑1 (MQ=255)
                        gCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCt                                                                                                      >  2:586485/1‑140 (MQ=255)
                        gCCTACGTGCCGTTCCTCCCGCTCGCGGGTTTTACACAGCTACAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCt                                                                                                      <  1:175204/140‑1 (MQ=255)
                        gCCTACGTGCCGTTCATCCCGCCCGGGGGTTTTACACAGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCt                                                                                                      <  1:276810/140‑1 (MQ=255)
                         ccTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCGGGCGCGGTTGTTACCGCGTTCACCGAATATTTTGCTg                                                                                                     >  1:260420/1‑140 (MQ=255)
                                  cGTTCTTCCCGCTCGAGGGCTTTACCCCGTTGCAATCATCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAggg                                                                                            <  1:186644/140‑1 (MQ=255)
                                           cGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCg                                                                                   >  1:487118/1‑140 (MQ=255)
                                                  gggCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAGTATATTGCTGATGCCAGGTACGTCTTCGGTTgcgc                                                                            >  2:302249/1‑140 (MQ=255)
                                                                               tttCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGaaaaa                                               <  2:393471/140‑1 (MQ=255)
                                                                                 tCCGATGTTGCTGCGTGCCTGGGTGCAACTCCAATGCTGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAAtt                                             <  1:302249/140‑1 (MQ=255)
                                                                                 tCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAAtt                                             >  1:402797/1‑140 (MQ=255)
                                                                                   cGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGc                                           >  2:329289/1‑140 (MQ=255)
                                                                                    gATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCa                                          >  2:510445/1‑140 (MQ=255)
                                                                                     aTGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCAt                                         <  2:402797/140‑1 (MQ=255)
                                                                                          gctgcGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTc                                    <  1:510445/140‑1 (MQ=255)
                                                                                               gAGCCTGGGTGCAACCCCAATGGAGGTGGCGCTGGCGGGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAgg                               <  1:328667/140‑1 (MQ=255)
                                                                                               gAGCCTGGGTGCAACAGCAATGCAGGTGGCGGTGGTGCGACTGTTACACCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAgg                               <  2:295709/140‑1 (MQ=255)
                                                                                                          cAACACCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCt                    <  2:342566/140‑1 (MQ=255)
                                                                                                             cacCAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACg                 <  1:55251/140‑1 (MQ=255)
                                                                                                               ccAATGCAGGTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGtt               <  2:610952/140‑1 (MQ=255)
                                                                                                                       ggTGGCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTTCTGAAAAAATGCATCTATCTGAGGAAGTGTGGTCTACGTTGGATGGTa       >  2:268693/1‑140 (MQ=255)
                                                                                                                          ggCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACg                 <  1:354242/127‑1 (MQ=255)
                                                                                                                          ggCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACg                 >  2:354242/1‑127 (MQ=255)
                                                                                                                          ggCGCTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTAttt    >  1:466941/1‑140 (MQ=255)
                                                                                                                            cgcTGGCGCGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTcg  >  1:134623/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCG  >  NC_000913/1474809‑1475072

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