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%	142.9 / NA	44	*246* (TAG→TAA)	fliP	flagellar biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/19);  total (25/19)

TGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAGAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAGTATTTTGCAGG  >  NC_000913/2022476‑2022733                                                                                                                                   |                                                                                                                                tgtTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCaa                                                                                                                         <  2:404544/139‑1 (MQ=255)                              cccccAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATgg                                                                                            >  2:1368816/1‑139 (MQ=255)                                          aTTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTg                                                                                      <  1:966708/133‑1 (MQ=255)                                          aTTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTg                                                                                      >  2:966708/1‑133 (MQ=255)                                              ctctGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGaa                                                                            >  1:916299/1‑139 (MQ=255)                                                   cccTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGt                                                                                                       >  2:1171517/1‑107 (MQ=255)                                                   cccTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGt                                                                                                       <  1:1171517/107‑1 (MQ=255)                                                   cccTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTcg                                                                       >  1:771141/1‑139 (MQ=255)                                                             tGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACt                                                             >  1:196139/1‑139 (MQ=255)                                                             tGATGCTGTTTGTACTGGAGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACt                                                             >  1:897273/1‑139 (MQ=255)                                                               aTGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      >  1:768979/1‑80 (MQ=255)                                                               aTGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      <  2:768979/80‑1 (MQ=255)                                                               aTGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTgg                                                           <  2:916299/139‑1 (MQ=255)                                                                      ttGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGccc                                                     >  2:1192825/1‑138 (MQ=255)                                                                      ttGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGccc                                                     <  1:1192825/138‑1 (MQ=255)                                                                       tGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCa                                                                                                                          >  2:306086/1‑68 (MQ=255)                                                                       tGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCa                                                                                                                          <  1:306086/68‑1 (MQ=255)                                                                              gtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGc                                                                                                                           <  1:339621/60‑1 (MQ=255)                                                                              gtggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGc                                                                                                                           >  2:339621/1‑60 (MQ=255)                                                                               tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTcgcg                                                                     <  1:481628/113‑1 (MQ=255)                                                                               tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTcgcg                                                                     >  2:481628/1‑113 (MQ=255)                                                                               tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCg                                                   <  1:904629/131‑1 (MQ=255)                                                                               tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCg                                                   >  2:904629/1‑131 (MQ=255)                                                                               tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAttgtt                                           >  1:1542570/1‑139 (MQ=255)                                                                               tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAttgtt                                           >  1:125343/1‑139 (MQ=255)                                                                               tggATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTAttgtt                                           >  2:300245/1‑139 (MQ=255)                                                                                   tggctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTa                                       >  1:633038/1‑139 (MQ=255)                                                                                    ggctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGccc                                                     <  1:1298967/124‑1 (MQ=255)                                                                                    ggctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGccc                                                     >  2:1298967/1‑124 (MQ=255)                                                                                      ctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      >  1:136800/1‑57 (MQ=255)                                                                                      ctggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      <  2:136800/57‑1 (MQ=255)                                                                                        ggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      >  1:960638/1‑55 (MQ=255)                                                                                        ggcAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      <  2:960638/55‑1 (MQ=255)                                                                                                   gTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      >  2:316958/1‑44 (MQ=255)                                                                                                   gTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAAt                                                                                                                      <  1:316958/44‑1 (MQ=255)                                                                                                             tGGCGCAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAg             <  2:897273/139‑1 (MQ=255)                                                                                                                  cAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCAc                            <  2:583986/119‑1 (MQ=255)                                                                                                                  cAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCAc                            >  1:583986/1‑119 (MQ=255)                                                                                                                  cAGAGCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAGTAttt        >  2:861939/1‑139 (MQ=255)                                                                                                                     agCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGCTGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCAc                            <  1:940126/116‑1 (MQ=255)                                                                                                                     agCTTTTACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCAc                            >  2:940126/1‑116 (MQ=255)                                                                                                                        ttttACAGCTAAAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAGTATTTTGCAgg  <  1:861939/139‑1 (MQ=255)                                                                                                                               gCTAAAGAGGCAAAATGACACCTGAATCGGTCatgatg                                                                                                <  1:1184536/38‑1 (MQ=255)                                                                                                                               gCTAAAGAGGCAAAATGACACCTGAATCGGTCatgatg                                                                                                >  2:1184536/1‑38 (MQ=255)                                                                                                                                   |                                                                                                                                TGTTGATGGCATTGGGGATGATGATGGTTCCCCCAGCCACCATTGCTCTGCCCTTTAAACTGATGCTGTTTGTACTGGTGGATGGCTGGCAATTGCTGGTCGGTTCGCTGGCGCAGAGCTTTTACAGCTAGAGAGGCAAAATGACACCTGAATCGGTCATGATGATGGGGACTGAAGCGATGAAAGTCGCGCTGGCACTGGCTGCCCCGCTATTGTTGGTAGCGTTGGTCACGGGCCTTATCATCAGTATTTTGCAGG  >  NC_000913/2022476‑2022733

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