Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,835,958 G→A P490S (CCG→TCG)  hypF ← carbamoyltransferase HypF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,835,9580GA98.0% 161.8 / NA 49P490S (CCG→TCG) hypFcarbamoyltransferase HypF
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/23);  minor base C (0/1);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGGCCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACG  >  NC_000913/2835828‑2836082
                                                                                                                                  |                                                                                                                            
gCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATg                                                                                                                     <  1:840426/140‑1 (MQ=255)
  aGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGct                                                                                                                   >  2:1050072/1‑140 (MQ=255)
       aGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGcc                                                                                                              <  2:553977/140‑1 (MQ=255)
              cAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACAcg                                                                                                       >  1:332633/1‑140 (MQ=255)
                 gTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACAcgccg                                                                                                    <  1:271649/140‑1 (MQ=255)
                    cTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTg                                                                                                 >  1:152372/1‑140 (MQ=255)
                     tCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGc                                                                                                <  1:1050072/140‑1 (MQ=255)
                     tCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGc                                                                                                <  1:832672/140‑1 (MQ=255)
                       acaTTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCAt                                                                                              <  1:945990/140‑1 (MQ=255)
                          ttCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATggg                                                                                           <  1:628222/140‑1 (MQ=255)
                           tCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCAt                                                                                              >  2:120763/1‑136 (MQ=255)
                           tCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCAt                                                                                              <  1:120763/136‑1 (MQ=255)
                            cgcgATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCa                                                                                         >  1:517875/1‑140 (MQ=255)
                               gATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCatga                                                                                      <  2:1143416/140‑1 (MQ=255)
                                            cAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCa                                                                                         >  2:231792/1‑124 (MQ=255)
                                            cAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCa                                                                                         <  1:231792/124‑1 (MQ=255)
                                            cAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCg                                                                         >  2:536711/1‑140 (MQ=255)
                                              ggCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGtt                                                                       <  2:152372/140‑1 (MQ=255)
                                                         ccccACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAg                                                            >  1:129805/1‑140 (MQ=255)
                                                         ccccACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAg                                                            >  2:712533/1‑140 (MQ=255)
                                                          cccACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGa                                                           >  1:104902/1‑140 (MQ=255)
                                                                    gcCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCAc                                                 >  1:106817/1‑140 (MQ=255)
                                                                      cGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCAcgc                                               >  2:161162/1‑140 (MQ=255)
                                                                           tcCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCa                                          <  2:832101/140‑1 (MQ=255)
                                                                           tcCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCa                                          <  2:517875/140‑1 (MQ=255)
                                                                            cccccATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCAc                                         <  2:1103023/140‑1 (MQ=255)
                                                                                  taccgataccgCCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGCCCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActggctg                                   <  2:1201557/140‑1 (MQ=255)
                                                                                   accgataccgTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActggctgg                                  >  2:1122726/1‑140 (MQ=255)
                                                                                       ataccgTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActggctgg                                  <  2:196777/136‑1 (MQ=255)
                                                                                       ataccgTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActggctgg                                  >  1:196777/1‑136 (MQ=255)
                                                                                       ataccgTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGAc                              >  2:868161/1‑140 (MQ=255)
                                                                                       ataccgTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGAc                              >  2:1012029/1‑140 (MQ=255)
                                                                                            gTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCa                                          <  2:476704/123‑1 (MQ=255)
                                                                                            gTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCa                                          >  1:476704/1‑123 (MQ=255)
                                                                                            gTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGc                         >  1:157602/1‑140 (MQ=255)
                                                                                             tCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActggctg                                   >  1:74174/1‑129 (MQ=255)
                                                                                             tCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActggctg                                   <  2:74174/129‑1 (MQ=255)
                                                                                             tCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActg                                       <  2:1157694/125‑1 (MQ=255)
                                                                                             tCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCActg                                       >  1:1157694/1‑125 (MQ=255)
                                                                                             tCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGcc                        >  1:974385/1‑140 (MQ=255)
                                                                                               gAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGTATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCg                      >  2:560023/1‑140 (MQ=255)
                                                                                               gAGCGTCAAAGCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCg                      >  1:199835/1‑140 (MQ=255)
                                                                                                         gCAATGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGTATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTc            <  1:560023/140‑1 (MQ=255)
                                                                                                           aaTGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCAt          <  2:106817/140‑1 (MQ=255)
                                                                                                           aaTGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCAt          <  2:129805/140‑1 (MQ=255)
                                                                                                             tGACATCACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGc        <  1:536711/140‑1 (MQ=255)
                                                                                                                   cACCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACg  >  2:639310/1‑140 (MQ=255)
                                                                                                                    aCCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGt             >  1:414796/1‑128 (MQ=255)
                                                                                                                    aCCGCCATCCAGCGACCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGt             <  2:414796/128‑1 (MQ=255)
                                                                                                                                  |                                                                                                                            
GCAGGCAAGCCGCCCAGGTGCTCACATTCGCGATAGTTCACCCGCAGGCACTCGCCGCCCCACAAAGCGCCGTTCTCCCCCATACCGATACCGTCGAGCGTCAAAGCAATGACATCACCGCCATCCAGCGGCCACTGATGCTCTGCCAGACACGCCGCTGCATGGGCATGATGATGCAGTACCGTTTGCGTCGGCAGATTCATTTCACGCGCCCACTGGCTGGAGACATAGCCCGGATGCGCGTCATGCACAACG  >  NC_000913/2835828‑2836082

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