Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	459,705	Δ1 bp	coding (137/276 nt)	LY180_02175 ←	transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	459,702	0	C	.	91.3%	84.7 / ‑1.3	23	coding (140/276 nt)	LY180_02175	transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (12/9);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCGCCCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGCTAATCGGCAGG  >  CP006584/459441‑459991                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  gCACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCa                                                                                                                                                                                                                                                              >  1:150541/1‑298 (MQ=255)                                    tCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCAt                                                                                                                                                                                                                           >  1:189622/1‑298 (MQ=255)                                          aTTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGCACCGTTGGAGTATGGCATGGCATTCGGCAACCCCCTCCACCCCCCCTTCCAGGGCCTCAATCt                                                                                                                                                                                                                     >  2:165999/1‑298 (MQ=255)                                             tCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGcc                                                                                                                                                                                                                  <  1:165999/298‑1 (MQ=255)                                                            ggAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCGGAGCATCAATCTGCCCTCGAATCCGACaaa                                                                                                                                                                                                   >  1:282228/1‑298 (MQ=255)                                                                                          aTCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCg                                                                                                                                                                     >  1:83739/1‑298 (MQ=255)                                                                                             tttctcttttCCTCAAGTAT‑GTAATAGATTCAGCGACTTAAGATAGCGTCGTACCAGTTCAAGCGTGTCGTCTACGTATTCGCTGCCTGCGCGGCGGTAGCAGGCATGTCGGGCAAGCCTGTCCCGGATACGGCGTTCACGCACTTCTTGCATCTGCGCATTAGCCGCGCCCCGAACGGCAGCGGTCTGTTGGAGTATGGCACAGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCGCTTCCGGa                                                                                                                                                                   <  2:28720/290‑1 (MQ=255)                                                                                                   tctTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTgg                                                                                                                                                            <  1:128738/298‑1 (MQ=255)                                                                                                    ctTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTgga                                                                                                                                                           <  1:192841/298‑2 (MQ=255)                                                                                                               atGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGGAAGGACCTTTTTCTTCTCTTTCGGAGTACTGGGCATTTCGCAc                                                                                                                                                 >  1:282592/1‑298 (MQ=255)                                                                                                                   tAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCAt                                                                                                                                                                                                                                                      <  2:79606/193‑1 (MQ=255)                                                                                                                   tAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCAt                                                                                                                                                                                                                                                      >  1:79606/1‑193 (MQ=255)                                                                                                                          ttCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTg                                                                                                                                      <  1:10575/298‑1 (MQ=255)                                                                                                                                     aaGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAAtata                                                                                                                           >  1:154592/1‑298 (MQ=255)                                                                                                                                                                                  gACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGta                                                                              <  1:7706/298‑1 (MQ=255)                                                                                                                                                                                    cTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGtata                                                                            >  1:252626/1‑298 (MQ=255)                                                                                                                                                                                               tAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATCCGGCAGCACCCTCCCGCGCCCCTTCCCGCGCATCAATCTGCCCCCGAATACCCCGCACTCGCGTACGGACCTTTTTTCTCCCTTCCGTAGTACTGGGCACTAAGCCGATCCACATACGCACACAATATACAACAGGCGCACATGCAATTTGTCAATCAATACCCCCACACAGTTTacaacaatagaa                                                                 >  2:166514/1‑292 (MQ=255)                                                                                                                                                                                                          tCGGTCAAGCGTTCCCCCGTAGTGCATTTCAAGCATTTTTGTCATTAGCCCTTTGTCCGC‑CCCCGGACTGCAGCGATCTGTTGTGGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGGTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTAGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTa                                                      <  2:150541/298‑1 (MQ=255)                                                                                                                                                                                                          tCGGTCAAACGTTTCCCGGATATGGCTTTCAAGAACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTa                                                      <  1:122004/298‑1 (MQ=255)                                                                                                                                                                                                           cGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTAt                                                     <  1:242737/298‑1 (MQ=255)                                                                                                                                                                                                                      ttCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCCTACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAAtgt                                          >  1:305906/1‑298 (MQ=255)                                                                                                                                                                                                                                                  gCCATCAGCCCGTTAGCCGCG‑CCCGAACGGCAGCGATCTGTTGCAGTATGGCACGGCATTCGGCATCACCCTCCAGCTACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAggggggg              <  1:167611/298‑1 (MQ=255)                                                                                                                                                                                                                                                              ttAGCCGCG‑CCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGCTAATCGGCAgg  >  1:187446/1‑298 (MQ=255)                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  GCACTTCACCTTTTTTCGGTGGTGCAACGTCAATTTCAACGATTTCCAGCGGTTTACCGGGAGCAAATGCAACGGCAGCACGTGATTTCATCTCTCGCTCTTCCTCAA‑TATGGTAATAGATTCAGCTATTTAAGATAGGCTCGAACCAGTTCAATAGTGTCGTCAACGGATTGGCTGACTTCGCGGCTGTAGCAGTCATTTCGGTCAAACGTTTCCCGGATATGGCTTTCAAGCACTTCTGCCATCAGCCCGTTAGCCGCGCCCCGAACGGCAGCGATCTGTTGGAGTATGGCACGGCATTCGGCATCACCCTCCAGCGACCGTTCCAGAGCATCAATCTGCCCTCGAATACGACGAACTCGAGTAAGGACCTTTTTCTTCTCTTCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCACGCAATATACTATAGGGGGGTATGCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGCTAATCGGCAGG  >  CP006584/459441‑459991

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