Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,485,872 C→T R27H (CGC→CAC)  pepA ← multifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,485,8720CT97.5% 153.1 / ‑6.1 40R27H (CGC→CAC) pepAmultifunctional aminopeptidase A: a cyteinylglycinase, transcription regulator and site‑specific recombination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/17);  minor base G (1/0);  total (23/17)
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

ATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAA  >  NC_000913/4485725‑4486016
                                                                                                                                                   |                                                                                                                                                
aTAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGAc                                                                                                                                               >  1:111817/1‑151 (MQ=255)
         aTTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGa                                                                                                                                      <  1:276127/151‑1 (MQ=255)
            cGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAg                                                                                                                                    <  1:192971/150‑1 (MQ=255)
              cTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACg                                                                                                                                 >  2:19499/1‑151 (MQ=255)
                        aCATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCacgacg                                                                                                                        <  2:56425/150‑1 (MQ=255)
                         cATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGAt                                                                                                                      >  1:538717/1‑151 (MQ=255)
                         cATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGAt                                                                                                                      >  1:424932/1‑151 (MQ=255)
                          aTTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATg                                                                                                                     >  2:14142/1‑151 (MQ=255)
                              ggAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAgg                                                                                                                 <  1:532652/151‑1 (MQ=255)
                               gAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                                                                                                >  2:246344/1‑151 (MQ=255)
                               gAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGc                                                                                                                <  2:632808/151‑1 (MQ=255)
                                       tGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCt                                                                                                        >  2:516016/1‑151 (MQ=255)
                                              acaaTGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTtctc                                                                                                 >  1:361568/1‑151 (MQ=255)
                                                   gTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGc                                                                                            >  2:314149/1‑151 (MQ=255)
                                                     cTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCt                                                                                           >  2:512161/1‑150 (MQ=255)
                                                      tGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTAc                                                                                         >  1:346020/1‑151 (MQ=255)
                                                       gCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCt                                                                                           <  1:644783/148‑1 (MQ=255)
                                                               tttCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCttttt                                                                                 <  2:347524/150‑1 (MQ=255)
                                                                  ccttccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac                                                                             >  1:557422/1‑151 (MQ=255)
                                                                  ccttccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTacac                                                                             <  2:245308/151‑1 (MQ=255)
                                                                      ccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaaa                                                                         <  1:362160/151‑1 (MQ=255)
                                                                      ccAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTaa                                                                          <  2:552523/150‑1 (MQ=255)
                                                                              ccccGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGc                                                                 >  1:424113/1‑151 (MQ=255)
                                                                                 cGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACt                                                              >  2:173269/1‑151 (MQ=255)
                                                                                                cTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAGGTGGACGTGGGTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACGACGCGCCTGAATCtt                                               >  2:115995/1‑151 (MQ=255)
                                                                                                  gATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTaa                                             <  2:84760/151‑1 (MQ=255)
                                                                                                   aTGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTaa                                             >  1:438109/1‑150 (MQ=255)
                                                                                                     gTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGTATCTTAAAg                                           <  1:520841/150‑1 (MQ=255)
                                                                                                       aCCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACa                                        <  1:4323/151‑1 (MQ=255)
                                                                                                          cATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAAcg                                     >  2:122380/1‑151 (MQ=255)
                                                                                                                  aTTTTATCGAGCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGc                             <  2:425968/151‑1 (MQ=255)
                                                                                                                            gCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGAt                      >  2:162677/1‑148 (MQ=255)
                                                                                                                            gCTGTTCTGCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAg                    <  1:449926/150‑1 (MQ=255)
                                                                                                                             cTGTTCTGCAATCGGAGAAAGGTGACGAGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGa                   >  2:573960/1‑150 (MQ=255)
                                                                                                                                    gCAATCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCCCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAg           <  2:340301/151‑1 (MQ=255)
                                                                                                                                      aaTCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCt         >  2:457297/1‑151 (MQ=255)
                                                                                                                                        tCGGAGAAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCttt       >  2:219841/1‑151 (MQ=255)
                                                                                                                                           gagaAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCttttg     >  1:428578/1‑148 (MQ=255)
                                                                                                                                            agaAAGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTa   <  2:84263/151‑1 (MQ=255)
                                                                                                                                              aaaGGTGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTaa  >  1:402834/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATAAGGAGAATTCGCTCGGAAAGTACATTCGGAACATGGTGCAGCAACAATGTCTGCCCCGGTTTTCCTTCCAGTTCGCCCCGACGTAGCAGGGCGCTGATGTACCCATCGCTGATTTTATCGAGCTGTTCTGCAATCGGAGAAAGGCGACGTGGTTCGAAGACGCCCACGACGATGCAGGCACTCCGCTGTTTCTCCGGGCTACCGCTTTTTACACTAAACTCCATGCACTACGCTCCTGAATCTTAAAGACAACGGCGGTGGCTACAGATAGAATTGCAAGCTTTCGTAA  >  NC_000913/4485725‑4486016

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