Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,759,782 G→A 7.8% A645A (GCC→GCT cntO_2 ← Metal‑pseudopaline receptor CntO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,759,7820GA 7.8% 160.5 / 4.3 51A645A (GCC→GCTcntO_2Metal‑pseudopaline receptor CntO
Reads supporting (aligned to +/- strand):  ref base G (15/32);  new base A (2/2);  total (17/34)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CGCTGAACGCCGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCTGCG  >  ACB122/2759652‑2759908
                                                                                                                                  |                                                                                                                              
cGCTGAACGCCGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCgg                                                                                                                            <  5:249213/135‑1 (MQ=255)
cGCTGAACGCCGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCgg                                                                                                                            <  6:133652/135‑1 (MQ=255)
cGCTGAACGCCGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCgg                                                                                                                            <  5:132044/135‑1 (MQ=255)
cGCTGAACGCCGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCgg                                                                                                                            <  7:203125/135‑1 (MQ=255)
      aCGCCGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCgg                                                                                                                      <  1:83447/135‑1 (MQ=255)
          cGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATa                                                                                                                  <  8:16390/135‑1 (MQ=255)
          cGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATa                                                                                                                  <  3:39404/135‑1 (MQ=255)
            aCTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAAc                                                                                                                <  5:25620/135‑1 (MQ=255)
            aCTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAAc                                                                                                                <  4:247327/135‑1 (MQ=255)
              tGTACATACTCGCGATCGAACAGGGTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGc                                                                                                              <  7:161367/135‑1 (MQ=255)
                 aCATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTAGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCt                                                                                                             <  5:178587/133‑1 (MQ=255)
                         gcgATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATa                                                                                                   <  5:66353/135‑1 (MQ=255)
                          cgATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATa                                                                                                   <  5:82499/134‑1 (MQ=255)
                          cgATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATATTCATAACGCTCATCCACATa                                                                                                   <  7:247529/134‑1 (MQ=255)
                             tCGCACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGc                                                                                                              <  8:171570/120‑1 (MQ=255)
                             tCGCACAGGTTGTTCAACAGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGc                                                                                                              >  7:171570/1‑120 (MQ=255)
                             tCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCgg                                                                                                                            >  8:48947/1‑106 (MQ=255)
                             tCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCgg                                                                                                                            <  7:48947/106‑1 (MQ=255)
                             tCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAgcgc                                                                                               >  5:116791/1‑135 (MQ=255)
                              cGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAgcgcg                                                                                              <  7:144029/135‑1 (MQ=255)
                               gAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGc                                                                                                              <  1:247581/118‑1 (MQ=255)
                               gAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGc                                                                                                              >  2:247581/1‑118 (MQ=255)
                                    ggTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCaa                                                                                        <  6:131004/135‑1 (MQ=255)
                                    ggTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCaa                                                                                        <  2:260321/135‑1 (MQ=255)
                                       tgttCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGc                                                                                     <  4:217532/135‑1 (MQ=255)
                                        gttCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCGGGGTATTGGCGGTATAGGCAAAACGCTCATCCACATAGCCCGCCCCAATGcc                                                                                    >  6:177262/1‑135 (MQ=255)
                                          tCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCgg                                                                                  >  6:97796/1‑135 (MQ=255)
                                                  aaCCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCAc                                                                           <  5:62234/134‑1 (MQ=255)
                                                    ccccAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACctg                                                                        >  4:247808/1‑135 (MQ=255)
                                                      ccAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACctgct                                                                      >  7:188913/1‑135 (MQ=255)
                                                            tCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCg                                                                 >  1:191561/1‑134 (MQ=255)
                                                                                cAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCaaa                                              <  5:177262/133‑1 (MQ=255)
                                                                                cAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCaaa                                              <  3:18390/133‑1 (MQ=255)
                                                                                  ggCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTAGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGtt                                          >  4:28305/1‑135 (MQ=255)
                                                                                  ggCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTAGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGtt                                          >  4:24395/1‑135 (MQ=255)
                                                                                                  aCCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAACCCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGc                          >  3:218081/1‑135 (MQ=255)
                                                                                                  aCCACGGTGTACCCGGGGGCGCTCTGGGTATTAGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGc                          <  3:28305/135‑1 (MQ=255)
                                                                                                  aCCACGGTGTACCCGGGGGCGCTCTGGGTATTAGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGc                          <  3:24395/135‑1 (MQ=255)
                                                                                                        gtgtACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTg                    <  6:116791/135‑1 (MQ=255)
                                                                                                                    gcgcTCTGGGTATTGGCGGTATCGGCGTAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCCACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTg        >  8:292733/1‑135 (MQ=255)
                                                                                                                      gcTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGcc      >  5:116122/1‑135 (MQ=255)
                                                                                                                       ctctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCt     <  5:84292/135‑1 (MQ=255)
                                                                                                                       ctctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCt     <  8:289088/135‑1 (MQ=255)
                                                                                                                       ctctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCt     <  2:32467/135‑1 (MQ=255)
                                                                                                                       ctctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCt     <  3:247808/135‑1 (MQ=255)
                                                                                                                       ctctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCt     <  3:188202/135‑1 (MQ=255)
                                                                                                                       ctctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGGGGCACGTTGGCCGGCCGGTTGCCt     <  5:97796/135‑1 (MQ=255)
                                                                                                                        tctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGc                                                                                     <  6:184728/54‑1 (MQ=255)
                                                                                                                        tctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGc                                                                                     >  5:184728/1‑54 (MQ=255)
                                                                                                                        tctGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCTg    >  1:180730/1‑135 (MQ=255)
                                                                                                                          tGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCTgcg  >  4:81780/1‑135 (MQ=255)
                                                                                                                                  |                                                                                                                              
CGCTGAACGCCGACTGTACATACTCGCGATCGAACAGGTTGTTCAACTGCAACCCCAACTTCACATCCGGCAACACCTGCCAGGCGATGTTGGCATCGACCACGGTGTACCCGGGGGCGCTCTGGGTATTGGCGGTATCGGCATAACGCTCATCCACATAGCGCGCCCCAATGCCGGCCTCCACCTGCTGGCCGAAACCTTTGCTCAACCAAAGGTTGGCGGTACGCCGTGGCACGTTGGCCGGCCGGTTGCCTGCG  >  ACB122/2759652‑2759908

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