Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,950,559 (C)5→6 intergenic (‑37/‑273) aspS ← / → yecD aspartyl‑tRNA synthetase/isochorismatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,950,5541.C97.7% 155.2 / ‑3.5 44intergenic (‑32/‑278)aspS/yecDaspartyl‑tRNA synthetase/isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (22/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTACC  >  NC_000913/1950418‑1950684
                                                                                                                                         |                                                                                                                                    
gAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑‑cc                                                                                                                                   >  1:511069/1‑139 (MQ=255)
 aaGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑‑cc‑c                                                                                                                                 <  2:1164255/139‑1 (MQ=255)
        ggCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑‑cc‑cccctgct                                                                                                                          <  2:1176404/139‑6 (MQ=255)
                   tCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGAc                                                                                                                >  2:787413/1‑139 (MQ=255)
                    cACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCAGCGGATG‑‑CCGCCCTGCTGCGCAGGTGAcc                                                                                                               <  1:877103/139‑1 (MQ=255)
                     acgacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCa                                                                                                              >  1:729704/1‑139 (MQ=255)
                       gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAg                                                                                                            >  2:325578/1‑139 (MQ=255)
                       gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCC‑‑‑‑CTTAGCTGCCGGCGGATG‑‑cc‑cccc                                                                                                                              <  1:245784/116‑2 (MQ=255)
                       gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCC‑‑‑‑CTTAGCTGCCGGCGGATG‑‑cc‑cccc                                                                                                                              >  2:245784/1‑115 (MQ=255)
                                       aCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑‑cc‑cccctg                                                                                                                            <  2:48347/106‑4 (MQ=255)
                                       aCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑‑cc‑cccctg                                                                                                                            >  1:48347/1‑103 (MQ=255)
                                                  tCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACt                                                                                 >  1:1065170/1‑139 (MQ=255)
                                                  tCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACt                                                                                 >  2:922353/1‑139 (MQ=255)
                                                  tCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACt                                                                                 >  1:581325/1‑139 (MQ=255)
                                                   cACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtg                                                                                >  1:634229/1‑139 (MQ=255)
                                                                      gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCGTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   >  2:927796/1‑133 (MQ=255)
                                                                      gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   <  1:927796/133‑1 (MQ=255)
                                                                      gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtat                                                             >  2:1013251/1‑139 (MQ=255)
                                                                                 cTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                     >  2:858346/1‑88 (MQ=255)
                                                                                 cTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTACCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                     <  1:858346/88‑1 (MQ=255)
                                                                                       acaATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGc                                            <  1:188157/139‑1 (MQ=255)
                                                                                         aaTATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCg                                          >  2:1074667/1‑139 (MQ=255)
                                                                                              tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   >  1:836205/1‑109 (MQ=255)
                                                                                              tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                   <  2:836205/109‑1 (MQ=255)
                                                                                                tGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                   <  1:1013251/139‑1 (MQ=255)
                                                                                                       aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                   <  1:406484/132‑1 (MQ=255)
                                                                                                       aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                   >  2:406484/1‑132 (MQ=255)
                                                                                                       aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                            >  2:207901/1‑139 (MQ=255)
                                                                                                          agaTATCCCTTAACTTAGCTGCCGGTGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGc                         >  1:3583/1‑139 (MQ=255)
                                                                                                           gaTATCCCTTAACTTAGCTGCCGGTGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGcc                        <  2:3583/139‑1 (MQ=255)
                                                                                                           gaTATCCCTTAACTTAGCTGCCGGCGGATGCCCC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGc                         <  2:267629/139‑1 (MQ=255)
                                                                                                           gaTATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGcc                        <  2:88608/139‑1 (MQ=255)
                                                                                                           gaTATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGcc                        <  2:511069/139‑1 (MQ=255)
                                                                                                           gaTATCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGcc                        <  1:1074667/139‑1 (MQ=255)
                                                                                                            atatCCCTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGccc                       >  2:467612/1‑139 (MQ=255)
                                                                                                                cccTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCg                                          >  2:247937/1‑116 (MQ=255)
                                                                                                                cccTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCg                                          <  1:247937/116‑1 (MQ=255)
                                                                                                                 ccTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                                <  1:803809/61‑1 (MQ=255)
                                                                                                                 ccTTAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                                >  2:803809/1‑61 (MQ=255)
                                                                                                                    tAACTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGa               <  1:325578/139‑1 (MQ=255)
                                                                                                                     aaCTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             >  1:924112/1‑60 (MQ=39)
                                                                                                                     aaCTTAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                             <  2:924112/60‑1 (MQ=39)
                                                                                                                         tAGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGcc          >  2:124612/1‑139 (MQ=255)
                                                                                                                          aGCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCt         >  2:793808/1‑139 (MQ=255)
                                                                                                                           gCTGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc        <  2:581325/139‑1 (MQ=255)
                                                                                                                             tGCCGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGt      <  1:467612/139‑1 (MQ=255)
                                                                                                                               ccGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa    <  2:729704/139‑1 (MQ=255)
                                                                                                                               ccGGCGGATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa    <  2:94832/139‑1 (MQ=255)
                                                                                                                                cggcggATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc   <  1:787413/139‑1 (MQ=255)
                                                                                                                                 ggcggATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAcc  <  1:922353/139‑1 (MQ=255)
                                                                                                                                 ggcggATGC‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAcc  >  2:1180133/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                    
GAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑‑CC‑CCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTACC  >  NC_000913/1950418‑1950684

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