Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,083,334	T→C	Y155C (TAT→TGT)	speB ←	agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,083,334	0	T	C	100.0%	175.6 / NA	48	Y155C (TAT→TGT)	speB	agmatinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/26);  total (22/26)

GGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCAGCAGCTTCTCGGCG  >  NC_000913/3083191‑3083474                                                                                                                                                |                                                                                                                                             ggCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATc                                                                                                                                       <  2:431278/151‑1 (MQ=255)       cGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGtgtg                                                                                                                                 >  2:651605/1‑151 (MQ=255)          ccAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTggg                                                                                                                               <  1:233934/150‑1 (MQ=255)          ccAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGc                                                                                                                              <  2:560363/151‑1 (MQ=255)                cGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCaaa                                                                                                                        <  2:629596/151‑1 (MQ=255)                cGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCaaa                                                                                                                        <  1:366298/151‑1 (MQ=255)                    aaaGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAgtgt                                                                                                                    >  1:330456/1‑151 (MQ=255)                    aaaGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAgtgt                                                                                                                    >  2:327305/1‑151 (MQ=255)                     aaGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAgtgt                                                                                                                    >  2:193081/1‑150 (MQ=255)                                   aTCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTg                                                                                                     <  1:801546/151‑1 (MQ=255)                                   aTCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTg                                                                                                     <  1:480853/151‑1 (MQ=255)                                    tCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGc                                                                                                    >  2:289033/1‑151 (MQ=255)                                         cTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGaa                                                                                                >  1:210364/1‑150 (MQ=255)                                           cGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATg                                                                                             <  1:654093/151‑1 (MQ=255)                                               aCGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTc                                                                                         <  2:62864/151‑1 (MQ=255)                                                  aaTACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGc                                                                                       >  1:825541/1‑150 (MQ=255)                                                      ccAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGa                                                                                   <  2:664192/150‑1 (MQ=255)                                                        aaTCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGCTCGCACAGGGATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGa                                                                                   >  2:56258/1‑148 (MQ=255)                                                              cACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACgcagc                                                                          <  1:697438/151‑1 (MQ=255)                                                                cGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGcagcag                                                                        <  2:521492/151‑1 (MQ=255)                                                                cGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGcagcag                                                                        <  2:830634/151‑1 (MQ=255)                                                                 ggAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCagcagc                                                                       <  2:548641/151‑1 (MQ=255)                                                                 ggAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCagcagc                                                                       >  2:228459/1‑151 (MQ=255)                                                                  gAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcg                                                                      <  2:438265/151‑1 (MQ=255)                                                                      gactccgttCGAACACACCTTCTTTCGGCGAGCTATAGAACATAGTGCCGTGGTCAAACTCCCGACCGTTCGCACAGCTATCTGTGTACGCGTCAAATTGCCCCAGCGCCATTTAGCCGAAATACTTCGCATGATCCCGCAgcagcggca                                                                   <  1:523262/142‑1 (MQ=255)                                                                        ttCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAgcagcggcagcg                                                                >  1:218636/1‑151 (MQ=255)                                                                          cGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTa                                                              <  1:332647/151‑1 (MQ=255)                                                                               cGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACaa                                                          >  1:43189/1‑150 (MQ=255)                                                                                    aGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGt                                                    >  1:560553/1‑151 (MQ=255)                                                                                     gACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTc                                                   >  2:175533/1‑151 (MQ=255)                                                                                       ccTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac                                                 >  1:800716/1‑151 (MQ=255)                                                                                       ccTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcac                                                 <  2:111354/151‑1 (MQ=255)                                                                                        cttcttTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTcacc                                                >  1:39907/1‑151 (MQ=255)                                                                                           cttTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCaccacc                                             >  1:667529/1‑151 (MQ=255)                                                                                               cGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGaa                                          <  1:472811/150‑1 (MQ=255)                                                                                                  cgcgGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAagag                                      <  1:610231/151‑1 (MQ=255)                                                                                                  cgcgGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAag                                        >  2:854365/1‑149 (MQ=255)                                                                                                            aaCATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCt                            <  2:359550/151‑1 (MQ=255)                                                                                                            aaCATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCGTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCt                            <  1:774307/151‑1 (MQ=255)                                                                                                             aCATAGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCtt                           >  2:46502/1‑151 (MQ=255)                                                                                                                 aGTGCCGTGGTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTAcc                        <  1:851240/150‑1 (MQ=255)                                                                                                                         ggTCAAATTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGcagccag               <  1:432606/151‑1 (MQ=255)                                                                                                                             aaaTTCACAACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCagcagc           >  1:496350/1‑151 (MQ=255)                                                                                                                                   acaACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCAGCAGCTTCTCg     <  1:505022/151‑1 (MQ=255)                                                                                                                                    caACCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCAGCAGCTTCTCgg    >  1:361002/1‑151 (MQ=255)                                                                                                                                     aaCCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCAGCAGCTTCTCGGc   <  2:353299/151‑1 (MQ=255)                                                                                                                                      aCCGTTCGCACAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCAGCAGCTTCTCGGCg  >  1:22538/1‑151 (MQ=255)                                                                                                                                      aCCGCTCGCACAGGTATCGGTGTGGGCGTCAAAGTGCACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCAGCAGCTTCTCGGc   >  1:184681/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                             GGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCAAATTCACAACCGTTCGCATAGGTATCGGTGTGGGCGTCAAAGTGTACCAGCGCCATTTTGCCGAAATGCTTCGCATGAGCACGCAGCAGCGGCAGCGTAACAAAGTGGTCACCACCGAAAGAGAGCATACGCTTACCGGCAGCCAGCAGCTTCTCGGCG  >  NC_000913/3083191‑3083474

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