Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,052,447 T→G K306Q (AAA→CAA)  yccM ← putative 4Fe‑4S membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,052,4470TG98.0% 178.2 / ‑4.5 50K306Q (AAA→CAA) yccMputative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (28/21);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGC  >  NC_000913/1052301‑1052593
                                                                                                                                                  |                                                                                                                                                  
cATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCttgttt                                                                                                                                                 <  2:512092/150‑1 (MQ=255)
 aTTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTa                                                                                                                                                >  2:70733/1‑150 (MQ=255)
    ggCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTAttt                                                                                                                                             <  2:547830/150‑1 (MQ=255)
     gcgcTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGc                                                                                                                                           <  2:442247/151‑1 (MQ=255)
       gcTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCg                                                                                                                                         <  1:824244/151‑1 (MQ=255)
           gtgGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAgg                                                                                                                                     <  2:182455/151‑1 (MQ=255)
            tgGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGt                                                                                                                                    >  1:173806/1‑151 (MQ=255)
             ggAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGtt                                                                                                                                   >  1:322474/1‑151 (MQ=255)
               aaTTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGtttt                                                                                                                                 <  2:322353/151‑1 (MQ=255)
                   aaGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGc                                                                                                                              >  1:532943/1‑150 (MQ=255)
                   aaGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCa                                                                                                                             <  2:375689/151‑1 (MQ=255)
                    aGCGTCAGTACAATACCTCCGGCGCAGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAAGCTCCGCAAACATAATCCCCAGAACCAGTAGCGTCATTATCCAGCCAGCCAACGCAAAGGCTTTTTAATTTGACACAGGTTTATGCAg                                                                                                                            <  2:840453/151‑1 (MQ=255)
                        tCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGa                                                                                                                         >  2:886567/1‑150 (MQ=255)
                                aTTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCa                                                                                                                 >  1:426940/1‑150 (MQ=255)
                                aTTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCa                                                                                                                 >  1:701639/1‑150 (MQ=255)
                                      ccGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAAACAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGtt                                                                                                            <  1:562735/149‑1 (MQ=255)
                                            ccGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTAc                                                                                                    >  1:371637/1‑151 (MQ=255)
                                              gggCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCg                                                                                                  >  2:538240/1‑151 (MQ=255)
                                               ggCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCg                                                                                                  <  1:60738/150‑1 (MQ=255)
                                                    tGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATg                                                                                            >  2:364564/1‑151 (MQ=255)
                                                                 cATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGt                                                                               <  2:580737/151‑1 (MQ=255)
                                                                  aTACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTc                                                                              >  2:805882/1‑151 (MQ=255)
                                                                     cATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCAt                                                                            >  2:64023/1‑150 (MQ=255)
                                                                     cATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATa                                                                           <  1:459086/151‑1 (MQ=255)
                                                                       ttGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATAc                                                                          >  2:852024/1‑150 (MQ=255)
                                                                                    cACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTc                                                             <  1:125791/150‑1 (MQ=255)
                                                                                     aCCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTcc                                                            >  1:118123/1‑150 (MQ=255)
                                                                                      ccGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCa                                                           <  1:438243/150‑1 (MQ=255)
                                                                                      ccGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCAc                                                          >  2:793256/1‑151 (MQ=255)
                                                                                              aTAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGt                                                   >  1:783448/1‑150 (MQ=255)
                                                                                                aaTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTAcc                                                >  1:406070/1‑151 (MQ=255)
                                                                                                 aTCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCt                                               <  1:746824/151‑1 (MQ=255)
                                                                                                   ccccAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTg                                              <  1:351625/150‑1 (MQ=255)
                                                                                                    cccAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGaa                                            <  2:618595/151‑1 (MQ=255)
                                                                                                        gAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGaatta                                         >  2:241711/1‑150 (MQ=255)
                                                                                                         aaCCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGaattaat                                       <  1:120904/151‑1 (MQ=255)
                                                                                                          aCCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGaattaatt                                      <  2:283026/151‑1 (MQ=255)
                                                                                                           ccAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAAttt                                     >  1:570192/1‑151 (MQ=255)
                                                                                                                 gCGTCATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAc                               >  1:23622/1‑151 (MQ=255)
                                                                                                                     cATTAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTTAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCgg                            <  2:790934/150‑1 (MQ=255)
                                                                                                                        tAACCAGCCAGACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGAt                         <  1:279150/150‑1 (MQ=255)
                                                                                                                            cagccagACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTg                    >  2:716299/1‑151 (MQ=255)
                                                                                                                               ccagACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATg                 <  2:307136/151‑1 (MQ=255)
                                                                                                                               ccagACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATg                 <  1:247658/151‑1 (MQ=255)
                                                                                                                                 agACAACGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATggg               >  2:396057/1‑151 (MQ=255)
                                                                                                                                      aCGCAAAGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAAtt          >  2:46520/1‑151 (MQ=255)
                                                                                                                                          aaaGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAAtttttt      >  2:353966/1‑151 (MQ=255)
                                                                                                                                           aaGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTg     >  1:33588/1‑151 (MQ=255)
                                                                                                                                            aGGCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTgc    >  1:787903/1‑151 (MQ=255)
                                                                                                                                             ggCTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTgcg   >  2:607119/1‑151 (MQ=255)
                                                                                                                                               cTTGTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTgcgc  >  1:540235/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGACCGGGATCCGTGATGGGCAATTTTTTGCGC  >  NC_000913/1052301‑1052593

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