Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 269,112 A→G 10.5% L281P (CTG→CCG)  desA ← delta‑9 fatty acid desaturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947269,1120AG10.5% 92.9 / 3.3 38L281P (CTG→CCG) desAdelta‑9 fatty acid desaturase
Reads supporting (aligned to +/- strand):  ref base A (15/19);  new base G (2/2);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

GACTCGATGCGCACGTGGTGGCGGTCTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTGGCCAGACGCAGCAGGCAGAACAGGCGGATC  >  NC_002947/268977‑269261
                                                                                                                                       |                                                                                                                                                     
gACTCGATGCGCACGTGGTGGCGGTCTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTgg                                                                                                                                         >  1:503721/1‑150 (MQ=255)
       tGCGCACGTGGTGGCGGTCTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCGCGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGa                                                                                                                                  <  2:829130/150‑1 (MQ=255)
          gcACGTGGTGGCGGTCTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCt                                                                                                                               <  1:36162/150‑1 (MQ=255)
               tggtggCGGTCTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAg                                                                                                                          >  2:452965/1‑150 (MQ=255)
                ggtggCGGTCTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGc                                                                                                                         <  2:48247/150‑1 (MQ=255)
                         cTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGTGTGCGGTACTGGGCCATGATCTGGAAGCGGttgttgc                                                                                                                >  1:262651/1‑149 (MQ=255)
                            ccAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGa                                                                                                             >  1:796767/1‑150 (MQ=255)
                              agcagGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTgg                                                                                                                                         <  2:52452/120‑1 (MQ=255)
                              agcagGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTgg                                                                                                                                         >  1:52452/1‑120 (MQ=255)
                              agcagGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAg                                                                                                               >  2:311354/1‑146 (MQ=255)
                              agcagGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCAGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAg                                                                                                               <  1:311354/146‑1 (MQ=255)
                                     tGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCtt                                                                                                                                                   >  2:153586/1‑103 (MQ=255)
                                     tGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCtt                                                                                                                                                   <  1:153586/103‑1 (MQ=255)
                                     tGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCg                                                                                                                                                <  2:796878/106‑1 (MQ=255)
                                     tGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCg                                                                                                                                                >  1:796878/1‑106 (MQ=255)
                                      ggTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGc                                                                                                   <  2:579560/150‑1 (MQ=255)
                                              gcgAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGtgtcca                                                                                           >  1:33139/1‑150 (MQ=255)
                                                       gcCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCGTGGCGGTGTCCATGTCCAGGc                                                                                  >  2:27995/1‑150 (MQ=255)
                                                         cGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTTGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTg                                                                                <  1:772998/150‑1 (MQ=255)
                                                                cACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTCGCCTTGc                                                                         <  2:33139/150‑1 (MQ=255)
                                                                  cggcggAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCg                                                                       <  1:471127/150‑1 (MQ=255)
                                                                   ggcggAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCgg                                                                      >  2:828873/1‑150 (MQ=255)
                                                                        aaGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCAcc                                                                 >  1:300667/1‑150 (MQ=255)
                                                                             gTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTg                                                            <  2:196795/150‑1 (MQ=255)
                                                                              tGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTgg                                                           <  2:631757/150‑1 (MQ=255)
                                                                                     ccGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATgg                                                    <  2:796767/150‑1 (MQ=255)
                                                                                      cGAGGCGTCGACACGGGCCAGTTCCTGCTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATggg                                                   <  1:817304/150‑1 (MQ=255)
                                                                                          gCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGcc                                               >  1:131989/1‑150 (MQ=255)
                                                                                                  acGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGTATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTg                                       <  1:828873/150‑1 (MQ=255)
                                                                                                    gggCCAGTTCCTGTTTCACCAGCGGGCCGATCACCGGCTTGCGGTACTGGGCCATGATCTgg                                                                                                                             >  1:680207/1‑62 (MQ=255)
                                                                                                    gggCCAGTTCCTGTTTCACCAGCGGGCCGATCACCGGCTTGCGGTACTGGGCCATGATCTgg                                                                                                                             <  2:680207/62‑1 (MQ=255)
                                                                                                    gggCCAGTTCCTGTTTCACCAGCGGGCCGATCACCGGCTTGCGGTACTGGGCCATGATCTgg                                                                                                                             <  2:679437/62‑1 (MQ=255)
                                                                                                    gggCCAGTTCCTGTTTCACCAGCGGGCCGATCACCGGCTTGCGGTACTGGGCCATGATCTgg                                                                                                                             >  1:679437/1‑62 (MQ=255)
                                                                                                         aGTTCCTGTTTCGCCCGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTg                                <  2:262651/150‑1 (MQ=255)
                                                                                                         aGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTg                                <  1:452965/150‑1 (MQ=255)
                                                                                                                           gggCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTGGCCAGACGCAgcaggcag              >  2:204885/1‑150 (MQ=255)
                                                                                                                                 aTCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTGGCCAGACGCAgcagg                 >  1:193713/1‑141 (MQ=255)
                                                                                                                                 aTCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTGGCCAGACGCAgcagg                 <  2:193713/141‑1 (MQ=255)
                                                                                                                                       aGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTGGCCAGACGCAGCAGGCAGAACAGGCGGATc  >  2:567499/1‑150 (MQ=255)
                                                                                                                                       |                                                                                                                                                     
GACTCGATGCGCACGTGGTGGCGGTCTTCCAGCAGGCTGGTTTCGCGCGAGAGCAGCCGTTTGGCACGGCGGAAGCGGTGGCGCACCGAGGCGTCGACACGGGCCAGTTCCTGTTTCACCAGCGGGCCGATCACCAGCTTGCGGTACTGGGCCATGATCTGGAAGCGGTTGTTGAGGATGGCCATGGCGGTGTCCATGTCCAGGCTGGCCTTGCCGGCCACCCGGTGGGCGATGGGCGCCACGCGCTGCACCTTGGCCAGACGCAGCAGGCAGAACAGGCGGATC  >  NC_002947/268977‑269261

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