Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,052,428	C→T	G312D (GGC→GAC)	yccM ←	putative 4Fe‑4S membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,052,428	0	C	T	95.1%	127.6 / ‑1.6	41	G312D (GGC→GAC)	yccM	putative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (18/21);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAC  >  NC_000913/1052291‑1052564                                                                                                                                          |                                                                                                                                         aGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAAccagtca                                                                                                                                        >  2:740418/1‑140 (MQ=255) aGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAAccagtca                                                                                                                                        >  1:95510/1‑140 (MQ=255)     aCCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACa                                                                                                                                    >  2:169301/1‑140 (MQ=255)           cATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCaa                                                                                                                              <  2:940779/140‑1 (MQ=255)           cATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCaa                                                                                                                              <  2:196407/140‑1 (MQ=255)           cATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCaa                                                                                                                              <  2:170403/140‑1 (MQ=255)                    ttGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCaaa                                                                                                                             <  1:206057/132‑1 (MQ=255)                    ttGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCaaa                                                                                                                             >  2:206057/1‑132 (MQ=255)                             tAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCg                                                                                                            <  1:125361/140‑1 (MQ=255)                             tAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCg                                                                                                            <  2:274164/140‑1 (MQ=255)                               aGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCa                                                                                                          <  1:262186/140‑1 (MQ=255)                                                tccGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGaaaa                                                                                         >  2:435479/1‑140 (MQ=255)                                                 ccGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAg                                                                                        >  1:392087/1‑140 (MQ=255)                                                   ggCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGt                                                                                      <  1:169301/140‑1 (MQ=255)                                                               tGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGc                                                                          >  1:876845/1‑140 (MQ=255)                                                                 ccATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTa                                                                        <  2:355629/140‑1 (MQ=255)                                                                    tACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGt                                                                                      <  1:736519/123‑1 (MQ=255)                                                                    tACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGt                                                                                      >  2:736519/1‑123 (MQ=255)                                                                         ccGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACAt                                                                >  1:26682/1‑140 (MQ=255)                                                                           gCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGa                                                              <  2:973723/140‑1 (MQ=255)                                                                              tACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGActc                                                           >  1:197859/1‑140 (MQ=255)                                                                                   tGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCTGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGc                                                      >  1:593420/1‑140 (MQ=255)                                                                                   tGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACa                                                                 <  2:908255/129‑1 (MQ=255)                                                                                   tGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACa                                                                 >  1:908255/1‑129 (MQ=255)                                                                                   tGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTCTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGc                                                      <  2:866940/140‑1 (MQ=255)                                                                                     cGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCaa                                                    <  1:435479/140‑1 (MQ=255)                                                                                             aTCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATa                                              >  1:433968/1‑138 (MQ=255)                                                                                             aTCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATa                                              <  2:433968/138‑1 (MQ=255)                                                                                             aTCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACa                                            >  2:309367/1‑140 (MQ=255)                                                                                                 ccGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCg                                        <  2:95510/140‑1 (MQ=255)                                                                                                    cAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCAt                                               <  1:220328/130‑1 (MQ=255)                                                                                                    cAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCAt                                               >  2:220328/1‑130 (MQ=255)                                                                                                    cAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTa                                     >  2:596526/1‑140 (MQ=255)                                                                                                    cAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTa                                     >  2:23255/1‑140 (MQ=255)                                                                                                    cAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTa                                     >  2:778966/1‑140 (MQ=255)                                                                                                       aCATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACAt                                  <  1:596526/140‑1 (MQ=255)                                                                                                       aCATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACAt                                  <  1:23255/140‑1 (MQ=255)                                                                                                       aCATAATCCCCAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACAt                                  <  1:778966/140‑1 (MQ=255)                                                                                                                ccAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGtt                         <  2:26682/140‑1 (MQ=255)                                                                                                                 cAGAACCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTc                        >  1:300155/1‑140 (MQ=255)                                                                                                                     aCCAGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATTACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTAc                    <  2:876845/140‑1 (MQ=255)                                                                                                                        aGTAGCGTCATTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTg                 >  1:174170/1‑140 (MQ=255)                                                                                                                                 aTTAACCAGTCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAAttt        >  1:686606/1‑140 (MQ=255)                                                                                                                                       cagtcagACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAc  <  2:392087/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         AGTGACCAATCATTGGCGCTTGTGGAATTAAGCGTCGGTACAATTCCTCCGGCACCGGGCTTTGCCATACTCCCGCATACATTGCGTAACCAATCACCGCAAACATAATCCCCAGAACCAGTAGCGTCATTAACCAGCCAGACAACGCAAAGGCTTTTTTATTTGCCGCAGGTTTTTGCAGTGAAAAGGTCAATGTTGAGGCTACCGGACATGACTCTACGCAAGTCATACAGCCGGTACATTCCACTGTTCGTACCTGAATTAATTTATCGAC  >  NC_000913/1052291‑1052564

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