Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 982,560 T→A 23.6% 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,5600TA23.6% 62.8 / 17.8 34intergenic (‑122/+481)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (13/13);  new base A (3/5);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTAC  >  CP009273/982416‑982704
                                                                                                                                                |                                                                                                                                                
gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAt                                                                                                                                              >  2:498724/1‑149 (MQ=255)
 ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAtt                                                                                                                                             >  2:472358/1‑149 (MQ=255)
         ttGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTGCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCa                                                                                                                                           >  1:493500/1‑143 (MQ=255)
         tggcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCCGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTaa                                                                                                                                     <  2:109289/147‑1 (MQ=255)
           gcgcTTCATCATTATTTATTACCCTCATGGTTTTTTTTTTTTTTTTTGCCACTGCCGTCAATAATTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCcaa                                                                                                                                         >  1:99532/1‑143 (MQ=255)
               ttCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTaa                                                                                                                                     <  1:31522/143‑1 (MQ=255)
               ttCATCATGATTTATTACCCTCCTGGTTTTTTTGATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTc                                                                                                                               <  2:28886/149‑1 (MQ=255)
                tcatcaTTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAAGGTTTCg                                                                                                                              >  2:167601/1‑149 (MQ=255)
                      ttatttatTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCg                                                                                                                        >  2:239007/1‑149 (MQ=255)
                                 ccTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCAGCTTTCCATTCAAACTAACGTTTCGCTACCGTGaaa                                                                                                                   <  1:352260/143‑1 (MQ=255)
                                     aTGGTTTTTTTTAGGACACCTGCCACTGCTGTCCATAAGTTCTGTCAAGCAAAATTTCCGGAACTATTCATGCGAGTTTGGTGGCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACCTTTCGCTACCGTGAAAGTGc                                                                                                               <  1:458516/143‑1 (MQ=255)
                                           tttttcATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaa                                                                                                         <  1:366829/143‑1 (MQ=255)
                                            tttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaag                                                                                                        <  1:447787/143‑1 (MQ=255)
                                                       ccTGCCACTGCCGTCAATAAGTTCTGTCAATATAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACGAACGTTTCGCTACCGTGAAAGTGCTACAAAGATATAGATTTGGTTTCa                                                                                       <  2:433928/149‑1 (MQ=255)
                                                        cTGCCACTGCCGTCAATAAGATCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTgg                                                                                            >  1:206557/1‑143 (MQ=255)
                                                             aCTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTAGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTGCACAGATAAAGATTTGGGTTGAAAAAgg                                                                                 >  2:438093/1‑148 (MQ=255)
                                                               tGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaa                                                                                     >  1:206818/1‑143 (MQ=255)
                                                                 ccGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCaaaaa                                                                                   <  1:524557/143‑1 (MQ=255)
                                                                  cGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaaaa                                                                            <  2:547978/149‑1 (MQ=255)
                                                                          aGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAata                                                                          >  1:147476/1‑143 (MQ=255)
                                                                                      aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAg                                                              >  1:411252/1‑143 (MQ=255)
                                                                                       aaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAtt                                                       >  2:14959/1‑149 (MQ=255)
                                                                                          aTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTaa                                                    <  2:395008/149‑1 (MQ=255)
                                                                                               cGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAAc                                               >  2:575193/1‑149 (MQ=255)
                                                                                                  aaCTATTGATGAGAGTTTGGTGTCTTTAGGGGGCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTCCCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACttt                                            <  2:67842/149‑1 (MQ=255)
                                                                                                  aaCTATTGAGGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACttt                                            <  2:39891/149‑1 (MQ=255)
                                                                                                     tATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAAc                                               >  1:453328/1‑143 (MQ=255)
                                                                                                                gTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac                                    <  1:573993/143‑1 (MQ=255)
                                                                                                                gTTTGGTGTCTTTATGTGCCTGCAGGCATCTTACCATTCAAACTTTCGTGTCTCTACCGTGAAAGTGATACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac                                    <  1:428148/143‑1 (MQ=255)
                                                                                                                        tCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAt                      <  2:459652/149‑1 (MQ=255)
                                                                                                                            tATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTcc                        <  1:423527/143‑1 (MQ=255)
                                                                                                                               gtgtCTGCAGGCATCTTTTCATTCAAACTAAAGTTTTGCTCCCGTGAAAGTGCTACAAAGATAAAGGTTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAAACTTGTGAACACCGTCATATTTcacta                     >  1:412620/1‑139 (MQ=255)
                                                                                                                                         gCATCTTTCCATTCAAACGCACGTTTCGCGACCGTGAAAGTGCGACAAAGATAAAGATTTGGTTTAAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATAtt     <  2:515772/149‑1 (MQ=255)
                                                                                                                                            tCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTAc  >  2:170340/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTAC  >  CP009273/982416‑982704

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: