Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 791,624 T→C S103P (TCG→CCG)  modB → molybdate ABC transporter permease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273791,6240TC96.5% 90.5 / ‑3.2 29S103P (TCG→CCG) modBmolybdate ABC transporter permease subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (14/14);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

ACCGCCCGTGGTCGTCGGTTACTTATTATTAGTTTCGATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGTCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTATCACGTTACCGCTGACCTTACC  >  NZ_CP009273/791488‑791757
                                                                                                                                        |                                                                                                                                     
aCCGCCCGTGGTCGTCGGTTACTTATTATTAGTTTCGATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCg                                                                                                                                     >  1:452202/1‑139 (MQ=255)
              tcgGTTACTTATTATTAGTTTCGATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGt                                                                                                                       <  2:78143/139‑1 (MQ=255)
                ggTTACTTATTATTAGTTTCGATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTg                                                                                                                      <  1:65755/138‑1 (MQ=255)
                            ttaGTTTCGATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTc                                                                                                         <  1:36143/139‑1 (MQ=255)
                                  tCGATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGgcg                                                                                                    <  1:464313/138‑1 (MQ=255)
                                    gATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCACGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTg                                                                                                 >  2:400502/1‑139 (MQ=255)
                                      tGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGa                                                                                               <  1:402197/139‑1 (MQ=255)
                                            ggCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGt                                                                                         <  2:243228/139‑1 (MQ=255)
                                              cgcgGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGtt                                                                                        >  1:210234/1‑138 (MQ=255)
                                                                      tATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGc                                                                 >  2:92963/1‑137 (MQ=255)
                                                                        tggcTGGTGTGGTATTACCTTCGCCTTTAGCTGGAGCGGCGCAGTTCTCGATGCCGCAGTGACGCCGTTTCCGCTGATGGTGCGGGCAATACGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGa                                                             <  2:236854/136‑1 (MQ=255)
                                                                              gTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTg                                                       >  2:538365/1‑139 (MQ=255)
                                                                                ttGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTggg                                                     >  2:570110/1‑139 (MQ=255)
                                                                                        aCCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGgc                                             >  1:536978/1‑139 (MQ=255)
                                                                                        aCCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGgc                                             >  1:507596/1‑139 (MQ=255)
                                                                                             cGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTgg                                        >  1:131387/1‑139 (MQ=255)
                                                                                              gCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGgc                                       >  2:302151/1‑139 (MQ=255)
                                                                                                cTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGgcgc                                     <  2:210234/139‑1 (MQ=255)
                                                                                                         gcgcggcgcgGTTCTCGCTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCt                              <  1:538365/137‑1 (MQ=255)
                                                                                                              gcgcggTTCTCGCTGCCGCCGTCATGCCGTTTCCGCGGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTAt                       <  2:429522/139‑1 (MQ=255)
                                                                                                                          cTGCCGCCGTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTATCACGTTACCGCt           <  1:92963/139‑1 (MQ=255)
                                                                                                                              cgccgTCATGTCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTATCACGTTACCGCTGAcc       >  2:536867/1‑139 (MQ=255)
                                                                                                                              cgccgTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTATCACGTTACCGCTGAcc       >  2:489198/1‑139 (MQ=255)
                                                                                                                              cgccgTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTATCACGTTACCGCTGAcc       >  2:536465/1‑139 (MQ=255)
                                                                                                                                ccgTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCt                              <  2:183495/114‑1 (MQ=255)
                                                                                                                                ccgTCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCt                              >  1:183495/1‑114 (MQ=255)
                                                                                                                                   tCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCt                              <  2:226757/111‑1 (MQ=255)
                                                                                                                                   tCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCt                              >  1:226757/1‑111 (MQ=255)
                                                                                                                                   tCATGCCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTATCACGTTACCGCTGACCTTAcc  <  1:216971/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
ACCGCCCGTGGTCGTCGGTTACTTATTATTAGTTTCGATGGGACGGCGCGGATTTATCGGTGAACGTCTGTATGACTGGTTTGGTATTACCTTCGCCTTTAGCTGGCGCGGCGCGGTTCTCGCTGCCGCCGTCATGTCGTTTCCGCTGATGGTGCGGGCAATTCGTCTGGCGCTGGAAGGGGTTGATGTCAAACTGGAACAGGCCGCAAGAACACTGGGGGCCGGGCGCTGGCGCGTTTTCTTTACTATCACGTTACCGCTGACCTTACC  >  NZ_CP009273/791488‑791757

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