Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,505,925	A→G	I20I (ATT→ATC)	LY180_07200 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,505,925	0	A	G	94.4%	53.7 / NA	18	I20I (ATT→ATC)	LY180_07200	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (8/9);  total (9/9)
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

CATTGTGGTTCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAATAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGAATTACGTCGTTGGATTTATCTTCAAAGAATAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCTTCAATTGG  >  CP006584/1505650‑1506026                                                                                                                                                                                                                                                                                    |                                                                                                      cATTGTGGTTATAAATTACAGTTTTTGTGTCCGTGCTATTTTAAGAATGGAGGGTGGAAGGTCATTGGGGATTATCTTGATCAGCTCCCCGGGGTTGCAGGGGAGGGGCAAGATTAAACCAGTTCGTTCGTGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTGGGATTTATCtt                                                                                 <  2:354791/298‑1 (MQ=255)         ttCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTGATCTTCAAAGaac                                                                         >  1:145886/1‑297 (MQ=255)                ttACAGTTTTGTACCAAGGGCTATTTTGAGGATAGAGGATTAAAGGTGATTGGGGATGTTTTGCATCAGCTCCCCGTGGTTGCAGGGGGGTGGTAAGATTAAACCAGTTCGGTAGGGCAGGTTGCGGCTTTTTGCAGATTGCTTAAGTTTTGCAGGGTAGTCTGAGAAATACTGGTGAGAGATTCGGCTGGCAGGAATGCCTGGTGCTCGGTAAACAGCACGTGGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATGTATCTTCAAAGAACAGATcgc                                                                  <  2:232459/298‑1 (MQ=255)                     gTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGcgttcg                                                             <  1:204740/298‑1 (MQ=255)                       ttCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGcgttcgtt                                                           >  2:410417/1‑298 (MQ=255)                        tttGACTCAGGGCTATTTTAAGAAGAGAGGATGAAAGGTCATTGGGGTTTTGCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGGTCGTTCGGGCAGGGTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCgttcgttc                                                          <  2:26521/296‑1 (MQ=255)                          tGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTtctc                                                        <  1:375836/298‑1 (MQ=255)                              tCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTTGCAGGCAGACCGACGCCGAACTCCATAAACCTGACACACGAGCTTGGAATTATCTTCAAAGCAAACATCGCGTTCGTTCCCacaa                                                    >  2:400028/1‑295 (MQ=255)                                       ttttAAGAAAAGAGGATTAAAGGTCATTGTGAATGTGCTGAATCAGCTCCCCTGGGGTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGGGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATa                                           <  2:363182/298‑1 (MQ=255)                                                          aaGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCt                        >  1:371635/1‑298 (MQ=255)                                                             gggtttGGGGATTATCGGAGTCAGCTCCCCTGGGTTGCAGGGGGGCGGCAAGATTAAACCAGTTCGGTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGAt                     <  2:371635/294‑1 (MQ=255)                                                                   ggggATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCATACAGGCGACGGACTACGTCATCCTGGATCACGTCGTTGCATTTATCTTCAAAGAACAGATCCCCGTAGTTCTCATAAACGTCCATCACCAACCACCCCATTTTTTGATTTTCCa               >  2:350728/1‑298 (MQ=255)                                                                    gggATTGTCTGAATCAGCTCCCCTGGGTTGTAGGGGTGGGGAAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCCGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGGCAGGGATGCCTGGTGCCCGGTAAACAGAACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAg              <  2:60282/298‑1 (MQ=255)                                                                     ggATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGCGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTACCCCGCACGTTGTGGCAGGCAGAAAGGCGACGGTATACGTCATCCTGGCTCCCGTCGTTCGATTGTTCGTCAAAGAACAGATCGCGTTCGTTCTTCTACACGTCACTCCCCAAACAACCAATTTTCTGATTTTTCacc             >  2:54691/1‑296 (MQ=255)                                                                     ggATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAgc             >  1:382977/1‑298 (MQ=255)                                                                        ttATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCt          >  1:108561/1‑298 (MQ=255)                                                                         tATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCTCGTTCGTTCTCATAAACGTCCATACCCAACGAACCACTTTTCTGATTTTTTCgcgcat         >  2:172398/1‑296 (MQ=255)                                                                                ataaGCTGCCATGGTTTGCAGGGGAACGGCAAGATTAAACCAGTTCGTTCGGGTAGGTTTCGCCTTTTTCCAGATTGCTTAAGGTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGATTCTGCTGGCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCTTCAATTgc  <  2:373529/295‑2 (MQ=255)                                                                                                                                                                                                                                                                                    |                                                                                                      CATTGTGGTTCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAATAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGAATTACGTCGTTGGATTTATCTTCAAAGAATAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCTTCAATTGG  >  CP006584/1505650‑1506026

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