Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	449,135	C→A	V503L (GTA→TTA)	cyoB ←	cytochrome bo3 ubiquinol oxidase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	449,135	0	C	A	97.6%	115.9 / NA	42	V503L (GTA→TTA)	cyoB	cytochrome bo3 ubiquinol oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (22/19);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAAC  >  NC_000913/449000‑449274                                                                                                                                        |                                                                                                                                            gAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCg                                                                                                                                      >  1:939784/1‑143 (MQ=255)   ggCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCgct                                                                                                                                    <  1:716089/143‑1 (MQ=255)      ggggAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGCGAATACCCAGCGCCATCAGTAATGCCCCgctggc                                                                                                                                 >  1:116331/1‑143 (MQ=255)        ggAAGAGGTGGCCCTCTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATCGAAACGTACATCTGAATACCGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCgctggctg                                                                                                                               <  2:677096/143‑1 (MQ=255)         gAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGTCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGc                                                                                                                              >  2:766592/1‑143 (MQ=255)                    ccACTCCGGCGTACGGCCACCCCACGGGTCGGCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCa                                                                                                                   <  2:604576/143‑1 (MQ=255)                    ccACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCa                                                                                                                   >  1:53198/1‑143 (MQ=255)                              gTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGACTAACGAGGCAGAGAATACCCAGCGCAGTCAGTAATGCACCGCTGGCGGCAAGCCTCAGCATGGtgtgtag                                                                                                         >  1:966155/1‑140 (MQ=255)                               tACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAAc                                                                                                        <  1:220458/143‑1 (MQ=255)                               tACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAACACGAACTGCATCGCGGGCTGCAATCATCAGCATGGTGTGGAAc                                                                                                        >  1:1008180/1‑143 (MQ=255)                                aCGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACt                                                                                                       >  2:711542/1‑143 (MQ=255)                                          cACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATAAGAATAAGGAGGGAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGGACGTCGGGTCGAt                                                                                             >  1:246324/1‑143 (MQ=255)                                           aCGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATc                                                                                            >  1:22955/1‑143 (MQ=255)                                              ggTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATctgc                                                                                         <  1:345233/143‑1 (MQ=255)                                                 caccagccagCACACGGTTCTGGTAGCGGTCGCGAATAGAACCGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTgg                                                                                      <  1:789763/141‑1 (MQ=255)                                                 cGCCGGTCAGGTCCCGGTTACCGGCGCGGTCGCGAATAGAAACGTACAGCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTgg                                                                                      <  2:692122/143‑1 (MQ=255)                                                   ccAGTCCGGTGACGGTTCTGGTCGCGGTCGCGAATACAACCGTCAACCTGAATAGCGAGGCAGGGAATACCCGGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCt                                                                                    <  2:271131/143‑1 (MQ=255)                                                                 gTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAACCATCAGCATGGTGTGGGACTGCGGGACGATCTGCTGGCTCAACACACGGGTCa                                                                      >  2:888231/1‑143 (MQ=255)                                                                 gTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGAACCGCTGGCTGCAATCATCAGCATGGTGTGGAACCGCGGGTCAATCTGTTGGAGCAAACGACGGGTCa                                                                      >  2:919042/1‑143 (MQ=255)                                                                           cggtcgcgAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGaa                                                            <  1:298065/143‑1 (MQ=255)                                                                           cggcggcgAACAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGaa                                                            <  2:189141/138‑1 (MQ=255)                                                                                gcgAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCa                                                       <  1:711542/143‑1 (MQ=255)                                                                                gcgAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACCGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCa                                                       >  2:432980/1‑143 (MQ=255)                                                                                  gAATAGAAACGTACATCTGAATAACGAGGGAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATTATCAGCATGGTGTGGAACTGCGGGTCAATATGCTCGCTCAAACGACGGGACACGCCCGTGAAcccacgc                                                     >  2:915192/1‑139 (MQ=255)                                                                                  gAATAGAAACGTACAACTGCATAGCGAGGCAGAGAATACCGAGCGCAAGCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGCGTGGAACTGCGGGTGAATCTACTGGCTCAGACGAAAGGCGAACCCCATGAAGccaaca                                                     >  2:132516/1‑139 (MQ=255)                                                                                     tAGAAACGTACAGCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCa                                                  <  1:223838/143‑1 (MQ=255)                                                                                       gAAGCGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATa                                                <  1:315920/143‑1 (MQ=255)                                                                                        aaaCGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATAc                                               >  1:1032860/1‑143 (MQ=255)                                                                                             tACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTgg                                          <  1:473651/143‑1 (MQ=255)                                                                                                aTCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTCGCAt                                       >  1:679597/1‑143 (MQ=255)                                                                                                   tGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATaaa                                    >  1:875423/1‑143 (MQ=255)                                                                                                     aaTAACGAGGCAGAGAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAgg                                  >  1:335986/1‑143 (MQ=255)                                                                                                               cAGGGAGTACCCAGCGCAGTCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAAGGACGGGTCATGCCCATGAAGCGCAGCGCATACAGTGGCATAAAGGCAACgaagaa                        <  2:825913/143‑1 (MQ=255)                                                                                                                agagAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGaaggag                       >  1:102352/1‑143 (MQ=255)                                                                                                                agagAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGaagaag                       >  1:124736/1‑143 (MQ=255)                                                                                                                agagAATACCCAGCGCAATCAGTAATGCACCGCTGGCTGCAAGCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGCCCCCAACCACGGGTCTTGCCCAGGAAGCCCAGCGAATACCGTCGCGTAAAGGCAACGgaggag                       >  2:372025/1‑143 (MQ=255)                                                                                                                      tgcccAGCTCCACCAGTAATGCTCCGCTCGCTGCAATCATCACCATGGTGTGGCACTGCGGCTCTATCTGCCTCCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGACGCCgatg                 <  2:864775/141‑1 (MQ=255)                                                                                                                       aCCCAGCGCAATCAGCAATGCACCGCTGGCTGCAATCATCAGCATGGCGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCgatga                <  1:569524/143‑1 (MQ=255)                                                                                                                         ccggcgcAATCAGTAAGGGCCCGGTGGCTGCAATCAGCAGCATGGTGTGGAACGGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGGGGCATAAAGGCAACGAAGAAGCCGATGATc              <  2:894574/140‑1 (MQ=255)                                                                                                                         ccAGCGCAATCAGTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACAGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATACACAAGAAGCACAGCGCACACAGTGGCATAAAGGCAACGACGACGCCGATGATc              >  2:668239/1‑143 (MQ=255)                                                                                                                         ccAGCGCAATCAGGAATGCTCCGCTGGCTGCAATCATCAGCATGGTGCGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATc              <  2:484275/143‑1 (MQ=255)                                                                                                                                     gTAATGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATcgagaaccccaac  >  2:847375/1‑138 (MQ=255)                                                                                                                                        |                                                                                                                                            GAGGCGGGGAAGAGGTTGCCCACTCCAGCGTACGGCCACCCCACGGGTCGCCAGTCAGGTCACGGTTCTGGTCGCGGTCGCGAATAGAAACGTACATCTGAATAACGAGGCAGAGAATACCCAGCGCAATCAGTACTGCACCGCTGGCTGCAATCATCAGCATGGTGTGGAACTGCGGGTCAATCTGCTGGCTCAAACGACGGGTCATGCCCATGAAGCCCAGCGCATACAGTGGCATAAAGGCAACGAAGAAGCCGATGATCCAGAACCAGAAC  >  NC_000913/449000‑449274

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