Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	451,552	A→G	V20A (GTA→GCA)	cyoA ←	cytochrome o ubiquinol oxidase subunit II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	451,552	0	A	G	100.0%	131.2 / NA	34	V20A (GTA→GCA)	cyoA	cytochrome o ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/14);  total (20/14)

GTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAA  >  NC_000913/451404‑451697                                                                                                                                                     |                                                                                                                                                  gTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCa                                                                                                                                                 <  1:279705/151‑1 (MQ=255)   aCTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAg                                                                                                                                                <  1:557416/150‑1 (MQ=255)   aCTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGt                                                                                                                                               <  1:67356/151‑1 (MQ=255)        cAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTg                                                                                                                                          >  2:454564/1‑151 (MQ=255)        cAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTg                                                                                                                                          >  1:233847/1‑151 (MQ=255)         aGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGc                                                                                                                                         <  1:162468/151‑1 (MQ=255)          ggCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCa                                                                                                                                        <  2:284872/151‑1 (MQ=255)                  cAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGa                                                                                                                                <  1:294645/151‑1 (MQ=255)                            cAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATccc                                                                                                                      >  2:341332/1‑151 (MQ=255)                            cAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATcc                                                                                                                       <  1:591862/150‑1 (MQ=255)                               gATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCaaa                                                                                                                   <  1:487309/151‑1 (MQ=255)                                        aaTAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTg                                                                                                          >  1:296606/1‑151 (MQ=255)                                           aaCGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTAt                                                                                                       >  1:137002/1‑151 (MQ=255)                                              gACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTc                                                                                                    >  1:314851/1‑151 (MQ=255)                                                    cAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGt                                                                                              >  1:654644/1‑151 (MQ=255)                                                     aaCATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGtc                                                                                             >  1:675056/1‑151 (MQ=255)                                                      aCATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGtct                                                                                            <  1:302982/151‑1 (MQ=255)                                                            ggCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCAttt                                                                                       >  2:416612/1‑150 (MQ=255)                                                                      cGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCa                                                                             >  1:570366/1‑150 (MQ=255)                                                                            tATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCAc                                                                      >  2:544182/1‑151 (MQ=255)                                                                            tATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCAc                                                                      >  1:453181/1‑151 (MQ=255)                                                                                    aaCGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCtg                                                              >  1:529162/1‑151 (MQ=255)                                                                                     aCGTTGCACCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCtgg                                                             >  1:552751/1‑151 (MQ=255)                                                                                          gCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCg                                                        >  1:261147/1‑151 (MQ=255)                                                                                                 aCCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAg                                                 <  2:293794/151‑1 (MQ=255)                                                                                                         gTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGtgtgtg                                           <  1:565090/149‑1 (MQ=255)                                                                                                          tCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGc                                        >  2:133567/1‑151 (MQ=255)                                                                                                           ccTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGc                                        >  1:1784/1‑150 (MQ=255)                                                                                                                    tCTAACAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACggg                              >  2:304637/1‑151 (MQ=255)                                                                                                                       aaCAGCGCAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGaa                            <  1:178337/150‑1 (MQ=255)                                                                                                                              cAGAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACAt                    >  1:470315/1‑151 (MQ=255)                                                                                                                                gAATTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTa                  <  1:260061/151‑1 (MQ=255)                                                                                                                                 aaTTACAGCCACTGAGCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTAc                 >  2:181440/1‑151 (MQ=255)                                                                                                                                                gCAATGCAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTaa  <  1:448456/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GTACTTCCAGGCGAAACCAACAGCCATCAAGATTGCGGGAATAACGACAATCAACATCAGGCCAAATGCCGTCAGTATCAGTGAACGTTGCTCCAGACCAATCTGTCCTTTGGGATCTAACAGCGCAGAATTACAGCCACTGAGCAATACAGTGCCTGCAAATAATGACAACCATCCCAAACTTTTATTGTATTTCCTGAGTCTCATTTAACGACCTCAATTCCACGGGATCTGGTGGCGTTTAAAGTGTGTGGGCATTTTACGGGAAGGTTACATTACTGTAAACATGATTAA  >  NC_000913/451404‑451697

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