Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,505,565 A→C N436H (AAC→CAC)  ydcR → PLP‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,505,5650AC75.6% 95.7 / 28.7 41N436H (AAC→CAC) ydcRPLP‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (3/7);  new base C (15/16);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

AATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTTTAATTTAATTCCTCTTAGA  >  NZ_CP009273/1505432‑1505692
                                                                                                                                     |                                                                                                                               
aaTTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTgg                                                                                                                            >  1:734137/1‑139 (MQ=255)
          aaTGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGttttt                                                                                                                  >  2:80221/1‑139 (MQ=255)
                            ctctGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACcgcg                                                                                                 <  1:342534/138‑1 (MQ=255)
                            ctctGGTGGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACcgcg                                                                                                 <  1:744721/138‑1 (MQ=255)
                               tGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTCAATACCGCGTgg                                                                                              >  2:102388/1‑138 (MQ=255)
                                        cTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGgag                                                                                    <  2:249702/139‑1 (MQ=255)
                                           cccGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTc                                                                                                                <  2:725515/108‑1 (MQ=255)
                                           cccGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTc                                                                                                                >  1:725515/1‑108 (MQ=255)
                                                     tAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACcgcg                                                                                                 >  1:225188/1‑113 (MQ=255)
                                                     tAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACcgcg                                                                                                 <  2:225188/113‑1 (MQ=255)
                                                       gATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAgg                                                                     <  2:430366/139‑1 (MQ=255)
                                                           ccGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCgg                                                                  >  2:72917/1‑138 (MQ=255)
                                                            cGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGt                                                                 >  2:637021/1‑138 (MQ=255)
                                                              gCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTaa                                                               >  1:181190/1‑138 (MQ=255)
                                                                  aTTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACaa                                                          >  2:601665/1‑139 (MQ=255)
                                                                            gcggcACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTg                                                                                               <  2:533968/92‑1 (MQ=255)
                                                                            gcggcACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTg                                                                                               >  1:533968/1‑92 (MQ=255)
                                                                            gcggcACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCaaa                                                 <  1:477914/138‑1 (MQ=255)
                                                                            gcggcACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAAc                                                >  1:550562/1‑139 (MQ=255)
                                                                                      cGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGa                                      <  2:550562/139‑1 (MQ=255)
                                                                                      cGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGa                                      <  1:333139/139‑1 (MQ=255)
                                                                                        catcatATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAAc                                    <  1:80221/139‑1 (MQ=255)
                                                                                        catcatATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAAc                                    <  1:601665/139‑1 (MQ=255)
                                                                                          tcatATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACgg                                  >  2:71499/1‑139 (MQ=255)
                                                                                                     ttGCGCCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCg                       >  2:126123/1‑139 (MQ=255)
                                                                                                        cgcCGGGTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGtt                     >  1:313760/1‑138 (MQ=255)
                                                                                                        cgcCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTaa                                                               <  2:286552/96‑1 (MQ=255)
                                                                                                        cgcCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTaa                                                               >  1:286552/1‑96 (MQ=255)
                                                                                                             ggTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTg                                                                              >  1:438566/1‑76 (MQ=255)
                                                                                                             ggTAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTg                                                                              <  2:438566/76‑1 (MQ=255)
                                                                                                               tAAAATGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGtttaatttaa             >  2:564856/1‑139 (MQ=255)
                                                                                                                 aaaTGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTttaatttaat            <  1:126123/138‑1 (MQ=255)
                                                                                                                 aaaTGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTttaatttaat            <  1:71499/138‑1 (MQ=255)
                                                                                                                 aaaTGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTttaatttaat            <  2:344994/138‑1 (MQ=255)
                                                                                                                 aaaTGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTttaatttaat            <  1:849009/138‑1 (MQ=255)
                                                                                                                 aaaTGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTTtaatttaatt           <  1:578149/139‑1 (MQ=255)
                                                                                                                 aaaTGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTttaatttaat            <  1:556026/138‑1 (MQ=255)
                                                                                                                  aaTGTTTTCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTTTAATTTAATTc          <  2:734242/139‑1 (MQ=255)
                                                                                                                      ttttCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTTTAATTTAATTCctct      <  2:181190/139‑1 (MQ=255)
                                                                                                                       tttCTACCGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTTTAATTTAATTCctct      >  1:723123/1‑138 (MQ=255)
                                                                                                                             ccGGTGAACACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTTTAATTTAATTCCTCTTAGa  <  1:72917/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
AATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATATCAGTATTGCGCCGGGTAAAATGTTTTCTACCGGTGAAAACTGGTCACGTTTTTTCCGTTTTAATACCGCGTGGCAGTGGGGAGAGCGTGAAGAACAGGCGGTAAAACAATTAGGCAAACTTATTCAAGAACGGCTGTAATAGCGTTTAATTTAATTCCTCTTAGA  >  NZ_CP009273/1505432‑1505692

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: ‑