Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	510,152	Δ1 bp	coding (164/459 nt)	ybbJ ←	NfeD family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	510,151	0	C	.	94.4%	137.8 / 0.5	36	coding (165/459 nt)	ybbJ	NfeD family protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (16/18);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTC  >  NZ_CP009273/510064‑510233                                                                                        |                                                                                   gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                                  >  2:312580/1‑88 (MQ=255) gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                                  >  1:2176901/1‑88 (MQ=255) gtgCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcac                                                                                  >  2:1783551/1‑88 (MQ=255)     ttttGGCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTccc                                                                             <  1:290726/90‑1 (MQ=255)      tttGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCa                                                                            >  2:2575791/1‑90 (MQ=255)             gcgcACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGc                                                                     >  2:3014979/1‑90 (MQ=255)             gcgcACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGc                                                                     >  1:1538915/1‑90 (MQ=255)             gcgcACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGc                                                                     >  1:1696251/1‑90 (MQ=255)                  cccGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCAc                                                                <  2:1538915/90‑1 (MQ=255)                  cccGCCGCGAAAACCATTTCCCCCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCAc                                                                <  1:1324296/90‑1 (MQ=255)                      ccgcGACAACCATTTCCACCACAGCCAGG‑GGCGAGCCGCGTCAGGAAGGCAAACATTAACCCCTG‑CACTCCCCACCCCGCGGCAcccgcc                                                           >  2:2417557/1‑90 (MQ=255)                         cgACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCAccagccag                                                         <  1:1140217/90‑1 (MQ=255)                         cgACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCAccagccag                                                         <  1:546529/90‑1 (MQ=255)                          gACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGa                                                        <  1:598663/90‑1 (MQ=255)                                   ttCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CTCTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                               <  1:1476215/90‑1 (MQ=255)                                   ttCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                               >  2:455785/1‑90 (MQ=255)                                   ttCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGc                                               <  1:2868032/90‑1 (MQ=255)                                    tCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGcc                                              <  1:443000/90‑1 (MQ=255)                                        ccacAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTa                                          >  2:2494835/1‑90 (MQ=255)                                         cacaGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTaa                                         >  1:603701/1‑90 (MQ=255)                                         cacaGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTaa                                         >  2:2688573/1‑90 (MQ=255)                                           caGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATc                                       <  1:1783551/90‑1 (MQ=255)                                              ccAGGCGGCGAGCAGCGTCAGGATGGCAAACATTCCCCCTTGCC‑CTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                    <  1:1388739/90‑1 (MQ=255)                                              ccAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACGCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                    <  2:1801090/90‑1 (MQ=255)                                              ccAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                    <  1:286607/90‑1 (MQ=255)                                                    ggcgAGCAGCGTCAGGATGGCAAACATTCCCCCTTG‑CACCCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGcc                              <  1:5005/90‑1 (MQ=255)                                                        agcagcGTCAGGATGGCAAACATTCCCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgc                          <  1:2575791/90‑1 (MQ=255)                                                         gcagcGTCAGGATGGCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCAcgcc                         <  1:2688573/90‑1 (MQ=255)                                                                      ggCAACCATTCCCCCTTG‑CACCCCCACCCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa            <  1:2626491/90‑1 (MQ=255)                                                                      ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa            >  2:2054287/1‑90 (MQ=255)                                                                      ggCAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACaaa            >  2:2994865/1‑90 (MQ=255)                                                                        cAAACATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATa          >  1:1428494/1‑90 (MQ=255)                                                                         aaaCATTCCCCCTTGCC‑CTCCCAACCCAGCGGCCCCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATaa         <  2:1352434/90‑1 (MQ=255)                                                                           aCATTCCCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAAcc       <  2:1628864/90‑1 (MQ=255)                                                                            cATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCa      >  2:1218626/1‑90 (MQ=255)                                                                            cATTACCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCa      >  2:677787/1‑90 (MQ=255)                                                                                aCCCCTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTc  >  1:982911/1‑90 (MQ=255)                                                                                   ccTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                    <  2:2734654/53‑1 (MQ=255)                                                                                   ccTTG‑CACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCAcc                                    >  1:2734654/1‑53 (MQ=255)                                                                                        |                                                                                   GTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCAGCGGCACCAGCCAGACCACCAGGCCGGTAATCACCGCTGCCACGCCGCTCCACAACAAATAACCATTTC  >  NZ_CP009273/510064‑510233

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