Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 374,775 T→G 8.5% L22R (CTG→CGG)  PP_0312 → electron transfer flavoprotein subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947374,7750TG 8.5% 135.7 / 5.2 47L22R (CTG→CGG) PP_0312electron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base G (2/2);  total (24/23)
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.12e-01

TGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCACCGGCG  >  NC_002947/374626‑374921
                                                                                                                                                     |                                                                                                                                                   
tggtggCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGACGAGTGGATCGCCCGTAACCGTCTGCACCCGCt                                                                                                                                                     >  7:340206/1‑149 (MQ=255)
   tggCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGCGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGc                                                                                                                                                   <  3:188760/149‑1 (MQ=255)
     gCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCAt                                                                                                                                                 <  8:231333/149‑1 (MQ=255)
        gACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGcc                                                                                                                                              <  4:49855/149‑1 (MQ=255)
             gCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGAt                                                                                                                                        >  1:21154/1‑149 (MQ=255)
                 aTCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATg                                                                                                                                       >  7:180866/1‑146 (MQ=255)
                     aggaggACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCATCGCGGCATGCAGCGGGTGCAGACgg                                                                                                                                <  4:228590/149‑2 (MQ=255)
                             gCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCcgc                                                                                                                                            <  1:251144/130‑1 (MQ=255)
                               aTCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAg                                                                                                                      >  6:141254/1‑149 (MQ=255)
                                tCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGc                                                                                                                     >  7:91333/1‑149 (MQ=255)
                                  cACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCt                                                                                                                    <  8:91333/148‑1 (MQ=255)
                                   aCTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTg                                                                                                                   >  3:88321/1‑148 (MQ=255)
                                   aCTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCATCGCGGCATGCAGCGGGTGCAGACggttac                                                                                                                            <  3:42194/139‑6 (MQ=37)
                                   aCTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCATCGCGGCATGCAGCGGGTGCAGACggttac                                                                                                                            >  4:42194/1‑134 (MQ=37)
                                     tGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGAt                                                                                                                <  8:180866/149‑1 (MQ=255)
                                      gTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATg                                                                                                               >  4:8843/1‑149 (MQ=255)
                                       tCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATgg                                                                                                              <  2:21154/149‑1 (MQ=255)
                                          cAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGAGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGc                                                                                                           >  4:278205/1‑149 (MQ=255)
                                           aaaGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATgggg                                                                                                            >  4:100751/1‑147 (MQ=255)
                                           aaaGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATgggg                                                                                                            <  3:100751/147‑1 (MQ=255)
                                                  cAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGc                                                                                                   <  6:69797/149‑1 (MQ=255)
                                                        gCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATcc                                                                                            >  4:190459/1‑149 (MQ=255)
                                                          tAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGCCCAGCTGGATGGGGCCCAACGGCATCATCCg                                                                                           >  6:93197/1‑149 (MQ=255)
                                                            aaCGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGc                                                                                          <  4:88321/148‑1 (MQ=255)
                                                                      tGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCcaca                                                                               >  6:359992/1‑149 (MQ=255)
                                                                        cGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                                                            <  3:190459/149‑1 (MQ=255)
                                                                          gAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCg                                                                           <  5:141254/149‑1 (MQ=255)
                                                                                            aCATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCCCATGCGATTGCCGCAGGCTTTGTc                                                         >  8:261968/1‑149 (MQ=255)
                                                                                                tATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcggcc                                                     >  4:124393/1‑149 (MQ=255)
                                                                                                  tCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCCTCATCCGCAAGACCCCACATGCGATTGCCGCAGGCTTTGTcggcccg                                                   >  4:38081/1‑149 (MQ=255)
                                                                                                    cgccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcggcccggc                                                 >  8:306955/1‑149 (MQ=255)
                                                                                                      ccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCccggccg                                               <  4:350297/149‑1 (MQ=255)
                                                                                                              cACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAg                                       <  8:340206/149‑1 (MQ=255)
                                                                                                               aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgc                                      <  3:278205/149‑1 (MQ=255)
                                                                                                               aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgc                                      <  4:176553/149‑1 (MQ=255)
                                                                                                               aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgc                                      >  3:176553/1‑149 (MQ=255)
                                                                                                                         tGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCg                            >  1:183931/1‑149 (MQ=255)
                                                                                                                           gATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCa                          <  8:83535/149‑1 (MQ=255)
                                                                                                                                 ccGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAgcggcg                    >  4:18581/1‑149 (MQ=255)
                                                                                                                                        cGTCTGCATCGCGGCATGCAGCGGGTGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCcag               >  7:193614/1‑147 (MQ=255)
                                                                                                                                             gCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGgcgc                   >  3:32948/1‑138 (MQ=255)
                                                                                                                                             gCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGgcgc                   <  4:32948/138‑1 (MQ=255)
                                                                                                                                              cACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                                                            >  6:108847/1‑80 (MQ=255)
                                                                                                                                              cACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                                                            <  5:108847/80‑1 (MQ=255)
                                                                                                                                              cACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCAGCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCAc       <  3:38081/149‑1 (MQ=255)
                                                                                                                                               aCCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCAcc      >  3:286049/1‑149 (MQ=255)
                                                                                                                                                 ccGCTGCATGCCGCGATGAAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCAGCAAGAACCCACATGCGATTGCCCCAGGCTTTGTCGGCCCGGCCGGCATCAAGCCCATCGAACGCAGCGGCCCCCAGCACGGCACCgg    >  7:140028/1‑149 (MQ=255)
                                                                                                                                                   gCTGCATGCCGCGATGCAGACGCAACAGACCAGATGGATGGGGCCCATCGGCATCATCCGCAAGAACCCTCATGCGTTTGCCGCAGGCTTTGTCGGCCCGGGCGGCATCAAGCGCATCGACCGCAGAGGCGCCCAGCAGGGCACCTgcg  <  8:140028/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGGTGGCCGACGTGCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGCGATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCACCGGCG  >  NC_002947/374626‑374921

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