Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,335 G→C 11.1% intergenic (+50/‑180) nadE → / → cho NAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,823,3350GC11.1% 101.0 / 7.8 45intergenic (+50/‑180)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (16/23);  new base C (3/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

AACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCT  >  NC_000913/1823191‑1823477
                                                                                                                                                |                                                                                                                                              
aaCAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgc                                                                                                                                            <  2:251891/149‑1 (MQ=255)
 acagtgcGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATgcg                                                                                                                                           <  2:756775/143‑1 (MQ=255)
  cAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTACACGCCGGATgcgg                                                                                                                                          >  1:413443/1‑149 (MQ=255)
   aGGTCGCCATAACAATAGTGAACTGGTATCTGAAAACCGAACATAAACTCCGTCCGCCAATTACCGTTTTCGATGATTTCTTGAAAAAGTAATAATTTGCACATTTTGGATTGTTCGAAAAAGAGTAATTTGTTCACGCCGGATgcggc                                                                                                                                         <  2:1009450/149‑1 (MQ=255)
     gTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCGGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGGACTTATTGGATTGTGCGAAAAAGAGTAATTTGTTAACGACGGATGCGGCGt                                                                                                                                       >  1:907730/1‑149 (MQ=255)
        gCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTAACGCCGGATGCGGCGTGaa                                                                                                                                    >  1:868905/1‑149 (MQ=255)
        gCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACAAAAACGCCGGCCGCCAATAACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGAGTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGaa                                                                                                                                    <  2:357951/149‑1 (MQ=255)
         ccAGAACAATAGACAACTGGTATCAGAAAACCGACCATACACGCCGTCCGCCAATTACCGCTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAAc                                                                                                                                   <  2:286605/149‑1 (MQ=255)
          cAGAACAATAGAGAACTGGTATCTGAAAACCTAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACg                                                                                                                                  <  1:744363/149‑1 (MQ=255)
                aaTAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAt                                                                                                                            >  1:535716/1‑149 (MQ=255)
                     agaACTGGTATCTGATAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                     <  2:148772/135‑1 (MQ=255)
                     agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGAc                                                                                                                       <  1:897980/149‑1 (MQ=255)
                     agaACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGCATCCGGCGTGa                                                                                                                                     >  1:148772/1‑135 (MQ=255)
                                cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                            >  1:329826/1‑149 (MQ=255)
                                cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                            >  2:766939/1‑149 (MQ=255)
                                cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTACACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                            <  2:413443/149‑1 (MQ=255)
                                       ccGAACATACACGCCCTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGAGTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAt                                                                                                     <  1:187341/149‑1 (MQ=255)
                                       ccGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACTTATTGGATTGTGCGAAAAAGATTATTTTTTTCACCCCGTATGCGCCTTGAACGCCTTATTCTACATATAAAAATATGCAAAt                                                                                                     >  2:541671/1‑149 (MQ=255)
                                           aCATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCaa                                                                                                 <  2:423354/149‑1 (MQ=255)
                                                  cGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATTATTTCCACATTTTGGATTGTGCGAAAAAGAGTAATTTGTTCACGGCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTg                                                                                          <  2:165507/149‑1 (MQ=255)
                                                    ccGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCa                                                                                        <  2:805794/149‑1 (MQ=255)
                                                        ccGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATTTTGGATTGTGCGAAAAAGCGTAATTTGTTCACGCGGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAg                                                                                    <  2:978052/149‑1 (MQ=255)
                                                         cGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGt                                                                                   <  2:826876/149‑1 (MQ=255)
                                                          gacAATTACCGTTTTCTATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                  <  1:543428/147‑1 (MQ=255)
                                                          gCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTc                                                                                  >  2:592370/1‑149 (MQ=255)
                                                                  ccGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGc                                                                          <  2:284828/149‑1 (MQ=255)
                                                                  ccGTTTTCGATGATTTCTGGAAAAAGTAATAATTAGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGAAGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTAAATAAATTGAAGGAGTCGAATAgtc                                                                          >  1:20951/1‑147 (MQ=38)
                                                                   cGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTTTTCACGCCTGATGCGGAGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTAAATACATTGCAGGACTCGAATAGGcc                                                                         >  2:827894/1‑149 (MQ=255)
                                                                      tttCGATGATTTCTGGAAAAAGTAATAAATAGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGa                                                                      <  1:213553/149‑1 (MQ=255)
                                                                           atgatTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGc                                                                 <  1:138351/149‑1 (MQ=255)
                                                                              atTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTa                                                              <  2:888186/149‑1 (MQ=255)
                                                                                     gAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTAAATACATTCCAGTAGTCGAATAGGCCTGACAGGCGTAGCACGTc                                                       >  2:788245/1‑149 (MQ=255)
                                                                                        aaaGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTAATAGACCTATAAAACTATGAAAATTCAATACAATGCAGGAGAAGAATAGGACTGAAAGGCATCGCACTTCAgg                                                    >  1:1029720/1‑148 (MQ=37)
                                                                                                     gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATTCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                       <  1:178194/149‑1 (MQ=255)
                                                                                                     gCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCttt                                       <  1:425596/149‑1 (MQ=255)
                                                                                                       acaTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGAGTAGCACGTCAGACGGAGTAACCTTTGt                                     >  1:698884/1‑149 (MQ=255)
                                                                                                        caTATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTc                                    >  2:784052/1‑149 (MQ=255)
                                                                                                              ggATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGAc                              >  2:885607/1‑149 (MQ=255)
                                                                                                                        aaaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCttt                    <  1:595120/149‑1 (MQ=255)
                                                                                                                         aaaaGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGGCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCtttt                   <  1:215481/149‑1 (MQ=255)
                                                                                                                                  aTTTGTTCACGCAGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGAGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCt          >  2:932981/1‑149 (MQ=255)
                                                                                                                                      gTTCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAAATCAATACAGTGCAGGAGTCGAATAGGCCTGACAGGCGTAGACCGTCAGACGGTGAAACGGTTGTCATCGCCCCGCCTCTGTTTTAATCGCTTccc      >  2:583684/1‑149 (MQ=255)
                                                                                                                                       ttCACGCCGCATCCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCGTTGTCATAGACCCGCCTCTATTTTAATCGCTTCCCg     >  2:3055/1‑149 (MQ=39)
                                                                                                                                         cACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGcc   >  2:392820/1‑149 (MQ=255)
                                                                                                                                          aCGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCt  <  2:735684/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
AACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCT  >  NC_000913/1823191‑1823477

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