Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 478,502 C→A 100% D342E (GAC→GAA USA300HOU_0457 → NADH dehydrogenase (ubiquinone)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730478,5020CA97.7% 63.8 / ‑6.2 22D342E (GAC→GAAUSA300HOU_0457NADH dehydrogenase (ubiquinone)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (7/14);  minor base T (1/0);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGTTCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTATCTTGGAATCA  >  CP000730/478364‑478639
                                                                                                                                          |                                                                                                                                         
aGTTCAAAGTTCATTGGGTGTATACTCAGCAGCGATTATTCATTTAATATTGCGCGGTATTTGTAAAGCAACATTATTTTTACAACCAGGTTCTACCGGGAATCGACTCAATATTCCAAAACAAGCATCTGCTAAAGAAGc                                                                                                                                         <  2:124835/141‑1 (MQ=255)
  ttCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCtt                                                                                                                                       <  2:34663/141‑1 (MQ=255)
                   gtATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATg                                                                                                                      <  1:22282/141‑1 (MQ=255)
                          cagcagCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGCCgtg                                                                                                               >  2:1256/1‑141 (MQ=255)
                               gcGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGGGTATTa                                                                                                          >  2:90996/1‑141 (MQ=255)
                                                atTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCAtt                                                                                         <  1:19813/141‑1 (MQ=255)
                                                   gCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGAATCCACGGAATATGCTTATGGCTGGATTGTCATGGGACGGGTATTAGCTATTATCGTGGCATTTCt                                                                                      >  1:181421/1‑141 (MQ=255)
                                                     cgtgTATTTTTAAAGCAACATTATGGTTACAATCAGGTTCTATCGTGACGCGACTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTAt                                                                                    <  1:75625/138‑1 (MQ=255)
                                                     aaggTAATTTAGAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTAt                                                                                    <  2:237658/139‑1 (MQ=255)
                                                      cGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTAtt                                                                                   <  1:216793/141‑1 (MQ=255)
                                                                    acaCATTTTTTTTACACTCAGGTTCTATCGTGAAGTGATTCAATATTCCAAAACAAGCATCTGGTAAAGAAGCTTATGGCTGGACTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTg                                                                     <  2:19262/139‑1 (MQ=255)
                                                                      aCATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTACTTTGGATGAGTAGTGAc                                                                   >  1:66642/1‑141 (MQ=255)
                                                                             ttttACAATCAGGTTCTATCTTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTg                                                            >  2:241207/1‑141 (MQ=255)
                                                                                   aaTCAGGTTCTATCGTAAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTGGCTATTATCGTGGCATCTCTATTTTGGATGAGTAGTGACAGAAGTGCATATg                                                      <  2:176059/141‑1 (MQ=255)
                                                                                            cTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTaa                                             >  1:91045/1‑141 (MQ=255)
                                                                                                         aTTCTATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAATTGCTAAGTGCACTCATTCt                                <  1:139518/141‑1 (MQ=255)
                                                                                                         aTTCAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTTGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCt                                <  1:144789/141‑1 (MQ=255)
                                                                                                            cAATATTCCAAAACAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGGGCACTCATTCTAgg                             >  1:31892/1‑141 (MQ=255)
                                                                                                                      aaaCAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTa                   >  2:201180/1‑141 (MQ=255)
                                                                                                                       aaCAAGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTAc                  <  1:241207/141‑1 (MQ=255)
                                                                                                                          aaGCATCTGCTAAAGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTg               <  2:130049/141‑1 (MQ=255)
                                                                                                                                       aGAAGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCGAGGATGGTCATTACTTGTATCTTGGAATCa  <  2:144083/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AGTTCAATGTGCATTGGGTGTATATTCAGCAGCGATTATTCATTTAATATTGCACGGTATTTTTAAAGCAACATTATTTTTACAATCAGGTTCTATCGTGAAGCGATTCAATATTCCAAAACAAGCATCTGCTAAAGACGCTTATGGCTGGATTGTCATGGGACGTGTATTAGCTATTATCGTGGCATTTCTATTTTGGATGAGTAGTGACAGAAGTGCATATGAAGTGTTAAGTGCACTCATTCTAGGATGGTCATTACTTGTATCTTGGAATCA  >  CP000730/478364‑478639

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: