Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,446,157	G→A	W682* (TGG→TGA)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,446,157	0	G	A	100.0%	194.1 / NA	53	W682* (TGG→TGA)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/25);  total (28/25)

AGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACACTGCTGGCTTAA  >  NC_000913/4446079‑4446234                                                                               |                                                                              aGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATc                                                                             <  1:918990/81‑1 (MQ=255)  ggAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCg                                                                            >  2:919044/1‑81 (MQ=255)       ttttGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGAc                                                                       <  2:918991/81‑1 (MQ=255)          tGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCcg                                                                    >  2:919042/1‑81 (MQ=255)           gATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCcgc                                                                   >  1:919045/1‑81 (MQ=255)            aTCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCcgc                                                                   >  2:919053/1‑80 (MQ=255)             tCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGa                                                                 <  2:918992/81‑1 (MQ=255)              cGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGa                                                                 >  1:919047/1‑80 (MQ=255)               gCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGata                                                               >  2:919035/1‑81 (MQ=255)               gCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGata                                                               >  1:919063/1‑81 (MQ=255)                 aaaGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCTGGTTGGCCCGTGATCGCTGACCCGCGATAtt                                                             >  2:919059/1‑81 (MQ=255)                 aaaGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATAtt                                                             <  2:918993/81‑1 (MQ=255)                    gacgaaCGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTgcg                                                          <  2:918995/78‑1 (MQ=255)                    gaagaaCGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTgcg                                                          <  2:918994/81‑1 (MQ=255)                      agaaCGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCt                                                        <  1:918996/81‑1 (MQ=255)                       gaaCGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTg                                                       >  2:919058/1‑81 (MQ=255)                       gaaCGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTg                                                       <  2:918997/81‑1 (MQ=255)                        aaCGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTgg                                                      <  2:918998/81‑1 (MQ=255)                               gAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTAcc                                                <  2:918999/80‑1 (MQ=255)                               gAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTAc                                                 <  1:919000/79‑1 (MQ=255)                               gAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCg                                               >  2:919050/1‑81 (MQ=255)                                aaGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCg                                               >  1:919057/1‑80 (MQ=255)                                aaGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCg                                               >  2:919049/1‑80 (MQ=255)                                  gggAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCaa                                            >  2:919037/1‑81 (MQ=255)                                    gAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATa                                          <  1:919001/81‑1 (MQ=255)                                     aaCTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATaa                                         <  1:919002/81‑1 (MQ=255)                                            aGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGc                                    <  2:919003/79‑1 (MQ=255)                                            aGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCa                                   >  2:919071/1‑80 (MQ=255)                                            aGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCa                                   >  2:919039/1‑80 (MQ=255)                                              cAATACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCaaa                                 >  1:919065/1‑80 (MQ=255)                                               aaTACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCaaaaa                               <  1:919004/81‑1 (MQ=255)                                               aaTACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCaaaa                                >  2:919054/1‑80 (MQ=255)                                                aTACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCaaaaaa                              >  2:919077/1‑81 (MQ=255)                                                 tACCCGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAAt                             <  2:919005/81‑1 (MQ=255)                                                   cccGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCa                           <  2:919006/81‑1 (MQ=255)                                                    ccGCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAg                          <  2:919007/81‑1 (MQ=255)                                                      gCTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCa                        >  2:919062/1‑81 (MQ=255)                                                       cTTCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCAt                       <  2:919008/81‑1 (MQ=255)                                                         tCCAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCg                     >  1:919070/1‑81 (MQ=255)                                                          ccAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATc                      <  1:919009/79‑1 (MQ=255)                                                          ccAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCgg                    >  2:919068/1‑81 (MQ=255)                                                           cAGACGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCgg                    <  2:919010/80‑1 (MQ=255)                                                              aCGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAGAATCAGCATCGGGCca                >  1:919055/1‑81 (MQ=255)                                                               cGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCcac               <  1:919011/81‑1 (MQ=255)                                                               cGCCGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCca                >  1:919052/1‑80 (MQ=255)                                                                  cGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACAct            >  1:919094/1‑81 (MQ=255)                                                                  cGGTTGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACAct            >  2:919066/1‑81 (MQ=255)                                                                     ttGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACActgct         >  2:919061/1‑81 (MQ=255)                                                                     ttGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACActgct         <  1:919012/81‑1 (MQ=255)                                                                      tGGCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACActgctg        >  1:919064/1‑81 (MQ=255)                                                                       ggCCCGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACActgctg        <  2:919013/80‑1 (MQ=255)                                                                          ccGTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACACTGCTGGCtt    <  2:919014/81‑1 (MQ=255)                                                                            gTGATCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACACTGCTGGCTTaa  >  2:919056/1‑81 (MQ=255)                                                                               |                                                                              AGGAAGTTTTGATCGCAAAGAAGAACGCTGGAAGGGAACTCTTAGCAATACCCGCTTCCAGACGCCGGTTGGCCCGTGGTCGCTGACCCGCGATATTGCGCTGGATTACCGCAATAAGGAGCAAAAAATCAGCATCGGGCCACACTGCTGGCTTAA  >  NC_000913/4446079‑4446234

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