Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,892,532 A→C G82G (GGT→GGG ygcU ← FAD‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,892,5320AC89.8% 159.1 / 7.7 49G82G (GGT→GGGygcUFAD‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (23/21);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCA  >  NZ_CP009273/2892399‑2892649
                                                                                                                                     |                                                                                                                     
ccAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCgg                                                                                                                   >  1:261958/1‑138 (MQ=255)
ccAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCgg                                                                                                                   >  2:583868/1‑138 (MQ=255)
ccAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCgg                                                                                                                   >  2:340045/1‑138 (MQ=255)
 cAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGt                                                                                                                  <  2:668870/138‑1 (MQ=255)
         caccacATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGcac                                                                                                         <  1:10672/139‑1 (MQ=255)
              caTTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGt                                                                                                    <  2:813597/139‑1 (MQ=255)
                   cgcCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAAtt                                                                                               >  2:166241/1‑139 (MQ=255)
                        tCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGt                                                                                                                  >  1:2639/1‑115 (MQ=255)
                        tCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGt                                                                                                                  <  2:2639/115‑1 (MQ=255)
                               cATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCAt                                                                                     <  1:764600/137‑1 (MQ=255)
                                      tcaatatcaataTTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAAt                                                                            >  1:710260/1‑139 (MQ=255)
                                            tcaatatTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGca                                                                      <  1:825431/139‑1 (MQ=255)
                                            ccaatatTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGca                                                                      <  1:647583/138‑1 (MQ=255)
                                               atatTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCacac                                                                   <  1:583868/139‑1 (MQ=255)
                                               atatTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCacac                                                                   <  1:415309/139‑1 (MQ=255)
                                                tatTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACg                                                                  >  1:482906/1‑139 (MQ=255)
                                                  ttAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACgg                                                                 <  1:340045/138‑1 (MQ=255)
                                                     aTGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACgg                                                                 >  2:200180/1‑135 (MQ=255)
                                                     aTGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACgg                                                                 <  1:200180/135‑1 (MQ=255)
                                                      tgatttgattCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATAc                                                            <  1:684073/139‑1 (MQ=255)
                                                         tttgattCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTg                                                         <  2:160222/139‑1 (MQ=255)
                                                          ttgattCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGc                                                        >  1:800927/1‑139 (MQ=255)
                                                          ttgattCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGc                                                        >  1:296327/1‑139 (MQ=255)
                                                                  tGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGct                                                       >  1:72398/1‑132 (MQ=255)
                                                                  tGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGct                                                       <  2:72398/132‑1 (MQ=255)
                                                                  tGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTgg                                                 >  1:19320/1‑138 (MQ=255)
                                                                   ggCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGaa                                               >  1:329511/1‑139 (MQ=255)
                                                                   ggCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGaa                                               >  1:861528/1‑139 (MQ=255)
                                                                       gAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAAcc                                             >  1:197524/1‑137 (MQ=255)
                                                                        aGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCg                                            >  2:510924/1‑137 (MQ=255)
                                                                        aGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCg                                            <  1:510924/137‑1 (MQ=255)
                                                                             tCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGcac                                                                     >  1:223855/1‑107 (MQ=255)
                                                                             tCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGcac                                                                     <  2:223855/107‑1 (MQ=255)
                                                                                    ccaccGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGc                              >  1:108601/1‑139 (MQ=255)
                                                                                       ccGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTg                                                         >  1:670910/1‑109 (MQ=255)
                                                                                       ccGAGTTGTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTg                                                         <  2:670910/109‑1 (MQ=255)
                                                                                               tCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCg                   >  1:680134/1‑139 (MQ=255)
                                                                                               tCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCg                   >  2:709826/1‑139 (MQ=255)
                                                                                                 tACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCg                                                                                                                    >  1:69543/1‑40 (MQ=255)
                                                                                                 tACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCg                                                                                                                    <  2:69543/40‑1 (MQ=255)
                                                                                                   caacaGTTTCCAGCCCCCCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCaa               <  1:338839/139‑1 (MQ=255)
                                                                                                   caacaGTTTCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCaa               <  2:680134/139‑1 (MQ=255)
                                                                                                         tttCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGtata         <  2:800927/139‑1 (MQ=255)
                                                                                                         tttCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGtat          >  1:326462/1‑138 (MQ=255)
                                                                                                         tttCCAGCCCACCTTCGGTGGCGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGtat          >  1:250991/1‑138 (MQ=255)
                                                                                                           tctatcccACCTTCGGTGGCGGAAGCACCGATACGCGGCACACCGTTAATTTTGTGCGCATTCATAAATTTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACAGCTGCTGGAATCGGCGGAGTATaaa       >  2:438949/6‑139 (MQ=255)
                                                                                                                cccACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCa  >  2:371014/1‑139 (MQ=255)
                                                                                                                              cGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGAc                                 >  2:760125/1‑94 (MQ=255)
                                                                                                                              cGGAAGCCCCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGAc                                 <  1:760125/94‑1 (MQ=255)
                                                                                                                                     aCCGGTACGCGGCACACCGTTAATTTTGTGCGCAt                                                                                     <  1:259490/35‑1 (MQ=255)
                                                                                                                                     aCCGGTACGCGGCACACCGTTAATTTTGTGCGCAt                                                                                     >  2:259490/1‑35 (MQ=255)
                                                                                                                                     |                                                                                                                     
CCAGCGGAACACCACATTGCGCCGTCGCCTGCATATTCTCAATATCAATATTAATGATTTGATTCATGGCGGAGCCGTCGAGCACCACCGAGTTTTCTACAACAGTTTCCAGCCCACCTTCGGTGGCGGAAGCACCGGTACGCGGCACACCGTTAATTTTGTGCGCATTCATAAAATTCAGCACACGGGATACTTGCTCTGTGGAACCGAGTTTTACGACCGCTGCCGGAATCGGCAAAGTATAAATGCCA  >  NZ_CP009273/2892399‑2892649

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