Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,521,245 C→T E69K (GAG→AAG)  hofO ← DNA utilization protein HofO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,521,2450CT95.5% 142.9 / ‑6.2 45E69K (GAG→AAG) hofODNA utilization protein HofO
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (16/26);  minor base G (2/0);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGC  >  NC_000913/3521117‑3521378
                                                                                                                                |                                                                                                                                     
tGCCGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAg                                                                                                                            <  2:34765/140‑1 (MQ=255)
tGCCGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAg                                                                                                                            <  1:473427/140‑1 (MQ=255)
 gCCGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGc                                                                                                                           >  2:720528/1‑140 (MQ=255)
                 cccACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGa                                                                                                           >  2:183754/1‑140 (MQ=255)
                   cacaGCGTTTTCAACGCCAACTCGCCTCCCTGCGCTGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATa                                                                                                         <  2:519831/140‑1 (MQ=255)
                   cacaGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATa                                                                                                         <  2:439310/140‑1 (MQ=255)
                        cGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCGGCGTTTTTGCCTGGCTGAAGGGAGCGGTGTGTACCAGGCGATACAGTg                                                                                                    >  2:245340/1‑140 (MQ=255)
                        cGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTAcc                                                                                                                 >  1:722735/1‑127 (MQ=255)
                        cGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTAcc                                                                                                                 <  2:722735/127‑1 (MQ=255)
                        cGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTg                                                                                                    <  1:134819/140‑1 (MQ=255)
                        cGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTg                                                                                                    >  1:805088/1‑140 (MQ=255)
                        cGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTg                                                                                                    <  1:586445/140‑1 (MQ=255)
                        cGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTGCCTGGCTGAAGGGAGCGGTGTATACCAGGCGATACAGTg                                                                                                    >  2:312025/1‑140 (MQ=255)
                          ttttCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGAGGTGTCTACGAGGCGATACAGTGCg                                                                                                  >  1:352540/1‑140 (MQ=255)
                          ttttCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGAGGAGTCTACTAGGCGATACAGTGCg                                                                                                  >  1:266513/1‑140 (MQ=255)
                           tttCAACGCCAACTCTCCTCCCTGCCCGGTTGGATGCTAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCgg                                                                                                 <  1:294423/140‑1 (MQ=255)
                                aCGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCAc                                                                                            <  2:635849/140‑1 (MQ=255)
                                  gcaaaCTCGCCTCCCTGCGCGGCTGGATGTCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCActg                                                                                          <  1:291343/137‑1 (MQ=255)
                                            cTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGtgatga                                                                                >  2:387244/1‑140 (MQ=255)
                                                cTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTg                                                                            >  2:54769/1‑140 (MQ=255)
                                                  gcgcGGATGGATGCCAGGAAACCAGGTGCACGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGcc                                                                          <  1:312025/140‑1 (MQ=255)
                                                   gagacagggatgCCAGGGAACCTGTTGCGCGCCGGGTAATTGAAAATCAAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGccc                                                                         <  2:351134/133‑1 (MQ=255)
                                                    gcgGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACtt                                                                             >  2:826308/1‑135 (MQ=255)
                                                    gcgGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACtt                                                                             <  1:826308/135‑1 (MQ=255)
                                                    gcgGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCg                                                                        <  1:54769/140‑1 (MQ=255)
                                                         tggatgCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTc                                                                   <  1:720528/140‑1 (MQ=255)
                                                               ccAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTa                                                             <  1:183754/140‑1 (MQ=255)
                                                                  ggAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATg                                                          >  2:549742/1‑140 (MQ=255)
                                                                    aaaCCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCg                                                        <  2:805088/140‑1 (MQ=255)
                                                                     aaCCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGCCTCGCCCGTTGTCGAAttctgtctc                                                       >  1:351134/1‑133 (MQ=255)
                                                                        cAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTcgc                                                    <  2:662191/140‑1 (MQ=255)
                                                                         aGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTcgcg                                                   <  2:93659/140‑1 (MQ=255)
                                                                                cgcCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCa                                            <  1:549742/140‑1 (MQ=255)
                                                                                             aaaaTCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATg                               >  1:822494/1‑140 (MQ=255)
                                                                                             aaaaTCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATg                               <  1:387244/140‑1 (MQ=255)
                                                                                               aaTCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGAc                             >  1:781739/1‑140 (MQ=255)
                                                                                                aTCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTTCCCGTAGTGGAATTAAGGCGTCGCGTTCTTCATGGTGTGTCGATGACa                            >  2:295763/1‑140 (MQ=255)
                                                                                                      tGGCGAAAAGGGCAGCGTTTTTGCCATGCGGAAGGGAACGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGaaaaa                      <  1:407395/140‑1 (MQ=255)
                                                                                                            aaaGGGCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCg                <  2:439782/140‑1 (MQ=255)
                                                                                                               gggCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGtaa             <  1:566542/140‑1 (MQ=255)
                                                                                                                ggCAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGtaac            >  2:774560/1‑140 (MQ=255)
                                                                                                                  cAGCGTTTTTTCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTAGAGGACAGAAAAATGAGCGtaacta          >  1:350332/1‑140 (MQ=255)
                                                                                                                       ttttttCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATc     <  2:253/140‑1 (MQ=255)
                                                                                                                       ttttttCCTTGCTGAAGAGAGCAGTGTCTACTAGGAGATCCAGTGCGGACCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATc     <  2:266513/140‑1 (MQ=255)
                                                                                                                          tttCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGc  >  1:781437/1‑140 (MQ=255)
                                                                                                                          tttCCTTGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGc  <  1:295763/140‑1 (MQ=255)
                                                                                                                                |                                                                                                                                     
TGCCGACGGCACTGCTTCCCACAGCGTTTTCAACGCCAACTCGCCTCCCTGCGCGGATGGATGCCAGGAAACCAGTTGCGCGCCGGATAACTGAAAATCCAGTGGCGAAAAGGGCAGCGTTTTTTCCTCGCTGAAGGGAGCGGTGTCTACCAGGCGATACAGTGCGGCCCACTGCTGGTGATGACTTGCCCGTAGTCGAATTAATGCGTCGCGCTCTTCATGGTGTGTCGATGACAGAAAAATGAGCGTAACTAACATCAGC  >  NC_000913/3521117‑3521378

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑