Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,735,910	C→A	37.3%	V66L (GTA→TTA)	U712_19255 ←	Quinol oxidase subunit 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,735,910	0	C	A	37.3%	42.4 / 5.9	40	V66L (GTA→TTA)	U712_19255	Quinol oxidase subunit 2
Reads supporting (aligned to +/- strand):  ref base C (9/12);  new base A (12/2);  total (26/14)
Fisher's exact test for biased strand distribution p-value = 1.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCTGTACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGTAAATAGAACAAAGACGACCCCGACGATAAACAGCATAAATCCGATAGATAACAGGATTAAATCACTTTGCTGTTCAGCTACAGGTCCTTTCGGATCTAGTACT  >  CP006881/3735773‑3736025                                                                                                                                           |                                                                                                                    cGTCTGTACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTAc                                                                                                                     <  6:44246/138‑1 (MQ=255)   tCTGTACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGGATTTTACtt                                                                                                                   >  3:75026/1‑137 (MQ=255)     tGTACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGGTTTTTTATaaa                                                                                                                 >  6:260179/1‑138 (MQ=255)       tACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTTCTTCTTAAAACGTGTTACGGTGAATTTCTGGGTTATAAGATCCTTTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAatta                                                                                                              <  6:350920/139‑1 (MQ=255)       tACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTTATAAaat                                                                                                               >  5:152961/1‑138 (MQ=255)       tACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTaaaa                                                                                                                <  8:98340/137‑1 (MQ=255)       tACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAaata                                                                                                              <  5:260179/139‑1 (MQ=255)       tACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAaata                                                                                                              <  1:155009/139‑1 (MQ=255)          tGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTAGC‑GGGTTTTTTATTAaataatt                                                                                                           >  1:207730/1‑138 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTCCCTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTAGCTGGTTTTTTATTAAATAAATg                                                                                                           >  2:68984/1‑138 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAATCGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATg                                                                                                           <  7:114225/137‑1 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTTTTTTTATAAaaaaattgg                                                                                                         >  1:130374/1‑136 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTTCTTAAATAATgg                                                                                                          >  7:211890/1‑138 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATg                                                                                                           >  8:367650/1‑137 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATg                                                                                                           <  7:367650/137‑1 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGt                                                                                                         <  6:266029/139‑1 (MQ=255)            ttGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACATCTAAAAACGTGTTACCGTGAAATACTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTAAATTAAATAAATAAAGGt                                                                                                         >  4:168576/1‑139 (MQ=255)              ggCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTTCTTAAATAATTGTa                                                                                                        >  8:62794/1‑138 (MQ=255)              ggCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTAc                                                                                                                     >  6:337341/1‑125 (MQ=255)              ggCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTAc                                                                                                                     <  5:337341/125‑1 (MQ=255)               gCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGTaaa                                                                                                      <  1:109659/139‑1 (MQ=255)               gCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGTaaa                                                                                                      <  5:324179/139‑1 (MQ=255)               gCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGTaaa                                                                                                      <  1:365157/139‑1 (MQ=255)               gCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGTaaa                                                                                                      <  3:300986/139‑1 (MQ=255)               gCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGTaaa                                                                                                      <  3:273694/139‑1 (MQ=255)                 acaGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGCGGGTGTTTTATAAAAAAATTGTTaaaa                                                                                                     >  1:213405/1‑138 (MQ=255)                 acaGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTTTTTTTCTAAAATAATGGTAAAt                                                                                                     >  8:264386/1‑138 (MQ=255)                 acaGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTTTTTTAATTAAATAATGGTAAATa                                                                                                    >  7:224922/1‑139 (MQ=255)                          gCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑AGTATTTTACTAAAATAATGGTAAATAGAACAAAga                                                                                           >  2:249153/1‑139 (MQ=255)                                   agtAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTTTTTTTGGTAAATAAATGTTTATATAAAAAATACGGCCCcgc                                                                                  >  6:98114/3‑138 (MQ=255)                                   aTTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGGTTTTTTCTATAATAATGGTTAATTGTACAAACCCGCCCCcgc                                                                                  >  6:300661/1‑138 (MQ=255)                                                                ttCTAAAGACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGGTTTTTTATAAAATAATTTTAAATATAACAAAGACGACCCCGACGATAAACAGCATAAATACGATATATaa                                                      >  7:275356/1‑138 (MQ=255)                                                                ttCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTTTTTTTATAAAATAATTTTTAAAATAAAAAAAACCACCCCGACGATTAAATGCATAAAGCCgatataata                                                      >  5:341279/1‑133 (MQ=255)                                                                ttCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGGTTTTTTTTAAAATAAATGTTAATATAAAAAAAACCCACCCGCACATAAAAAGCAAAAATACGAAATATaa                                                      >  3:64384/1‑138 (MQ=255)                                                                ttCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGGATTTTTATAAAATAATTTTAAATATAAAAAACAACACCCCCCCGCTAAACATCATAAATACCATATATaa                                                      >  1:35328/1‑138 (MQ=255)                                                                ttCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTATACGGTCGCGGGTATTTAAATAAAAAAAGGGT‑AATAAAACAAAGAAGCACCCGCCCAATAACAACATAAAAGAGATAGGTaa                                                      >  3:352538/1‑138 (MQ=255)                                                                taaaaaaacCGTGTCCCTGTCAATTTGTGGTTTATACTACCAATTATATTTCCTTTGCCGGGCGC‑TTTTTTTAAAAAAAATAATGAAAAAAACAACACCCACGACCCCGACGATAAACAGCATAAATCCGATAGATAAc                                                     <  4:344223/135‑1 (MQ=255)                                                                         cGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTTCTAAAATAAAGGTTAATATAACAAAGACGCCCCCCCCGCTAAACAGCATAAATAAGATAAATTACAGGATTaa                                             >  7:91453/1‑138 (MQ=255)                                                                            gtTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAAAAATGGTAAATAGAAAAAAGCAGCCCCCGCCGATAAACAGCAAAAATACGATAGATAACAGGATTAAATCa                                         >  4:351209/1‑139 (MQ=255)                                                                                      tttttTGGTTTATAATATCATTTATATTTCCCTTTCGGGTCTC‑TTTTTATAAAAAAAATTATTTTAAATACAACACCCACGACCCCGACGATAAACAGCATAAATCCGATAGATAACAGGATTAAATCACTTTGCTGtt                               <  6:152961/138‑1 (MQ=255)                                                                                          cTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGGTTTTTTATAAAATAATTGTTAATATAACAAAGACCACCCCGCCGATAAAAAGCAAAAATACGATAGAAAAAATGATTAAAAAACATTGGTGTTGAGc                           >  3:217941/1‑139 (MQ=255)                                                                                                                  tGCCTTTACGGTCGC‑GGGTTTTTTATAAAATAATTGTTAATAAAAAAAAAACGACCCCCCCGAAAAAAAGCAGAAATACGAAAGTATACAAGTTTTAATAACTTTTGTGTTCAGGATCAGGGCGTTTGGGAtggatta    >  2:118288/1‑132 (MQ=255)                                                                                                                   gCCTTTACGGTCGC‑GGTTTTTTTGTAAAATAATTGTTAATAGAACAAAGACGACCCCGCCGCTAAACAGCATAAATCCGATAGATAACAGGATTAAATCACTTTGGTGGTCAGCTTCAGGTCCTTTCGGGTCTAGTCCt  >  3:32655/1‑139 (MQ=255)                                                                                                                                           |                                                                                                                    CGTCTGTACTGTTGGCACAGAAAGTGCAATGACGATTAAAATTGGGATCACTGTCCAGACTACTTCTAAAAACGTGTTACCGTGAATTTCTGGGTTATAAGATCCATTATCTTTGCCTTTACGGTCGC‑GGTATTTTACTAAAATAATGGTAAATAGAACAAAGACGACCCCGACGATAAACAGCATAAATCCGATAGATAACAGGATTAAATCACTTTGCTGTTCAGCTACAGGTCCTTTCGGATCTAGTACT  >  CP006881/3735773‑3736025

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