Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,440,924	C→A	G36C (GGC→TGC)	rpoA ←	RNA polymerase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,440,924	0	C	A	100.0%	161.4 / NA	46	G36C (GGC→TGC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/23);  total (23/23)

CAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGA  >  NC_000913/3440793‑3441063                                                                                                                                    |                                                                                                                                            cAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAg                                                                                                                            >  1:96310/1‑149 (MQ=255)                    ccTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                         <  8:132524/149‑1 (MQ=255)                                      tCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                         >  7:159065/1‑131 (MQ=255)                                      tCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                         <  8:159065/131‑1 (MQ=255)                                      tCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTg                                                                                       >  2:86004/1‑149 (MQ=255)                                        aTGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGc                                                                                      <  1:86004/148‑1 (MQ=255)                                                 acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                 <  1:62085/144‑1 (MQ=255)                                                 acCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                 >  2:62085/1‑144 (MQ=255)                                                  ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                         >  6:177180/1‑119 (MQ=255)                                                  ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                         <  5:177180/119‑1 (MQ=255)                                                  ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                    >  7:49976/1‑140 (MQ=255)                                                  ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                    <  8:49976/140‑1 (MQ=255)                                                  ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACAGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                           >  4:45624/1‑149 (MQ=255)                                                   cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                           >  7:127909/1‑148 (MQ=255)                                                         tctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCa                                                                         >  8:9662/1‑144 (MQ=255)                                                         tctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCa                                                                         <  7:9662/144‑1 (MQ=255)                                                          ctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCCCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                   >  7:38187/1‑149 (MQ=255)                                                             aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtt                                                                >  3:1163/1‑149 (MQ=255)                                                                   ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACAGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAg                                                            <  3:45624/147‑1 (MQ=255)                                                                    gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGAGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGaaa                                                         >  7:61970/1‑149 (MQ=255)                                                                       aCCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                           >  6:32830/1‑128 (MQ=255)                                                                       aCCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                           <  5:32830/128‑1 (MQ=255)                                                                         cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAc                                                                                          <  7:71628/111‑1 (MQ=255)                                                                         cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAc                                                                                          >  8:71628/1‑111 (MQ=255)                                                                          gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGt                                                   <  2:96310/149‑1 (MQ=255)                                                                          gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGt                                                   <  8:127909/149‑1 (MQ=255)                                                                             cAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActct                                                     <  8:151056/144‑1 (MQ=255)                                                                             cAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTgcgat                                                >  7:151056/1‑145 (MQ=255)                                                                                      atggATGAGAGCAGAATACGGAGCAGTGCGTTACCCAGAGTATGTCAAAAGCCACGCTCTAAAGGCTCAAGGTTCACCTTGGCGTGCGTCGAACTCa                                                                                           <  8:117727/94‑1 (MQ=255)                                                                                      aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCa                                                                                           >  7:117727/1‑97 (MQ=255)                                                                                         gATGAGAGCAGAATACGGCGCTGTGCGTTACCCAGAGTATGGCAAAAGCCACGCGCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAccc                                         <  7:75800/144‑1 (MQ=255)                                                                                         gATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGACACGCTCTAAAGGCACAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGAAATCACCCAGGCGCCGATATAGAAACTCTGTCACAGAACCCtagaa                                    >  8:75800/1‑145 (MQ=255)                                                                                            gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCt                                                                                                             >  1:123606/1‑73 (MQ=255)                                                                                            gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCt                                                                                                             <  2:123606/73‑1 (MQ=255)                                                                                             agagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgt                                 <  4:1163/148‑1 (MQ=255)                                                                                              gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                                                                                        <  7:29234/92‑1 (MQ=255)                                                                                              gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACtt                                                                                        >  8:29234/1‑92 (MQ=255)                                                                                                    aaTACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                              <  2:32167/144‑1 (MQ=255)                                                                                                    aaTACGGCGCAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                              >  1:32167/1‑144 (MQ=255)                                                                                                            gcAGTGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTAc                 >  3:69819/1‑149 (MQ=255)                                                                                                                tGCGTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAg             >  2:77583/1‑149 (MQ=255)                                                                                                                   gTTACCCAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCtt          <  8:61970/149‑1 (MQ=255)                                                                                                                       cccAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt      >  6:82597/1‑149 (MQ=255)                                                                                                                         cAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg    <  5:51752/149‑1 (MQ=255)                                                                                                                         cAGAGTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg    <  4:69819/149‑1 (MQ=255)                                                                                                                           gagTATGGCAAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGa  <  8:38187/149‑1 (MQ=255)                                                                                                                                    |                                                                                                                                            CAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGA  >  NC_000913/3440793‑3441063

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