Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,743,603	Δ1 bp	coding (1171/1275 nt)	yiaN →	putative membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,743,603	0	C	.	100.0%	383.5 / NA	84	coding (1171/1275 nt)	yiaN	putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (34/50);  total (34/50)

TCGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATCGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATT  >  NC_000913/3743532‑3743672                                                                        |                                                                      tCGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGCGGCAAAACTCAAATTCGAT‑GATGc                                                                  >  2:1011003/1‑76 (MQ=255)  cGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCa                                                                 <  2:3926662/76‑1 (MQ=255)   ggCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAg                                                                <  1:587260/76‑1 (MQ=255)   ggCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAg                                                                >  2:3839732/1‑76 (MQ=255)     cTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTc                                                              >  2:5054459/1‑76 (MQ=255)      ttAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCa                                                             <  2:4284949/76‑1 (MQ=255)      ttAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCa                                                             >  2:3546800/1‑76 (MQ=255)      ttAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCa                                                             >  2:588793/1‑76 (MQ=255)       tAATTACTCCGCCTATCGGTAATGTACTTAGCGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCag                                                            <  2:5542915/76‑1 (MQ=255)       tAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCag                                                            <  2:2638233/76‑1 (MQ=255)        aaTTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCaga                                                           <  2:169149/76‑1 (MQ=255)         aTTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCagag                                                          <  1:1720580/76‑1 (MQ=255)         aTTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCagag                                                          <  1:3181080/76‑1 (MQ=255)           tACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGc                                                        <  1:4429940/76‑1 (MQ=255)            ccacCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCg                                                       <  2:2640657/75‑1 (MQ=255)            ccacCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCg                                                       <  2:1601786/75‑1 (MQ=255)             tacCGCCTATCGGTAACGTACTGAACGTTATTTCAGGGGAGGCAAAACTCAAATTAGAT‑GATGCAGTCAGAGGCGt                                                      >  1:705937/2‑76 (MQ=255)              cccGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGtt                                                     <  1:27291/75‑1 (MQ=255)              acCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGtt                                                     >  1:120958/1‑76 (MQ=255)               ccGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGttt                                                    <  2:1815965/76‑1 (MQ=255)               ccGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGttt                                                    <  2:4811728/76‑1 (MQ=255)                cGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGtttt                                                   <  2:1815337/76‑1 (MQ=255)                cGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGtttt                                                   <  1:3546800/76‑1 (MQ=255)                 gCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTc                                                  >  2:1448226/1‑76 (MQ=255)                    tATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCt                                               >  2:4511135/1‑76 (MQ=255)                     aTCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCtt                                              >  1:3166617/1‑76 (MQ=255)                      tCGGTAATGTACTTAACGTTATTTCAGGTGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTa                                             <  2:217479/76‑1 (MQ=255)                      tCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTa                                             <  2:1609398/76‑1 (MQ=255)                      tCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTa                                             <  2:1927531/76‑1 (MQ=255)                       cGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTAc                                            <  2:1708397/76‑1 (MQ=255)                       cGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTAc                                            <  2:1104982/76‑1 (MQ=255)                        ggTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACg                                           >  1:5342526/1‑76 (MQ=255)                        ggTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACg                                           <  1:4635704/76‑1 (MQ=255)                        ggTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACg                                           >  1:1666106/1‑76 (MQ=255)                         gTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGt                                          >  1:3599704/1‑76 (MQ=255)                          tAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTc                                         >  1:390744/1‑76 (MQ=255)                            aTGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCt                                       <  2:5479796/76‑1 (MQ=255)                            aTGTACTTAACGTTATTGTGGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCt                                       <  1:1386796/76‑1 (MQ=255)                             tGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTg                                      <  1:5033822/76‑1 (MQ=255)                              gTACTTAACGTTATTTCAGTTGTGGCAAAACTCAAATTCTAT‑GATGCAGTCAGCGGCGTTTTCCCTTCCGTCCTgg                                     <  1:231021/76‑1 (MQ=255)                              gTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATGCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTgg                                     <  1:237358/76‑1 (MQ=255)                               tACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGt                                    <  1:4197927/76‑1 (MQ=255)                               tACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGt                                    >  2:5254205/1‑76 (MQ=255)                                 cTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGttt                                  <  2:307083/76‑1 (MQ=255)                                 cTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGttt                                  >  2:4674832/1‑76 (MQ=255)                                 cTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGttt                                  >  2:5427898/1‑76 (MQ=255)                                    aaCGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTtat                               >  2:1673105/1‑76 (MQ=255)                                     aCGTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GCTGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTtata                              <  1:3351068/76‑1 (MQ=255)                                     aCGTTATTTCAGGGGTGGCAAAACTCAAATTAGAT‑GATGCAGTCAGAGGCGTTTTCCAGTACGTCCTGGTTTtata                              <  1:3366526/76‑1 (MQ=255)                                       gTTATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATAtt                            >  1:3032678/1‑76 (MQ=255)                                        ttATTTCAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTc                           <  2:3689776/76‑1 (MQ=255)                                           tttGAGGGGTGGCACATCTCAAGTTCGAT‑GATGGAGTCAGAGGCGTTTTCCCTTACGTGCTGGTTTTATATTCatt                        <  2:342504/76‑1 (MQ=255)                                              cAGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCattatt                     >  1:2431353/1‑76 (MQ=255)                                               aGGGGTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCAttatta                    <  2:3197265/76‑1 (MQ=255)                                                ggggTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAg                   <  1:2128522/76‑1 (MQ=255)                                                 gggTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAgc                  >  2:369865/1‑75 (MQ=255)                                                 gggTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGt                  >  2:4030260/1‑76 (MQ=255)                                                  ggTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTc                 >  1:3023447/1‑76 (MQ=255)                                                  ggTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTc                 >  2:114520/1‑76 (MQ=255)                                                   gTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCg                <  1:2953893/76‑1 (MQ=255)                                                   gTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCg                <  1:351828/76‑1 (MQ=255)                                                   gTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCg                >  1:362313/1‑76 (MQ=255)                                                   gTGGCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCg                >  2:5028336/1‑76 (MQ=255)                                                    ttgCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGt               <  1:856022/74‑1 (MQ=255)                                                     ggCAAAACTCAAATTCGCT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTa              <  2:705937/76‑1 (MQ=255)                                                     ggCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTa              >  2:910520/1‑76 (MQ=255)                                                     ggCAAAACTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTa              >  1:2730095/1‑76 (MQ=255)                                                        aaaaCTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTAttt           <  2:3954664/76‑1 (MQ=255)                                                         aaaCTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTg          <  1:3007476/76‑1 (MQ=255)                                                         aaaCTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTg          <  2:1269095/76‑1 (MQ=255)                                                         aaaCTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTg          <  1:4977826/76‑1 (MQ=255)                                                         aaaCTCAAATTCGAT‑GATGCAGTCAGAGGCGTGTTCCCTGACGTCCTGGTTTTATATTCATTATTAGTCGTATTTg          <  1:4989444/76‑1 (MQ=255)                                                           aCTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTggt        >  2:2946448/1‑74 (MQ=255)                                                           aCTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGtt        <  1:3331710/76‑1 (MQ=255)                                                            cTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGttt       <  1:1448226/76‑1 (MQ=255)                                                            cTCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGttt       <  2:3412889/76‑1 (MQ=255)                                                            cTCAAATTCGAT‑GAGGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGttt       <  2:28927/76‑1 (MQ=255)                                                             tCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGtttt      >  1:3575871/1‑76 (MQ=255)                                                             tCAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGtttt      >  2:681297/1‑76 (MQ=255)                                                              cAAATTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGttttt     <  2:3304742/76‑1 (MQ=255)                                                               aaaTTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTa    >  2:5181685/1‑76 (MQ=255)                                                               aaaTTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTa    <  2:1491920/76‑1 (MQ=255)                                                               aaaTTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTa    >  1:2162119/1‑76 (MQ=255)                                                                 aTTCGAT‑GATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTAtt  >  2:5852067/1‑76 (MQ=255)                                                                        |                                                                      TCGGCTTAATTACACCGCCTATCGGTAATGTACTTAACGTTATTTCAGGGGTGGCAAAACTCAAATTCGATCGATGCAGTCAGAGGCGTTTTCCCTTACGTCCTGGTTTTATATTCATTATTAGTCGTATTTGTTTTTATT  >  NC_000913/3743532‑3743672

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