Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,823,338 G→C 11.5% intergenic (+53/‑177) 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,3380GC11.5% 92.8 / 5.0 42intergenic (+53/‑177)nadE/choNAD synthetase, NH3/glutamine‑dependent/endonuclease of nucleotide excision repair
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base C (3/2);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACT  >  NC_000913/1823191‑1823485
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                    >  2:766939/1‑149 (MQ=255)
                                cTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTa                                                                                                                    >  1:329826/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)
                                                                                                                                                  tGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACAGTGCAGGAGTCGAATAGGCCTGACAGGCGGAGCACGTCAGACGGGGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACt  <  2:79882/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACAGGTCGCCAGAACAATAGAGAACTGGTATCTGAAAACCGAACATAAACGCCGTCCGCCAATTACCGTTTTCGATGATTTCTGGAAAAAGTAATAATTTGCACATATTGGATTGTGCGAAAAAGAGTAATTTGTTCACGCCGGATGCGGCGTGAACGCCTTATTCGACCTATAAAACTATGCAAATTCAATACATTGCAGGAGTCGAATAGGCCTGACAGGCGTAGCACGTCAGACGGTGTAACCTTTGTCATCGACCCGCCTCTTTTTTAATCGCTTCCCGCCTGTTACACT  >  NC_000913/1823191‑1823485

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: ‑