Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 408,719 T→C 16.6% L491L (TTA→TTG metE ← 5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine S‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730408,7190TC16.6% 63.8 / 5.5 30L491L (TTA→TTGmetE5‑methyltetrahydropteroyltriglutamate‑‑ homocysteine S‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (15/10);  new base C (2/3);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CTGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGTTTCATA  >  CP000730/408592‑408856
                                                                                                                               |                                                                                                                                         
cTGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTa                                                                                                                              >  1:187079/1‑141 (MQ=255)
 tGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCCCCCAACCGAATTTAGTTACTAAGAAACCTTGTaa                                                                                                                             >  2:30855/1‑141 (MQ=255)
    tAAGCTTTGTGCATAAACTGGTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAAttt                                                                                                                          <  1:43470/141‑1 (MQ=255)
         tttGTGCATAAACTGTTTCATCAACAGGTAAAGGCGCTGTCCATTTTACATCACCATAATTGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGACACCTTGTAATTTTactc                                                                                                                     >  2:191787/1‑141 (MQ=255)
            gtgCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCga                                                                                                                  >  1:44473/1‑141 (MQ=255)
              gCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCgaag                                                                                                                >  2:225929/1‑141 (MQ=255)
              gCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCgaag                                                                                                                >  2:247258/1‑141 (MQ=255)
              gCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCgaag                                                                                                                >  2:247245/1‑141 (MQ=255)
                aTAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCgaagaa                                                                                                              <  2:280813/141‑1 (MQ=255)
                    aCTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCa                                                                                                          <  1:286711/141‑1 (MQ=255)
                      tGTTTCATCAACAGTTAAAGGCGCTGTCCCTTTTACATCGCCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAAc                                                                                                        <  1:194012/141‑1 (MQ=255)
                       gTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAAcc                                                                                                       >  1:38811/1‑141 (MQ=255)
                                                             aCCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATc                                                                 >  2:256809/1‑141 (MQ=255)
                                                                   aaTGATTGGTTGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGcc                                                           <  1:50072/141‑1 (MQ=255)
                                                                   aaTGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAGCCATGTCATTACGTTCAAATTCCCCGTGAACTAATACATCTACGcc                                                           >  1:86434/1‑141 (MQ=255)
                                                                    aTGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAGATTCACCGTGAACTAATACATCTAAGCCa                                                          >  1:220927/1‑141 (MQ=255)
                                                                          ggtggtTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTTTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCTTGGACTAATACCTCTACGCCAACGtct                                                    >  2:206665/1‑141 (MQ=255)
                                                                                    cGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTTAACTAATACATCTAAGCCAATGTCTTCTTGAAttt                                          <  2:246130/141‑1 (MQ=255)
                                                                                        ccGTGAACCATATGATTGCACCCAACCCAACTCAGTCACCAACAAGCCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAAt                                      <  1:52747/140‑1 (MQ=255)
                                                                                          gTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATcc                                    >  1:101229/1‑141 (MQ=255)
                                                                                                         gCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAATTTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAt                     >  2:102437/1‑141 (MQ=255)
                                                                                                            cccAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttt                  >  2:159611/1‑141 (MQ=255)
                                                                                                              cAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttttt                <  1:225929/141‑1 (MQ=255)
                                                                                                              cAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttttt                <  1:247245/141‑1 (MQ=255)
                                                                                                              cAACCGAATTTAGTTACCAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCAtttttt                <  1:247258/141‑1 (MQ=255)
                                                                                                               aaCCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTa               <  1:188577/141‑1 (MQ=255)
                                                                                                                aCCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTaa              <  2:38811/141‑1 (MQ=255)
                                                                                                                    aaTTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTCCGTTCAAATTCACCGAGAACTACTACATCTAAGCCAATGTCATCTTGAAATTTAATCCATCGAGCAATATCAtttatgcagact          >  1:273818/1‑132 (MQ=255)
                                                                                                                      tttAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGt        >  2:123121/1‑141 (MQ=255)
                                                                                                                            tACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGTTTCata  <  1:71982/141‑1 (MQ=255)
                                                                                                                               |                                                                                                                                         
CTGTTAAGCTTTGTGCATAAACTGTTTCATCAACAGTTAAAGGCGCTGTCCATTTTACATCACCATAAATGATTGGTGGTTTTACGGCACGTGAACCATATGATTGCACCCAACCGAATTTAGTTACTAAGAAACCTTGTAATTTTTCTCCGAAGAATTCAACCATGTCATTACGTTCAAATTCACCGTGAACTAATACATCTAAGCCAATGTCTTCTTGAATTTTAATCCATCGAGCAATTTCATTTTTTAAGAATGTTTCATA  >  CP000730/408592‑408856

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