Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,447,365	A→G	V25A (GTC→GCC)	rplE ←	50S ribosomal subunit protein L5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,447,365	0	A	G	100.0%	220.9 / NA	60	V25A (GTC→GCC)	rplE	50S ribosomal subunit protein L5
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/34);  total (26/34)

CTTTGGTGATCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGACAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACTTGAAGAAACGGACTTTTTTA  >  NC_000913/3447231‑3447498                                                                                                                                       |                                                                                                                                      cTTTGGTGATCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGaa                                                                                                                                  >  2:781554/1‑140 (MQ=255)  tttGGTGATCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAAt                                                                                                                                 <  1:664535/140‑1 (MQ=255)   ttGGTGATCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAAtt                                                                                                                                <  1:45469/140‑1 (MQ=255)    tGGTGATCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTg                                                                                                                               <  1:768817/140‑1 (MQ=255)    tGGTGATCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTg                                                                                                                               >  1:714660/1‑140 (MQ=255)         aTCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGtt                                                                                                                          <  1:192188/140‑1 (MQ=255)          tCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTa                                                                                                                         >  1:955261/1‑140 (MQ=255)           cAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTaa                                                                                                                        <  2:710950/140‑1 (MQ=255)            aGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTaaa                                                                                                                       >  2:141014/1‑140 (MQ=255)              cGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACt                                                                                                                     >  1:1019337/1‑140 (MQ=255)              cGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACt                                                                                                                     >  1:1157489/1‑140 (MQ=255)              cGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTAAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACt                                                                                                                     <  1:528040/140‑1 (MQ=255)                      gACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATg                                                                                                             <  1:152332/140‑1 (MQ=255)                      gACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATg                                                                                                             <  1:859369/140‑1 (MQ=255)                        ccGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAg                                                                                                           <  2:443900/140‑1 (MQ=255)                         cGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGt                                                                                                          <  2:900592/140‑1 (MQ=255)                               ttGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTa                                                                                                    <  2:1157489/140‑1 (MQ=255)                                      cAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACtt                                                                                             <  2:40678/140‑1 (MQ=255)                                      cAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACtt                                                                                             <  2:343341/140‑1 (MQ=255)                                      cAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACtt                                                                                             <  1:141014/140‑1 (MQ=255)                                      cAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACtt                                                                                             <  2:1019337/140‑1 (MQ=255)                                      cAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACtt                                                                                             <  1:999718/140‑1 (MQ=255)                                      cAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACtt                                                                                             <  1:781554/140‑1 (MQ=255)                                        ggTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCATCGATCGCTTCACCAACACCCATGTTCAGGGCGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGATTTTTAACTACTTCg                                                                                           >  1:272388/1‑140 (MQ=255)                                         gTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGt                                                                                          <  1:991443/140‑1 (MQ=255)                                           ctgctgctgCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCt                                                                                        <  1:1024456/140‑1 (MQ=255)                                           ctgctgctgCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCt                                                                                        <  2:492422/140‑1 (MQ=255)                                            tgctgctgcGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCtt                                                                                       >  2:560099/1‑140 (MQ=255)                                              ctgctgcGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGCTTTTTAACTACTTCGTCTTTg                                                                                     >  2:1058141/1‑140 (MQ=255)                                               tgctgcGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTgt                                                                                    >  2:990506/1‑140 (MQ=255)                                                   gcGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTgtagta                                                                                >  2:655677/1‑140 (MQ=255)                                                   gcGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTgtagta                                                                                >  2:1079538/1‑140 (MQ=255)                                                          ccAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGc                                                                         <  1:774838/140‑1 (MQ=255)                                                                gTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTc                                                                   <  1:655677/140‑1 (MQ=255)                                                                  tttttGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGc                                                                 <  2:955261/140‑1 (MQ=255)                                                                  ttcttGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGc                                                                 <  2:714660/137‑1 (MQ=255)                                                                    tttGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCa                                                               <  1:1079538/140‑1 (MQ=255)                                                                     ttGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCAt                                                              <  2:810938/140‑1 (MQ=255)                                                                        tCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGt                                                           >  1:993985/1‑140 (MQ=255)                                                                                 gCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCa                                                  <  1:560099/140‑1 (MQ=255)                                                                                          acacCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTg                                         >  1:67694/1‑140 (MQ=255)                                                                                           cacCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGa                                        >  1:1013135/1‑140 (MQ=255)                                                                                             cccATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCaaa                                                <  1:1041512/130‑1 (MQ=255)                                                                                             cccATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCaaa                                                >  2:1041512/1‑130 (MQ=255)                                                                                             cccATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATa                                      >  2:115213/1‑140 (MQ=255)                                                                                                aTGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGtt                                   >  1:344581/1‑140 (MQ=255)                                                                                                aTGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGtt                                   <  2:67694/140‑1 (MQ=255)                                                                                                  gTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCa                                                  <  2:635420/123‑1 (MQ=255)                                                                                                  gTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCa                                                  >  1:635420/1‑123 (MQ=255)                                                                                                        ggTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACtt                                          >  1:360074/1‑125 (MQ=255)                                                                                                        ggTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACtt                                          <  2:360074/125‑1 (MQ=255)                                                                                                        ggTGATCTTCTCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGtt                           >  2:719784/1‑140 (MQ=255)                                                                                                             tcttctCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTc                                 >  1:760239/1‑129 (MQ=255)                                                                                                             tcttctCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTc                                 <  2:760239/129‑1 (MQ=255)                                                                                                             tcttctCGACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACt                      >  1:569763/1‑140 (MQ=255)                                                                                                                    gACCCGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACTTgaagaa               >  1:618638/1‑140 (MQ=255)                                                                                                                       ccGAGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACTTGAAGAAACg            <  2:818745/140‑1 (MQ=255)                                                                                                                          aGGGACTTGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACTTGAAGAAACGGAc         >  1:193076/1‑140 (MQ=255)                                                                                                                                ttGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACTTGAAGAAACGGACtttttt   >  1:1075573/1‑140 (MQ=255)                                                                                                                                 tGCATGGCAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACTTGAAGAAACGGACTTTTTTa  <  2:993985/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      CTTTGGTGATCAGCGGTTTTTGACCGGAGATTGCTGCCAGGTCTGCTGCTGCGTTATCCAGCAGTTTTTTGTCAGCGATCGCTTCACCAACACCCATGTTCAGGGTGATCTTCTCGACCCGAGGGACTTGCATGACAGAATTGTAGTTAAACTCAGTCATGAGTTTTTTAACTACTTCGTCTTTGTAGTAATCATGCAGTTTCGCCATCGTACTACTCCAAATTACTTGATAGTTTCGCTGTTAGACTTGAAGAAACGGACTTTTTTA  >  NC_000913/3447231‑3447498

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: ‑