Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,022,606	G→A	*246* (TAG→TAA)	fliP →	flagellar biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,022,606	0	G	A	100.0%	180.6 / NA	55	*246* (TAG→TAA)	fliP	flagellar biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/28);  total (27/28)

GCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAGAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAGTATTTTGC  >  NC_000913/2022473‑2022730                                                                                                                                      |                                                                                                                             gCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAgg                                                                                                                         >  2:523542/1‑139 (MQ=255) gCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAgg                                                                                                                         >  1:972897/1‑139 (MQ=255) gCGTGTTGATGGCATTGGGGATGATGATGATTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTGTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAgg                                                                                                                         >  2:1081256/1‑139 (MQ=255) gCGTGTTGATGGCATTGGGGATGATGATAGTTCCACCAGCCACCATTGCACTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAgg                                                                                                                         >  2:734714/1‑139 (MQ=255)              aTTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCt                                                                                                            >  2:484858/1‑139 (MQ=255)                  gggATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAAt                                                                                                        <  2:972897/139‑1 (MQ=255)                  gggATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAAt                                                                                                        <  1:7441/139‑1 (MQ=255)                                  ccccAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAGGAGGCAAAATGACACCTGAATCGGTCATGATGATggg                                                                                        >  2:120129/1‑139 (MQ=255)                                   cccAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATg                                                                                                                  <  1:230185/112‑1 (MQ=255)                                   cccAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATg                                                                                                                  >  2:230185/1‑112 (MQ=255)                                            caTTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCg                                                                              >  1:725498/1‑139 (MQ=255)                                                      cccTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTcg                                                                    >  1:1004895/1‑139 (MQ=255)                                                         tttAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCatgatgat                                                                                           >  2:972329/1‑113 (MQ=255)                                                         tttAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCatgatgat                                                                                           <  1:972329/113‑1 (MQ=255)                                                            aaaCTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTgg                                                              <  2:500054/139‑1 (MQ=255)                                                                tGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                   <  2:637106/82‑1 (MQ=255)                                                                tGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                   >  1:637106/1‑82 (MQ=255)                                                                        tttGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGccc                                                  <  1:484858/139‑1 (MQ=255)                                                                          tGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAGGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGTTGCCCCg                                                <  1:120129/139‑1 (MQ=255)                                                                          tGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCg                                                <  1:235049/139‑1 (MQ=255)                                                                          tGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGAAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCg                                                <  1:774111/139‑1 (MQ=255)                                                                             aCTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGc                                                                                                                        >  1:940019/1‑64 (MQ=255)                                                                             aCTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGc                                                                                                                        <  2:940019/64‑1 (MQ=255)                                                                              cTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAt                                            <  2:1004895/139‑1 (MQ=255)                                                                               tggtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAtt                                           <  1:523542/139‑1 (MQ=255)                                                                                ggtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTg                                                     <  2:575031/128‑1 (MQ=255)                                                                                ggtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTg                                                     >  1:575031/1‑128 (MQ=255)                                                                                ggtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAttg                                          >  1:556698/1‑139 (MQ=255)                                                                                 gtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGc                                                                                                                        <  2:807676/60‑1 (MQ=255)                                                                                 gtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGc                                                                                                                        >  1:807676/1‑60 (MQ=255)                                                                                 gtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAttgt                                         <  1:829128/139‑1 (MQ=255)                                                                                  tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTa                                             <  2:509429/134‑1 (MQ=255)                                                                                  tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTa                                             >  1:509429/1‑134 (MQ=255)                                                                                  tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAttgtt                                        >  2:857397/1‑139 (MQ=255)                                                                                  tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAttgtt                                        >  1:364779/1‑139 (MQ=255)                                                                                      tggctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                   <  2:893870/60‑1 (MQ=255)                                                                                      tggctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                   >  1:893870/1‑60 (MQ=255)                                                                                         ctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGaa                                                                                                         <  1:756289/67‑1 (MQ=255)                                                                                         ctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGaa                                                                                                         >  2:756289/1‑67 (MQ=255)                                                                                           ggcAATTGCTGGTCGGTTCGCTGGCGCAGAGTTTTTACAGCTAAAGAGGc                                                                                                                        <  2:630601/50‑1 (MQ=38)                                                                                           ggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGc                                                                                                                        >  1:630601/1‑50 (MQ=255)                                                                                           ggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                   >  1:34341/1‑55 (MQ=255)                                                                                           ggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                   <  2:34341/55‑1 (MQ=255)                                                                                           ggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGtt                               >  1:1049693/1‑139 (MQ=255)                                                                                                ttGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCa                          <  1:857397/139‑1 (MQ=255)                                                                                                 tGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCAc                         <  2:725498/139‑1 (MQ=255)                                                                                                  gCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACg                        <  1:263627/139‑1 (MQ=255)                                                                                                   cTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACgg                       <  2:1049693/139‑1 (MQ=255)                                                                                                           ttCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGtt                               <  2:604206/123‑1 (MQ=255)                                                                                                           ttCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGCCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGtt                               >  1:604206/1‑123 (MQ=255)                                                                                                                   cgcAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAAt                                                                                                        <  2:598008/42‑1 (MQ=255)                                                                                                                   cgcAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAAt                                                                                                        >  1:598008/1‑42 (MQ=255)                                                                                                                     cAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCg                                                                                                      <  1:393055/42‑1 (MQ=255)                                                                                                                     cAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCg                                                                                                      >  2:393055/1‑42 (MQ=255)                                                                                                                        agCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAGTATTTTGc  >  1:351257/1‑139 (MQ=255)                                                                                                                                      |                                                                                                                             GCGTGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAGAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAGTATTTTGC  >  NC_000913/2022473‑2022730

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