Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 982,560 T→A 100% intergenic (‑122/+481) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273982,5600TA93.7% 40.2 / ‑3.1 16intergenic (‑122/+481)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (6/9);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Not supported by required number of reads on both strands.

CAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATT  >  CP009273/982418‑982701
                                                                                                                                              |                                                                                                                                             
cAGCATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGAAGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTc                                                                                                                                         <  1:134057/149‑1 (MQ=255)
cAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGAAACCTGCCCCTGCCCTCCATAAGTTCTGTAAATAAAAATTTACGGAAATATTGATAAGAGTCTGGTTTCTTTATGTTTCTGCAGACCTCTTACCATTc                                                                                                                                         >  2:142915/1‑149 (MQ=255)
cAGAATATCGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTc                                                                                                                                         <  2:370941/149‑1 (MQ=255)
              tcctcatTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCg                                                                                                                           <  1:170558/146‑1 (MQ=255)
                     tatttattACCCGCATGGTTTTTTTTATGACACCTGCCACTGCCGTCCATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGt                                                                                                                    <  1:294080/149‑1 (MQ=255)
                                                    tccTGCCACTGCGGTCAATAAGTGCGGTTAATAAAAATTTAGGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGGCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTc                                                                                     <  1:38578/148‑1 (MQ=255)
                                                      cTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAATTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaa                                                                                   >  2:304253/1‑149 (MQ=255)
                                                                   cAATAAGGTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAatat                                                                      <  1:216753/149‑1 (MQ=255)
                                                                                   tAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTaaaa                                                      <  1:115197/149‑1 (MQ=255)
                                                                                   tAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTaaaa                                                      >  1:226709/1‑149 (MQ=255)
                                                                                      aaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCGGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAttt                                                   >  1:76398/1‑149 (MQ=255)
                                                                                                   tacgaatgaGAGTTTTGGCTCCTTATCTTTTTCGAGCCATCTTACCATTCAAACTAACGACTAGCTACCGTAAAAGTGCCAAACACAAAAAAATATTTATACAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgtg                                      <  2:284254/144‑1 (MQ=255)
                                                                                                          gagaGTTGGGTGTCTTTATGCGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTAAAAAAAAAAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCg                               <  2:454437/149‑1 (MQ=255)
                                                                                                             agTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCa                            >  1:46121/1‑149 (MQ=255)
                                                                                                                    tgtCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTcc                     <  1:316925/149‑1 (MQ=255)
                                                                                                                                       gCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGCGACGTGATGATAtt  >  2:443189/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                             
CAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATT  >  CP009273/982418‑982701

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