Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,626,460 G→T 7.8% E38D (GAG→GAT)  PP_2298 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,626,4600GT 7.8% 159.1 / 3.9 52E38D (GAG→GAT) PP_2298hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/27);  new base T (2/2);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

AAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCATG  >  NC_002947/2626325‑2626598
                                                                                                                                       |                                                                                                                                          
aaaGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGCGGATCTCCTCGAGcaa                                                                                                                                <  5:140430/148‑3 (MQ=255)
aaaGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTgg                                                                                                                               <  2:52807/149‑1 (MQ=255)
    cTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATgc                                                                                                                           <  3:328964/149‑1 (MQ=255)
        tAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGTGAGCTGGATGCGCGt                                                                                                                       >  2:266890/1‑149 (MQ=255)
        tAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGt                                                                                                                       >  6:163562/1‑149 (MQ=255)
         aGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCGCGAGGAGATCCGCGAGCTGGATGcgcggg                                                                                                                      >  6:127938/1‑147 (MQ=255)
           aTTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAg                                                                                                                    <  4:311538/149‑1 (MQ=255)
           aTTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAg                                                                                                                    >  2:40659/1‑149 (MQ=255)
             ttCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAg                                                                                                                    >  2:346126/1‑147 (MQ=255)
                     gTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAg                                                                                                           >  1:104135/1‑148 (MQ=255)
                     gTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCAGCGAGCTGGATGCGCGTGAGCAGGCCCAg                                                                                                           >  7:148570/1‑148 (MQ=255)
                       tGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCa                                                                                                         >  1:325897/1‑148 (MQ=255)
                         aaaCACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCaaa                                                                                                       <  5:163562/148‑1 (MQ=255)
                              ccGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGTGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGt                                                                                                 <  1:266890/149‑1 (MQ=255)
                              ccGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGt                                                                                                 <  4:292935/149‑1 (MQ=255)
                                ggCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTgg                                                                                               <  1:346126/149‑1 (MQ=255)
                                 gCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTgg                                                                                               >  3:382886/1‑148 (MQ=255)
                                 gCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTgg                                                                                               <  4:382886/148‑1 (MQ=255)
                                    cGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTg                                                                                                <  6:72363/144‑1 (MQ=255)
                                    cGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTg                                                                                                >  5:72363/1‑144 (MQ=255)
                                     gTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGCGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGtt                                                                                          >  8:353917/1‑149 (MQ=255)
                                     gTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGtt                                                                                          <  2:310301/149‑1 (MQ=255)
                                      tttGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTc                                                                                         <  7:171728/149‑1 (MQ=255)
                                         ggCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGa                                                                                       <  1:40659/148‑1 (MQ=255)
                                         ggCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGa                                                                                       <  7:167861/148‑1 (MQ=255)
                                                aaaGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATgccg                                                                               >  2:225355/1‑149 (MQ=255)
                                                 aaGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATgccgc                                                                              <  1:225355/149‑1 (MQ=255)
                                                  aGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATgccgcc                                                                             <  5:130525/149‑1 (MQ=255)
                                                            gatgatTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAgg                                                                   <  7:281904/149‑1 (MQ=255)
                                                            gatgatTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAgg                                                                   <  2:325897/149‑1 (MQ=255)
                                                                atTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTa                                                                <  1:322477/148‑1 (MQ=255)
                                                                  tGCAAAAGAAGCCCGTCAGGCCCTGGCGGTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTAtt                                                              >  1:223788/1‑148 (MQ=255)
                                                                    cAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGa                                                            <  7:353917/148‑1 (MQ=255)
                                                                     aaaaGAAGCCCGTCAGGCCCTGGCGCTTCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCTTGAGCAggaa                                                                                                              <  3:67356/97‑3 (MQ=255)
                                                                     aaaaGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGcc                                                                                                              >  4:67356/1‑97 (MQ=255)
                                                                     aaaaGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGc                                                          >  3:270509/1‑149 (MQ=255)
                                                                           agCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTggg                                                    <  1:177753/149‑1 (MQ=255)
                                                                            gCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTgg                                                     >  7:211535/1‑147 (MQ=255)
                                                                                gTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCCAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTg                                               <  2:104135/149‑1 (MQ=255)
                                                                                      cccTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTg                                         <  2:17784/149‑1 (MQ=255)
                                                                                       ccTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGa                                                                           <  8:265331/114‑1 (MQ=255)
                                                                                       ccTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGa                                                                           >  7:265331/1‑114 (MQ=255)
                                                                                           gcgcTGCCTCGGTTTGCCGAATACAATCAGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGTGCGTTCGAGGATGCCGCCGAGGAACAGTGTATTGAGCCTTGGGAGCTGGCGCTGCaa                                      <  8:328344/147‑1 (MQ=255)
                                                                                            cgcTGCAGCGCTTTGACGAAGCCAATCCGCAGTTGCTCGATTAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGTCGTTCGAGGATGCCTCCGAGGAACAGGGTATTGAGCCTTGTGAGCTGGCTCTGCAattg                                   <  8:211535/149‑1 (MQ=255)
                                                                                                           ccGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCg                    >  3:113104/1‑149 (MQ=255)
                                                                                                            cgcaGCCACTCCGCAGTTGCTCGAGGAGCTGCGCGAGCTGGATGGGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCgg                   <  4:301254/146‑1 (MQ=255)
                                                                                                                   aaTCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCt            >  7:147456/1‑149 (MQ=255)
                                                                                                                          aGTTGCTCGTGTACATCCGAGAGCTGGATGCGCGTGAGCAGGCCCAGCTAATCGTGTGGGCGTTCGAGGATGCCTCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCTGAGCAGCTGCAgg       <  8:148570/147‑1 (MQ=255)
                                                                                                                            ttGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCAt   >  3:274019/1‑149 (MQ=255)
                                                                                                                            cagCTCGCGGATCTCCTCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTCGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAattgatt                                >  7:180753/3‑120 (MQ=255)
                                                                                                                            cagCTCGCGGATCTCCTCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTCGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAattgatt                                <  8:180753/118‑1 (MQ=255)
                                                                                                                             agCTCGCGGATCTCCTCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCATg  >  8:47375/2‑149 (MQ=255)
                                                                                                                                       gaTCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCAt   <  2:74899/138‑1 (MQ=255)
                                                                                                                                       gaTCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCAt   >  1:74899/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                          
AAAGCTGTTAGATTTCAGCCGGTTGAAACACCGGCCCGTTTGGCCGGCAAAGGACACCACGATGATTGCAAAAGAAGCCCGTCAGGCCCTGGCGCTGCAGCGCTTTGCCGAAGCCAATCCGCAGTTGCTCGAGGAGATCCGCGAGCTGGATGCGCGTGAGCAGGCCCAGCAAATCGAGTGGGCGTTCGAGGATGCCGCCGAGGAACAGGGTATTGAGCCTTGGGAGCTGGCGCTGCAATTGATTGCCGAATCCCCGGAGCAGCTGCAGGCCATG  >  NC_002947/2626325‑2626598

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