Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,337,599 A→G 12.1% A103A (GCA→GCG)  PP_1163 → FAD‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,337,5990AG12.1% 79.1 / 3.5 33A103A (GCA→GCG) PP_1163FAD‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (18/11);  new base G (2/2);  total (20/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTAGAGCAGGCGTTTGGCAACTGCCCGCAGTACATCCAGTTGC  >  NC_002947/1337476‑1337747
                                                                                                                           |                                                                                                                                                    
ccAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTggg                                                                                                                            >  1:276255/1‑150 (MQ=255)
ccAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTggg                                                                                                                            >  1:468541/1‑150 (MQ=255)
 cAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTggg                                                                                                                            >  1:308670/1‑149 (MQ=255)
 cAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTggg                                                                                                                            <  2:308670/149‑1 (MQ=255)
       cGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGaa                                                                                                                     <  2:591943/150‑1 (MQ=255)
         gCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAATCTCGCCTGACCATCGTATGCTGACGATCCACAGCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGGCGGCCTGCTGGTGATCGAACt                                                                                                                   >  1:152912/1‑150 (MQ=255)
         gCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACGCTCGTATGCTGACGATCCACACCACGTGGCCAGCAGATGACCCCGCCACCCCCGGATTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACt                                                                                                                   >  1:593037/1‑150 (MQ=255)
         gCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACt                                                                                                                   >  2:331494/1‑150 (MQ=255)
         gCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGGCGGCCTGCTGGGGATCGAACt                                                                                                                   >  1:152950/1‑150 (MQ=255)
            ttGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGca                                                                                                                >  2:148491/1‑150 (MQ=255)
               ggCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCacac                                                                                                             >  2:603158/1‑150 (MQ=255)
                                   cAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCa                                                                                         >  2:179220/1‑150 (MQ=255)
                                     gAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTggg                                                                                                                            <  1:112212/113‑1 (MQ=255)
                                     gAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTggg                                                                                                                            >  2:112212/1‑113 (MQ=255)
                                                 aGCTCGCCTGCCCCTCGTATGCTGCCGCTACACACCCCGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGc                                                                           <  2:468541/150‑1 (MQ=255)
                                                  gCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGcc                                                                          >  2:327291/1‑150 (MQ=255)
                                                      gCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGGCCTTCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGc                                                                      >  1:375259/1‑150 (MQ=255)
                                                       ccTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACg                                                                                      >  2:98777/1‑133 (MQ=255)
                                                       ccTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACg                                                                                      <  1:98777/133‑1 (MQ=255)
                                                                cGTATGCTGACTATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCACCAACGGGCAGATCCGCCACGCTGCCCAAggg                                                            <  2:276255/150‑1 (MQ=255)
                                                                       tGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGc                                                     >  1:525917/1‑150 (MQ=255)
                                                                           gATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGt                                                 >  1:610163/1‑150 (MQ=255)
                                                                               cacaCCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGcc                                             <  2:610163/150‑1 (MQ=255)
                                                                                  accacGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTa                                          >  1:137549/1‑150 (MQ=255)
                                                                                    cacGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTaga                                        <  1:603158/150‑1 (MQ=255)
                                                                                                    ccggctgCCAACCCGGGGGTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTAGAGCAGGCGTTTGGCAAc                        >  2:1694/7‑150 (MQ=255)
                                                                                                       gctgCCAACCCGGGGGTGGCGGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAgg                                                             <  2:203492/107‑1 (MQ=255)
                                                                                                       gctgCCAACCCGGGGGTGGCGGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAgg                                                             <  2:203543/107‑1 (MQ=255)
                                                                                                       gctgCCAACCCGGGGGTGGCGGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAgg                                                             >  1:203543/4‑110 (MQ=255)
                                                                                                       gctgCCAACCCGGGGGTGGCGGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAgg                                                             >  1:203492/4‑110 (MQ=255)
                                                                                                       gctgCCAACCCGGGGGTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTAGAGCAGGCGTTTGGCAACTGc                     <  1:1694/147‑1 (MQ=255)
                                                                                                          gCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTAGAGCAGGCGTTTGGCAACTGCCCg                  <  2:525917/150‑1 (MQ=255)
                                                                                                                          cAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTAGAGCAGGCGTTTGGCAACTGCCCGCAGTACATCCAGTTGc  <  2:137549/150‑1 (MQ=255)
                                                                                                                           |                                                                                                                                                    
CCAGGGTCGGCCTTGGGCCACGCTGCTGGAAGGGCCAGAAGGGTTCGTAAGCTCGCCTGACCCTCGTATGCTGACGATCCACACCACGTTGCCAGCAGATGACCCCGCCACCCCCGGGTTGGCAGCCGGGCAGGCGGTCGGCCTGCTGGGGATCGAACTGCACACCCGTCGGCGCAATCGCATCAACGGGCAGATCCGCCACGCTGCCCAAGGGCAACTGCAAGTGGCCGTAGAGCAGGCGTTTGGCAACTGCCCGCAGTACATCCAGTTGC  >  NC_002947/1337476‑1337747

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