Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,569,373 A→G intergenic (‑209/+84) eutS ← / ← maeB ethanolamine utilization microcompartment protein EutS/NADP‑dependent oxaloacetate‑decarboxylating malate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,569,3730AG80.7% 58.9 / 13.7 31intergenic (‑209/+84)eutS/maeBethanolamine utilization microcompartment protein EutS/NADP‑dependent oxaloacetate‑decarboxylating malate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base G (16/9);  total (20/11)
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.40e-01

AAAGGTGAGGTAGATCACTAAGAATGGAGAGATATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGC  >  NZ_CP009273/2569246‑2569489
                                                                                                                               |                                                                                                                    
aaaGGTGAGGTAGATCACTAAGAATGGAGAGATATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGtaa                                                                                                            <  2:1192823/138‑1 (MQ=255)
                           agagaTATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTaaa                                                                                >  1:402937/1‑139 (MQ=255)
                            gagaTATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAAc                                                                               >  2:696993/1‑139 (MQ=255)
                                     gCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTAccc                                                                      >  1:258068/1‑139 (MQ=255)
                                         tggtgCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGaa                                                                  >  2:167422/1‑139 (MQ=255)
                                              cATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAat                                                                                              <  1:193512/106‑1 (MQ=255)
                                              cATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAat                                                                                              >  2:193512/1‑106 (MQ=255)
                                              cATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATAtt                                                              <  2:1150208/138‑1 (MQ=255)
                                                               aaCGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAatatat                                                                                          >  2:695829/1‑93 (MQ=255)
                                                               aaCGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAatatat                                                                                          <  1:695829/93‑1 (MQ=255)
                                                                 cGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGtt                                           >  1:895316/1‑138 (MQ=255)
                                                                                           gagaTATTATATAGACTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTg                 <  1:80517/138‑1 (MQ=255)
                                                                                                       aGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATAtt                                                              >  2:8708/1‑81 (MQ=255)
                                                                                                       aGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATAtt                                                              <  1:8708/81‑1 (MQ=255)
                                                                                                         aCTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGa                                                                   >  1:1258676/1‑74 (MQ=255)
                                                                                                         aCTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGa                                                                   <  2:1258676/74‑1 (MQ=255)
                                                                                                         aCTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  2:1194478/1‑139 (MQ=255)
                                                                                                         aCTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  2:324637/1‑139 (MQ=255)
                                                                                                          cTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  1:210380/1‑138 (MQ=255)
                                                                                                          cTGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  <  2:210380/138‑1 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGc                                                                            <  1:60592/63‑1 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGc                                                                            >  2:60592/1‑63 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  1:1069300/1‑137 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  2:745841/1‑137 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  1:697955/1‑137 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCtt                                                                                                     <  2:330898/38‑1 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCtt                                                                                                     >  1:330898/1‑38 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  1:317630/1‑137 (MQ=255)
                                                                                                           tGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGc  >  1:162382/1‑137 (MQ=255)
                                                                                                             cTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTg                 <  2:729910/120‑1 (MQ=255)
                                                                                                             cTTTTTTTTATTATTCGCGAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTg                 >  1:729910/1‑120 (MQ=255)
                                                                                                                               |                                                                                                                    
AAAGGTGAGGTAGATCACTAAGAATGGAGAGATATTCGCTGTGGTGCATAAACTTTCATATGCAACGCATGAATATAAATAAAATAATATCGAGATATTATATAGACTGCTTTTTTTTATTATTCGCAAACTCAGTAATAAGCTTTACTAATATATTGCAGCTAAACTGCTTACCCTGAATATTCAGGGTAAGCGTGAGAGTTAAAAAAAATTACAGCGGTTGGGTTTGCGCTTCTACCACGGC  >  NZ_CP009273/2569246‑2569489

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: