Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,337,867 A→C V246G (GTC→GGC)  yciT ← DNA‑binding transcriptional regulator YciT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,337,8670AC83.3% 87.8 / 13.9 36V246G (GTC→GGC) yciTDNA‑binding transcriptional regulator YciT
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base C (14/16);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

ATAATTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAGTATGAATGTACAGCACCAAATTTCGAGCTGT  >  NZ_CP009273/1337737‑1337999
                                                                                                                                  |                                                                                                                                    
ataatTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAAt                                                                                                                               <  1:836032/138‑1 (MQ=255)
 taatTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATgg                                                                                                                             <  2:353878/139‑1 (MQ=255)
 taatTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATgg                                                                                                                             <  1:719410/139‑1 (MQ=255)
       tataTCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGt                                                                                                                        <  2:708960/138‑1 (MQ=255)
              tGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACtgtg                                                                                                                >  1:630828/1‑139 (MQ=255)
               gagaGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACt                                                                                                                   <  2:866222/135‑1 (MQ=255)
                     gtTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTc                                                                                                             >  2:850586/1‑135 (MQ=255)
                                gTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCAt                                                                                                    >  2:728863/1‑133 (MQ=255)
                                gTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCAt                                                                                                    >  1:639866/1‑133 (MQ=255)
                                gTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCAt                                                                                                    >  1:816040/1‑133 (MQ=255)
                                   aaTAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATgg                                                                                                                             <  1:426627/105‑1 (MQ=255)
                                   aaTAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATgg                                                                                                                             >  2:426627/1‑105 (MQ=255)
                                        tCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATgg                                                                                                                             <  2:712806/100‑1 (MQ=255)
                                        tCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATgg                                                                                                                             >  1:712806/1‑100 (MQ=255)
                                        tCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGa                                                                                                         <  2:202985/120‑1 (MQ=255)
                                        tCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGa                                                                                                         >  1:202985/1‑120 (MQ=255)
                                        tCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGgcgc                                                                                      >  1:379527/1‑139 (MQ=255)
                                          ttAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGcgcg                                                                                     >  1:611697/1‑138 (MQ=255)
                                             aTTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAAt                                                                                                                               >  2:680692/1‑93 (MQ=255)
                                             aTTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATACGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAAt                                                                                                                               <  1:680692/93‑1 (MQ=255)
                                                        ggAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCgg                                                                       >  2:43223/1‑138 (MQ=255)
                                                                       cGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGaa                                                       >  2:568328/1‑139 (MQ=255)
                                                                         gatgatTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATc                                                     >  2:864986/1‑139 (MQ=255)
                                                                                 gCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAAc                                              <  2:40812/138‑1 (MQ=255)
                                                                                          cAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGTTCAAGATGCATCAGATCGCTGGGGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGa                                      <  2:660658/137‑1 (MQ=255)
                                                                                          cAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGa                                      <  1:864986/137‑1 (MQ=255)
                                                                                          cAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATg                                    <  2:639866/139‑1 (MQ=255)
                                                                                            aaTCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTACTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATg                                    >  1:304779/1‑137 (MQ=255)
                                                                                             atcatcCGTATAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAg                                 <  2:459540/139‑1 (MQ=255)
                                                                                             atcatcCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAg                                 <  2:816040/139‑1 (MQ=255)
                                                                                                atcCGTAGAGAAGGGAAATGGTTAAATGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAGTAt                              <  1:43223/139‑1 (MQ=255)
                                                                                                                       aaaTGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAg                                 <  2:812096/113‑1 (MQ=255)
                                                                                                                       aaaTGTCAACGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAg                                 >  1:812096/1‑113 (MQ=255)
                                                                                                                              aaCGCCGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAGTATGAATGTACAGCACCAAATTTCGAGCTgt  >  1:221324/1‑137 (MQ=255)
                                                                                                                              aaCGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCg                                       <  1:174830/100‑1 (MQ=255)
                                                                                                                              aaCGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGCATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCg                                       >  2:174830/1‑100 (MQ=255)
                                                                                                                                  |                                                                                                                                    
ATAATTTTATATCTTGAGAGTGTTAATAACAGGTAAATAGTCTTAATTATCAACCAGGAATCATCTTAGAGCGGATGATTTGCCAAACTGCAAATCATCCGTAGAGAAGGGAAATGGTTAAATGTCAACGACGTGAATGGTGAGTTTACTGTGCTCAAGATGCATCAGATCGCTGGCGCGTATTTTCGAATCGGTAATCACGCGATTGAATCGCTCAACGGGACCGATGGAGTATGAATGTACAGCACCAAATTTCGAGCTGT  >  NZ_CP009273/1337737‑1337999

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