Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,768,022 G→A 9.3% R222C (CGT→TGT)  ydiJ ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,768,0220GA 9.3% 127.7 / 3.4 43R222C (CGT→TGT) ydiJputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base A (2/2);  total (22/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 = 9.95e-01

ATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTT  >  NC_000913/1767874‑1768170
                                                                                                                                                    |                                                                                                                                                    
ataaAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACg                                                                                                                                                      >  3:79197/1‑149 (MQ=255)
         aGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCCGTCAGCAACgcca                                                                                                                                             <  4:243849/149‑4 (MQ=255)
          gCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGa                                                                                                                                            >  3:157943/1‑149 (MQ=255)
           cGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGAt                                                                                                                                           >  1:439290/1‑149 (MQ=255)
               ccTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGtt                                                                                                                                              >  1:393548/1‑142 (MQ=255)
               ccTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGtt                                                                                                                                              <  2:393548/142‑1 (MQ=255)
                  tCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCgt                                                                                                                                    >  2:450815/1‑149 (MQ=255)
                    tGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCgtgt                                                                                                                                  >  4:470162/1‑149 (MQ=255)
                     gAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGtt                                                                                                                                 <  2:57040/149‑1 (MQ=255)
                      aaCCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTta                                                                                                                                >  2:416773/1‑149 (MQ=255)
                       aCCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTtat                                                                                                                               <  1:241627/149‑1 (MQ=255)
                         ccatCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATaa                                                                                                                             >  4:346919/4‑149 (MQ=255)
                                   gcgcGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGa                                                                                                                   <  4:14761/149‑1 (MQ=255)
                                                     ggTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTAc                                                                                                 <  4:157943/149‑1 (MQ=255)
                                                       tcatctcatcGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCgtg                                                                                                                                   <  4:247726/113‑1 (MQ=255)
                                                       tcatctcatcGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCgtg                                                                                                                                   >  3:247726/1‑113 (MQ=255)
                                                          tctcatcGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGc                                                                                            >  3:39724/1‑149 (MQ=255)
                                                           ctcatcGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCg                                                                                           <  1:416773/149‑1 (MQ=255)
                                                              atcGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATaaa                                                                                                                            <  2:308080/113‑1 (MQ=255)
                                                              atcGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATaaa                                                                                                                            >  1:308080/1‑113 (MQ=255)
                                                                   tAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAA‑CATGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCa                                                                                  <  4:79197/149‑1 (MQ=255)
                                                                                   aTCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAg                                                                   <  4:166725/149‑1 (MQ=255)
                                                                                     cGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGt                                                                 <  2:439290/149‑1 (MQ=255)
                                                                                       tAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTAGTATTGGATTTACCCAGCGTTTCCGCCAATTAGACGGGTAAAGGTTg                                                               >  4:331088/1‑149 (MQ=255)
                                                                                        aaCCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATCACCCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGc                                                              <  2:163409/149‑1 (MQ=255)
                                                                                        aaCCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGc                                                              <  4:327519/149‑1 (MQ=255)
                                                                                            ggTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTAt                                                          >  2:130174/1‑149 (MQ=255)
                                                                                             gTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATc                                                         >  2:22382/1‑149 (MQ=255)
                                                                                               aaGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCga                                                       >  4:261810/1‑149 (MQ=255)
                                                                                                    aGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAata                                                  >  4:228025/1‑149 (MQ=255)
                                                                                                         ttAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGcc                                             >  2:235554/1‑149 (MQ=255)
                                                                                                          tAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGccc                                            <  2:361462/149‑1 (MQ=255)
                                                                                                           aaGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACAGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGcccc                                           >  4:61188/1‑149 (MQ=255)
                                                                                                             gTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGcccccc                                         <  4:39724/149‑1 (MQ=255)
                                                                                                                tGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCcaac                                      <  1:130174/149‑1 (MQ=255)
                                                                                                                 ggggAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCcaaca                                     >  3:63182/1‑149 (MQ=255)
                                                                                                                 ggggAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTAACCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCcaaca                                     >  3:195477/1‑149 (MQ=255)
                                                                                                                    gAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCACCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAacaaca                                  <  2:99559/149‑1 (MQ=255)
                                                                                                                                 aTTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGcc                     <  4:219189/149‑1 (MQ=255)
                                                                                                                                    aaCTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAg                  <  3:261810/149‑1 (MQ=255)
                                                                                                                                        aacGTTGCCGTCAGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTa                <  2:51750/145‑1 (MQ=255)
                                                                                                                                        aacGTTGCCGTCAGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTa                >  1:51750/3‑147 (MQ=255)
                                                                                                                                         gCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGTTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGAGTACGCCAAGTACGt             >  2:138699/1‑149 (MQ=255)
                                                                                                                                               cTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTg       >  2:81024/1‑149 (MQ=255)
                                                                                                                                                    ggCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGtt  >  3:128615/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTT  >  NC_000913/1767874‑1768170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑