Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,058,823	G→A	W139* (TGG→TGA)	bipA →	50S ribosomal subunit assembly factor BipA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,058,823	0	G	A	100.0%	52.9 / NA	19	W139* (TGG→TGA)	bipA	50S ribosomal subunit assembly factor BipA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/12);  total (7/12)

TACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCC‑ATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCACGTTCCTGCGCCGGACGTTGACCTTGACGGTCCGTTCCAGATGCAGATTTCTCAGCTCGATTACAACAGC  >  NC_000913/4058608‑4059054                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        tACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTTGGTGGTGAAGTTGAACGTGTAATGTCC‑ATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcc                                                                                                                                                                                                       <  1:25624/250‑1 (MQ=255) tACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACGTCGGCGGTGAAACTGAACGTGTACTGTTC‑ATGGAAGACTCAGTATTGCTGGAGGTTGACGAATTTGAAGGCCCGCTGCCGCCACCGCGCTTCGTAAACAATAAAGCGTATGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTAACCGCCCTGGCGCGAGTTCTGATCGAGTTGTGGTTTAGGGATTCGATCTGTTCGTTAACCt                                                                                                                                                                                                      >  1:434035/1‑251 (MQ=255)  aCCGTATCAACATCGTTGATTCCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCC‑ATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCGTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTAGCAACAAATTTGACCGCCCTTGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcat                                                                                                                                                                                                      <  1:414657/250‑3 (MQ=255)       aTCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCC‑ATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGAcgc                                                                                                                                                                                                <  1:157177/251‑1 (MQ=255)                                        ngcnngggtnAAGGTGAACGTGTAATGTCCCATG‑NCGACCGAGTGCTTCGGGTGGTTGTAGCATTGNANGGCCCGCAGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACNCCCTGAAGCCNATTNNTGTTNTNANCAAAGTTGACCGCCCTGGCGCGCGTCCTGNTTGAGTTGTGGATCAGGTATTCGNNCTNNNCGNNAACCt                                                                                                                                                                                                      <  2:339115/206‑1 (MQ=255)                                                                              aCTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACCCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTtg                                                                                                                                                                                                                                     >  1:194448/1‑144 (MQ=255)                                                                                           ggtggtTGACGCATTTGACGGCCCGATGCCGCAAACGCTCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCTATTGTTGTTATCAACAAAGTTGACCGCCCTGGCTCGCGTCCTGATTGATTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCgn                                                                                                             <  1:265333/251‑2 (MQ=255)                                                                                                                  cGATGCCGCAATCGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAATCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGTCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTAGCGCGGGTCTGGACCACGAAGATATGGCGGa                                                                                                            <  1:291012/229‑1 (MQ=255)                                                                                                                       ccGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTGCGCTTCTGCGCTGAACGGTATCGCGGTTCTGGACCACGAAGATATGGCGGAAGACATGACCGCGCTGTACCAGGCGAt                                                                                 <  1:159569/251‑1 (MQ=255)                                                                                                                       ccGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATCGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGCTTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcc                                                                                                                                                                                                       >  1:229670/1‑133 (MQ=255)                                                                                                                         gCAAACGCGCTTCGTAACCAAAAAAGCGTTGGCTTACGGCCTGAAGTTGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATAAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGagcagc                                                                                                                                                                                  <  1:230812/152‑1 (MQ=255)                                                                                                                          cAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGCTTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAAcc                                                                                                                                                                                                       >  1:208054/1‑130 (MQ=255)                                                                                                                                  cTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCACg                                                                      >  1:372247/1‑251 (MQ=255)                                                                                                                                                            aNGGNCTGAANCCGANTGTTGTTATCAACAANGTTGACCGCCCTGGCGCGCNNCNTGATTNAGTtg                                                                                                                                                                                                                                     <  2:194448/66‑1 (MQ=255)                                                                                                                                                              ggCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGCTCAGGTATTAGATCTGTTCGTTAACCTCGACGCGACCGACGAGCCGCTGGACTTCACGATCg                                                                                                                                                                   >  1:289722/1‑130 (MQ=255)                                                                                                                                                               gCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCACGTTCCTGCGCCGGa                                                         >  1:174173/1‑235 (MQ=255)                                                                                                                                                                                         aaaGTTGACCGCCCTGGCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTc                                                                                                                                  <  1:443821/136‑1 (MQ=255)                                                                                                                                                                                                        ggCGCGCGTCCTGATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCACGTTCCTGCGCCGGACGTTGACCTTGACGGTCCGTTCCAGATGCAGATTTCTCAGCTCGATTACAACAGc  <  1:344793/249‑1 (MQ=255)                                                                                                                                                                                                                    gATTGAGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAng                                                                                     <  1:349709/154‑1 (MQ=255)                                                                                                                                                                                                                         |                                                                                                                                                                                                                                        TACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGTCC‑ATGGTAGACTCAGTGCTGCTGGTGGTTGACGCATTTGACGGCCCGATGCCGCAAACGCGCTTCGTAACCAAAAAAGCGTTTGCTTACGGCCTGAAGCCGATTGTTGTTATCAACAAAGTTGACCGCCCTGGCGCGCGTCCTGATTGGGTTGTGGATCAGGTATTCGATCTGTTCGTTAACCTCGACGCGACCGACGAGCAGCTGGACTTCCCGATCGTTTACGCTTCTGCGCTGAACGGTATCGCGGGTCTGGACCACGAAGATATGGCGGAAGACATGACCCCGCTGTACCAGGCGATTGTTGACCACGTTCCTGCGCCGGACGTTGACCTTGACGGTCCGTTCCAGATGCAGATTTCTCAGCTCGATTACAACAGC  >  NC_000913/4058608‑4059054

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