Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,957,921	T→C	T322A (ACA→GCA)	ptrA ←	protease 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,957,921	0	T	C	100.0%	172.5 / NA	49	T322A (ACA→GCA)	ptrA	protease 3
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/27);  total (22/27)

CCTTTATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGTACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAA  >  NC_000913/2957803‑2958049                                                                                                                       |                                                                                                                                 ccTTTATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAAt                                                                                                             <  2:292504/140‑1 (MQ=255)  cTTTATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATc                                                                                                            >  1:832234/1‑140 (MQ=255)  cTTTATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATc                                                                                                            >  1:606650/1‑140 (MQ=255)   tttATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCa                                                                                                           <  1:204446/140‑1 (MQ=255)   tttATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCa                                                                                                           <  1:21702/140‑1 (MQ=255)     tATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCa                                                                                                               >  2:225016/1‑134 (MQ=255)     tATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCa                                                                                                               <  1:225016/134‑1 (MQ=255)     tATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGa                                                                                                         >  1:316495/1‑140 (MQ=255)      aTCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGAt                                                                                                        <  2:327027/140‑1 (MQ=255)              aaGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCa                                                                                                           >  2:161698/1‑129 (MQ=255)              aaGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCa                                                                                                           <  1:161698/129‑1 (MQ=255)                   gCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATAAGATAGGTAATCAATTc                                                                                           <  1:653109/140‑1 (MQ=255)                    cAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCa                                                                                          >  1:988008/1‑140 (MQ=255)                     agagATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCAt                                                                                         <  1:980204/140‑1 (MQ=255)                     agagATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCAt                                                                                         >  1:811170/1‑140 (MQ=255)                         aTCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGt                                                                                     >  1:570527/1‑140 (MQ=255)                             cTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTta                                                                                 >  1:829262/1‑140 (MQ=255)                                   cgccgcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACgg                                                                           >  1:914598/1‑140 (MQ=255)                                     ccgcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGaa                                                                         <  2:243096/140‑1 (MQ=255)                                      cgcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTtacta                                                                              <  2:292396/134‑1 (MQ=255)                                      cgcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTtacta                                                                              >  1:292396/1‑134 (MQ=255)                                       gcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACt                                                                       >  1:787567/1‑140 (MQ=255)                                       gcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACt                                                                       <  2:832234/140‑1 (MQ=255)                                       gcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACt                                                                       <  2:787567/140‑1 (MQ=255)                                       gcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACt                                                                       >  2:88663/1‑140 (MQ=255)                                       gcTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACt                                                                       <  2:606650/140‑1 (MQ=255)                                        cTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACtt                                                                      >  1:14086/1‑140 (MQ=255)                                          gTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCg                                                                    <  2:316495/140‑1 (MQ=255)                                            tGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCt                                                                  <  1:88663/140‑1 (MQ=255)                                            tGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCt                                                                  <  2:829262/140‑1 (MQ=255)                                             gCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCt                                                                  <  1:117060/139‑1 (MQ=255)                                             gCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCt                                                                  >  2:117060/1‑139 (MQ=255)                                              ccGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGa                                                                <  2:14086/140‑1 (MQ=255)                                               cGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAg                                                               <  2:338430/140‑1 (MQ=255)                                                gTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAgt                                                              >  1:834944/1‑140 (MQ=255)                                                   gACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAgttgt                                                           >  1:725091/1‑140 (MQ=255)                                                        tAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCg                                                      <  1:543216/140‑1 (MQ=255)                                                           gATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATg                                                   >  2:125808/1‑140 (MQ=255)                                                                           aaTGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcgcg                                     >  2:876303/1‑138 (MQ=255)                                                                           aaTGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAAcgcg                                     <  1:876303/138‑1 (MQ=255)                                                                            aTGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCaa                                  <  1:125808/140‑1 (MQ=255)                                                                                        aaTCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCgg                      >  2:50303/1‑140 (MQ=255)                                                                                         aTCCCTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGc                     >  1:973774/1‑140 (MQ=255)                                                                                             cTGCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATCGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAgcg                 <  2:914598/140‑1 (MQ=255)                                                                                               gCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCa               <  2:519288/140‑1 (MQ=255)                                                                                               gCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCa               <  2:297392/140‑1 (MQ=255)                                                                                               gCTTTTGCAGCCAGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCa               <  2:988008/140‑1 (MQ=255)                                                                                                           aGTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTaatga   >  1:251204/1‑140 (MQ=255)                                                                                                            gTCAGAAAGTGCACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTaatgaa  <  2:251204/140‑1 (MQ=255)                                                                                                                       |                                                                                                                                 CCTTTATCGGTTAAAGACGCAGAGATCGCTAATACGCCGCTGTTGCCGTTGACGATAGGATCGGAGTTGGCGCTAATGCCCTCAACTAATCCCTGCTTTTGCAGCCAGTCAGAAAGTGTACCTGGGCTGCGATTGCCAATCAGATAGGTAATCAATTCATCGGTTTTACTACGGAACTTCGCTGAGTTGTTATCGATGCGAAACTCAACGCGCAACACTTTACGCGGCAGCGCAGGGACGTAATGAA  >  NC_000913/2957803‑2958049

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