Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,469,323 T→A 8.4% G282G (GGT→GGA)  ydbC → putative NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,469,3230TA 8.4% 126.2 / 3.2 48G282G (GGT→GGA) ydbCputative NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base A (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

AATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTACCAGTA  >  CP009273/1469179‑1469466
                                                                                                                                                |                                                                                                                                               
aaTGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACGGAATATTGTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTAttt                                                                                                                                             >  2:639374/1‑149 (MQ=255)
 aTGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATCGCATCTTTCGGAGGAAGTGTTGTCTACGTTGGATGGTATTTc                                                                                                                                            >  2:21632/1‑149 (MQ=255)
                    ggCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCGACGTTGGATGGTATTTCGCGAGAATAACGAATATAc                                                                                                                         >  2:671545/1‑149 (MQ=255)
                       tGTTACAGCGTTCACCGAATATTTTGGTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACaaa                                                                                                                      <  1:625182/149‑1 (MQ=255)
                       ggTTACCGCGTTCACCGAAGATTTTGCTGATCCCAGGGACGTCTTCGGGTGGGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTGCGCTAGAATAACGACTATACaaa                                                                                                                      <  1:621416/148‑1 (MQ=255)
                        gTTACAGCGTCCCCCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATAAATACAAAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACaaaa                                                                                                                     <  1:513054/149‑1 (MQ=255)
                          tACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGGATATACAAAAgg                                                                                                                   >  2:840817/1‑149 (MQ=255)
                          tACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAgg                                                                                                                   <  1:440379/149‑1 (MQ=255)
                          tACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAgg                                                                                                                   <  2:809956/149‑1 (MQ=255)
                             aGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGaa                                                                                                                <  2:570895/149‑1 (MQ=255)
                                 ttCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTCCGTTGGATGGTATTTCGCGAGAATAACGGAGATACAAAAGGGAAAGAt                                                                                                            >  1:7815/1‑149 (MQ=255)
                                 ttCACCGAATATTTTGCTGATCCCAGGGACGTCTACGGTTGCGCATTTTCGGGAGAATATGGCAGATGAATGATAGTGTCTTTCAGAAGAATTGCTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGAt                                                                                                            <  2:854264/149‑1 (MQ=255)
                                  tCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAACTTGCATCTTTCTGAGGAAGTGTTGTCTACGTAGGATGGTATTTCGCGAGAATAACGAATATACAAAACGGAAAGATg                                                                                                           >  1:248906/1‑149 (MQ=255)
                                    aCCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCa                                                                                                         >  2:707657/1‑149 (MQ=255)
                                                          gggACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTaa                                                                                   <  2:58895/149‑1 (MQ=255)
                                                           ggACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGAATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGGATTTCCCTTTTTTTCTTTTTTAAt                                                                                  >  2:5937/1‑149 (MQ=255)
                                                           ggACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGAATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAAt                                                                                  >  2:829427/1‑149 (MQ=255)
                                                               ggcTTGGGCTGCGCATTTACGTGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCa                                                                              <  1:64098/147‑1 (MQ=255)
                                                                tCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGGTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAGTGGcgt                                                                             >  2:802760/1‑147 (MQ=255)
                                                                  ttCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCGTTTCCCTTTTTTTCTTTTTTAATGGCATgg                                                                           >  2:65808/1‑149 (MQ=255)
                                                                         gcgcATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTGCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCat                                                                    <  1:323846/149‑1 (MQ=255)
                                                                               ttACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTg                                                              <  2:829361/149‑1 (MQ=255)
                                                                                    ggAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGc                                                         >  1:344096/1‑149 (MQ=255)
                                                                                     gagaATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATCTGTTGAAGGcc                                                        >  2:434298/1‑149 (MQ=255)
                                                                                        aaTATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAgg                                                     >  2:542810/1‑149 (MQ=255)
                                                                                             gggttctgAAAAATTGCGTCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTg                                                <  1:205528/144‑1 (MQ=255)
                                                                                              gctgcgGAAAAATGGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGc                                               <  2:347225/149‑1 (MQ=255)
                                                                                                tgctgAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATTTTGAATGCCCGGAGTTACgg                                             >  1:460981/1‑149 (MQ=255)
                                                                                                  ctgAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATAGCATGGATTGCAGATGTTGACCGCAATGAGTGGCGGCt                                           >  1:19376/1‑149 (MQ=255)
                                                                                                  ctgAAAAATTGCATCCTACTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTTCCTTTTTTTCTTTTTGAATGGCATTGAGTGCATCTGTTGTAGGCCAGTGATTCCGGCt                                           >  1:251214/1‑149 (MQ=255)
                                                                                                    aaacaaTTGTATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACCAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTgg                                         <  1:616101/145‑1 (MQ=255)
                                                                                                       aaaTTGCATCTTTCTGAGGAAGTGTTGCCTACGTTGGATGGTCTTTCGCGAGAATAACGCATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCGTGAAGTGCACATGTTGAAGGCCGGGAATTGCGGCTGGGTa                                      >  1:268304/1‑149 (MQ=255)
                                                                                                       aaaTTGCATCTTTCTGAGGAAGTGTTGCCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATa                                      >  1:707414/1‑149 (MQ=255)
                                                                                                        aaTTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATAc                                     >  1:863227/1‑149 (MQ=255)
                                                                                                           tGCATCTTTTCGAGGAAGTGGTGTCGACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGtt                                  <  2:104213/149‑1 (MQ=255)
                                                                                                              atttttCTTAGCAAGTCTTGTCAACGTTGGATGGTAGTTCGCGAGAATAACGAATATACAAAAGGGAAATAGGCATTTCCCTTTTTTTCTTTTTTAATGGCATGAAGTGCATATGTTGAAGGCCAGGAATTGCGGCTTGATACGTTTAc                               <  1:724355/146‑1 (MQ=255)
                                                                                                              aTCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTAc                               <  1:844217/149‑1 (MQ=255)
                                                                                                                 tttCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGc                            >  2:232348/1‑149 (MQ=255)
                                                                                                                   tCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGACGGCATATCTTGCATGCCAGGAATTGCGGCTGGCTACGTGTACTcccc                          >  1:773967/1‑149 (MQ=255)
                                                                                                                           aGTGTTGTCTACGTTGGATGGAATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTAt                  <  1:829427/149‑1 (MQ=255)
                                                                                                                           aGTGTTGTCTACGTTGGATGGAATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTAt                  <  1:5937/149‑1 (MQ=255)
                                                                                                                                tggcTACTTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGa             <  2:217592/146‑1 (MQ=255)
                                                                                                                                     aCGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTa        >  1:512878/1‑149 (MQ=255)
                                                                                                                                     aCGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTa        >  2:265461/1‑149 (MQ=255)
                                                                                                                                      cggtGGATGGTCTTGCGCGGGAATAGCGGACATACAAAAGGGACAGGTGCCTTTCCCTCTTTTTCTTGTTTAAGGGCATGGCGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTAc       <  1:752825/146‑1 (MQ=255)
                                                                                                                                      cGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGGGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTAc       <  2:79610/149‑1 (MQ=255)
                                                                                                                                          ggATGTGATTTCGCGAGAATAACGAATATTCAAAAGGGAAAGATGCATTTCCCTTTTTTACTTTTGTAATGGCGTGGAGTGCATATCGTGACGGCCAGGAGTTGCGGCGGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTACCAGt   <  1:637534/149‑1 (MQ=255)
                                                                                                                                           gATGGTAGTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTACCAGTa  <  2:473876/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGAGAATATGGCTGCTGAAAAATTGCATCTTTCTGAGGAAGTGTTGTCTACGTTGGATGGTATTTCGCGAGAATAACGAATATACAAAAGGGAAAGATGCATTTCCCTTTTTTTCTTTTTTAATGGCATGGAGTGCATATGTTGAAGGCCAGGAATTGCGGCTGGATACGTTTACTGCCCCTGTTTATGCTGAGTCTACCAGTA  >  CP009273/1469179‑1469466

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