Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,221,741	C→T	R249R (CGG→CGA)	arpA ←	ankyrin repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,221,741	0	C	T	100.0%	212.4 / NA	56	R249R (CGG→CGA)	arpA	ankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/26);  total (30/26)

AATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGG  >  NC_000913/4221593‑4221883                                                                                                                                                     |                                                                                                                                               aaTCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCg                                                                                                                                              <  1:179695/151‑1 (MQ=255)               aaaTCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCAt                                                                                                                                <  2:6155/151‑1 (MQ=255)                  tCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGt                                                                                                                             >  1:768935/1‑151 (MQ=255)                   ccATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGtt                                                                                                                            <  1:368537/151‑1 (MQ=255)                     aTCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTAt                                                                                                                          >  2:179833/1‑151 (MQ=255)                        agagAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATcca                                                                                                                       <  1:452023/151‑1 (MQ=255)                          agaATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATccacc                                                                                                                     >  2:132194/1‑151 (MQ=255)                          agaATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATccacc                                                                                                                     >  2:216997/1‑151 (MQ=255)                           gaATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATccacc                                                                                                                     >  2:91777/1‑150 (MQ=255)                               aGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTa                                                                                                                 >  2:512840/1‑150 (MQ=255)                                   tttAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATg                                                                                                            >  1:216676/1‑151 (MQ=255)                                     tAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGc                                                                                                          >  1:273759/1‑151 (MQ=255)                                      aaTATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCt                                                                                                         >  2:654899/1‑151 (MQ=255)                                        tatTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGc                                                                                                       >  2:606061/1‑151 (MQ=255)                                          ttCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGcc                                                                                                      >  1:419964/1‑150 (MQ=255)                                           tctcTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCa                                                                                                     <  1:189293/150‑1 (MQ=255)                                           tctcTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCAt                                                                                                    >  2:669157/1‑151 (MQ=255)                                              ctTTGATTTCGCTCATAGTTCTTGTTTTATCGTTAGTTACCTTAGGATCATAAACCGATACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAATCGGACCGTCAAAGCCTGGTTATCCACCAGTAATATGGCTGcc                                                                                                      >  1:357567/1‑146 (MQ=255)                                                ttGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTATTGGGTACATCCTTCAGTTGTGATCTTTATTCTTAGTCGGACCGTCAAAGTATGGTTATCCACCAGTAATATGGCTGCCATAaca                                                                                                >  1:698759/1‑150 (MQ=255)                                                  gCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACt                                                                                              >  1:660367/1‑150 (MQ=255)                                                  gCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTg                                                                                             <  2:563505/151‑1 (MQ=255)                                                    tttCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGAt                                                                                           >  2:360184/1‑151 (MQ=255)                                                       cGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGcc                                                                                        >  1:628870/1‑151 (MQ=255)                                                        gCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGcc                                                                                        >  2:703194/1‑150 (MQ=255)                                                        gCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCa                                                                                       <  2:532422/151‑1 (MQ=255)                                                        gCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCa                                                                                       <  2:306270/151‑1 (MQ=255)                                                           cATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCTAtt                                                                                     <  1:132161/150‑1 (MQ=255)                                                           cATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCAt                                                                                      <  1:560194/149‑1 (MQ=255)                                                                ttCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTc                                                                               >  2:96975/1‑151 (MQ=255)                                                                  cTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTcc                                                                              >  1:333320/1‑150 (MQ=255)                                                                   ttATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCAt                                                                            >  1:507772/1‑151 (MQ=255)                                                                    tATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTcc                                                                              <  1:300148/148‑1 (MQ=255)                                                                    tATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCa                                                                             <  2:412185/149‑1 (MQ=255)                                                                          aTCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTc                                                                      >  2:590219/1‑150 (MQ=255)                                                                          aTCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCa                                                                     <  2:181694/151‑1 (MQ=255)                                                                              ttagttaCATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAat                                                                 >  2:459380/1‑151 (MQ=255)                                                                                 gttaCATTAGCATCATAAGCCGAAACAACATAATGAGTACATCCTTCAGTTGTGTTCTTTACGCTTAGTCGAACCGTCAAAGCATGGTTAGCCGCCAGAACTATGGCTGCAATAACACTGATGCCATTTTGTTCCAAATTGTCAAAat                                                                 <  1:483786/148‑1 (MQ=255)                                                                                    aCATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCt                                                            <  1:493291/150‑1 (MQ=255)                                                                                      aTTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTc                                                         <  1:258856/151‑1 (MQ=255)                                                                                      aTTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTc                                                         <  1:136362/151‑1 (MQ=255)                                                                                        tAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCa                                                        <  1:585289/150‑1 (MQ=255)                                                                                          ggATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATa                                                      >  2:78706/1‑150 (MQ=255)                                                                                               aTAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGc                                                <  1:606004/151‑1 (MQ=255)                                                                                                      gAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAAt                                         >  2:6923/1‑151 (MQ=255)                                                                                                      gAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAAt                                         >  1:12987/1‑151 (MQ=255)                                                                                                           cACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCg                                    <  2:666467/151‑1 (MQ=255)                                                                                                              gTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCg                                    >  1:762439/1‑148 (MQ=255)                                                                                                                aaTGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCgcaggc                               <  2:374790/151‑1 (MQ=255)                                                                                                                 aTGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCgcaggca                              <  2:581047/151‑1 (MQ=255)                                                                                                                    ggTACATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGgcagca                           <  2:447971/151‑1 (MQ=255)                                                                                                                       aCATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAgcaacg                        >  1:633736/1‑151 (MQ=255)                                                                                                                        cATCCTTCAGTTGTGTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAgcaacgc                       <  2:541627/151‑1 (MQ=255)                                                                                                                                   tgtgTTCTTTATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTc             <  2:177063/150‑1 (MQ=255)                                                                                                                                           ttATTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCt    >  1:657164/1‑151 (MQ=255)                                                                                                                                             aTTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTgg  >  2:236176/1‑151 (MQ=255)                                                                                                                                             aTTCTTAGTCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTg   <  1:10039/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               AATCTACATTCATAAAATCCATCAGAGAATAGTGTTTAATATTCTCTTTGCTTTCGCTCATAATTCTTATTTTATCGTTAGTTACATTAGGATCATAAACCGAAACCACGTAATGGGTACATCCTTCAGTTGTGTTCTTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGG  >  NC_000913/4221593‑4221883

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