Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 1,505,880 T→C L35L (CTA→CTG LY180_07200 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065841,505,8800TC92.9% 41.3 / ‑6.5 14L35L (CTA→CTGLY180_07200hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (9/4);  total (10/4)
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

CCGCCATCAGCAGGCTTAGCGCAACGAGCGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAATAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGAATTACGTCGTTGGATTTATCTTCAAAGAATAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCTTCAATTGGCGTAACGGC  >  CP006584/1505607‑1506035
                                                                                                                                                                                                                                                                                 |                                                                                                                                                           
ccGCCATCAGCAGGCTTAGCGCAACGAGCGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAgg                                                                                                                                     <  1:225681/298‑1 (MQ=255)
                  gcgaAACGGGTGTGGTTACTTTCTTCATTGTGTTTGTCAATTACAGTTTCGGACGCAGGGTTATTTTAAGAATAGAGGATGAAAGGGTATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGTGAGCGGCAAGATTAAACAAGTTGGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATcc                                                                                                                   <  2:336600/298‑1 (MQ=255)
                    gcAACGAGCGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGAAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTg                                                                                                                 >  1:89403/1‑298 (MQ=255)
                          gtggcggcgATTTTCTGATTTGTGGTTTTCAGTTTAAGTTGCGGACGGAGGGCTATTTTAAGAGTAGATGATTGAAGGTCATTGTGGAGTATCTGAATCAGATCCCCTGGGTTTCCGGGGAGCGGCAAGATTAAACCAGTTCGTTCGAGCAGGTTTCGCCTGTTTCCATATTGCTTAAGTTGTGCAGCGTTGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTTCCCGGTAAACAGCACGTTGTGGCATGCAGACAGGCGACGGAATACGTCAGCCTGGATCAc                                                                                                           <  2:319146/295‑1 (MQ=255)
                                     tttCTTCATTGTGGTTCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGAtt                                                                                                >  1:255872/1‑298 (MQ=255)
                                                                       cTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATa                                                              >  1:157384/1‑298 (MQ=255)
                                                                                     aaGACTAGAGGATGAAAGGTCATTGGGGATTATCTGAATCACCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAGCCAGTTCGTTCGGGCAGGTTTCGCCTTTGTCCAGATTGCTCAAGTTTTGCAGCGTAGTCTTAGAACTACTCGTCAGAGCTTCTGCTGTCAGGACTGCATGGTGCCCGGTAACTAGCAAGCTGTGGACGGCAGACAGGCGAGGGCATACGTCCTACGGGATCAGAGCGTTGGACTGATCTTAACAGAACAGATCGCTGTCGGTCGCACACACGTCAGTAccga                                                >  2:105404/1‑296 (MQ=255)
                                                                                         aTAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGa                                            >  1:221946/1‑298 (MQ=255)
                                                                                         aTAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGa                                            >  1:45424/1‑298 (MQ=255)
                                                                                                   aaaGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTc                                  <  1:73444/298‑1 (MQ=255)
                                                                                                                        gAATCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCCGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGACTACGTCATCCTGGATCAAGTCGTTGGATTTACCTTCAACGAAAATAACGCTTTCGTTCTCATCCACCTCACTCCCACACGACCCATATTTCTACTTTTTCCGCGCTCCCAtt             >  2:287954/1‑298 (MQ=255)
                                                                                                                          aTCAGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCTTCAATTgc           >  1:168358/1‑297 (MQ=255)
                                                                                                                             aGCTCCCCTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTGTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACCGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTACAAAGAACAAATCGCGTTCGGTCTCAGACACGTCCATACCCAACGAACCAATTTACTCATTTTTCAgcacctcccttgttcc        >  2:280713/1‑284 (MQ=255)
                                                                                                                                   ccTGGGTTGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAACAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGGATCACGTCGTTGGATTTATCTTCAAAGAACAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCTTCAATTgctgcctgaga  >  1:179397/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                           
CCGCCATCAGCAGGCTTAGCGCAACGAGCGCGGCTACTTTCTTCATTGTGGTTCTCAATTACAGTTTCTGACTCAGGGCTATTTTAAGAATAGAGGATGAAAGGTCATTGGGGATTATCTGAATCAGCTCCCCTGGAATGCAGGGGAGCGGCAAGATTAAACCAGTTCGTTCGGGCAGGTTTCGCCTTTTTCCAGATTGCTTAAGTTTTGCAGCGTAGTCTGAGAAATACTGGTCAGAGCTTCTGCTGTCAGGAATGCCTGGTGCCCGGTAAATAGCACGTTGTGGCAGGCAGACAGGCGACGGAATACGTCATCCTGAATTACGTCGTTGGATTTATCTTCAAAGAATAGATCGCGTTCGTTCTCATACACGTCCATACCCAACGAACCAATTTTCTGATTTTTCAGCGCTTCAATTGGCGTAACGGC  >  CP006584/1505607‑1506035

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