Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 982,560 T→A 24.0% 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,5600TA24.0% 86.9 / 25.9 50intergenic (‑122/+481)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (20/18);  new base A (9/3);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACA  >  CP009273/982415‑982705
                                                                                                                                                 |                                                                                                                                                 
tGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCa                                                                                                                                                <  2:1162973/149‑1 (MQ=255)
cgCCAAAATATTGCCCTTCACCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCa                                                                                                                                                <  2:960666/148‑1 (MQ=255)
 gCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAt                                                                                                                                               >  1:365964/1‑149 (MQ=255)
  ccAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAtt                                                                                                                                              >  2:53952/1‑149 (MQ=255)
  ccAGAATATTCCGCTTCATCATTATTTATTGCCCTCTTGGTTTTTTTTATGACACCTGCCGCTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCAtt                                                                                                                                              <  1:801178/149‑1 (MQ=255)
   cAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGGTTGGTGTCTTAATGTGTCTGCAGGCATCTTTCAATTc                                                                                                                                             >  2:826400/1‑149 (MQ=255)
         atTTCGCTTAATCCTTCTTTATTACCCTCAGGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTa                                                                                                                                       <  2:500794/149‑1 (MQ=255)
              gcTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTATCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGtt                                                                                                                                  >  1:749444/1‑149 (MQ=255)
                 tcgtcatTTTGTATTACCACCATGTTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCg                                                                                                                               <  2:43854/146‑1 (MQ=255)
                                 ccatcttgttttttctcATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGACTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCt                                                                                                               <  2:832141/140‑1 (MQ=255)
                                         gTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGCGTGTTTATGTGGCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACaaagat                                                                                                       <  1:1095552/149‑1 (MQ=255)
                                           tttttttATGACACCTGCCAATGCCGTCAATAACACGAATCCCTAAAAATTTACGTAACTTTTAATGAGAGATTGGTGTCTTTCTGTGTCTGCAGGCATCTTTACATTAAAACTAACGTTTCGATACCGTGAAAGTGCTACcaagataa                                                                                                     >  2:1134190/1‑149 (MQ=255)
                                            ttttttATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTCTCGCTACCGTGAAAGTGCTACaaagataaa                                                                                                    <  2:98605/149‑1 (MQ=255)
                                                  aTGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTg                                                                                              <  1:825696/149‑1 (MQ=255)
                                                  aTGACACCTGCCAATGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGTGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTTCCGTGAAAGTGCTACAACAATAAAAATTTg                                                                                              >  2:759960/1‑149 (MQ=255)
                                                    gaactccggcACGGCTGTCAATGAGGTCTGGCAATAAAAATTTACGGAACCATTGATGAGAGTTTGGTGTTTTTATGTGTCTTCAGGCATTTTTCCATTCAAACCAAGGTTTCGCTACCGTGAAAGTGCTATAAAGATAAAGATTTGGt                                                                                            <  2:325296/140‑1 (MQ=255)
                                                        ccTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCa                                                                                        >  2:937033/1‑149 (MQ=255)
                                                                   cGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGaaaaaa                                                                             >  1:713887/1‑149 (MQ=255)
                                                                    gTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTATGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTAGCTACGGTGCAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAat                                                                            >  2:133954/1‑149 (MQ=255)
                                                                     tCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTTGTTTCAAAAAGAAAAAata                                                                           >  2:227498/1‑149 (MQ=255)
                                                                      cAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAatat                                                                          >  1:1232238/1‑149 (MQ=255)
                                                                          aaGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTACTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTaa                                                                      >  2:1232067/1‑149 (MQ=255)
                                                                               cTGGCAATAAAAATTTACGTATCTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCCAACTAACGTTGTGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCaaa                                                                 <  1:826400/149‑1 (MQ=255)
                                                                                tGTCAATAAAAATTTACGGAACTATTGATGAGTGTTTGGTGTCGTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTTAAAGTTCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCaaaa                                                                <  2:184021/149‑1 (MQ=255)
                                                                                       aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAt                                                         >  1:928122/1‑149 (MQ=255)
                                                                                       aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAAt                                                         >  1:965971/1‑149 (MQ=255)
                                                                                       aaaaaTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATAGGTAACCAAAAGTAAAAt                                                         >  1:934135/1‑149 (MQ=255)
                                                                                           aTTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTaa                                                     <  1:53952/149‑1 (MQ=255)
                                                                                            tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAg                                                    >  1:1147823/1‑149 (MQ=255)
                                                                                            tttACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATGTAAg                                                    >  2:246418/1‑149 (MQ=255)
                                                                                               acagAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGaa                                                 <  1:618286/146‑1 (MQ=255)
                                                                                                      tATTGATAAGAGTTTGGTGTCTTTATGTTTCAGGAGGCATCCTACCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGAGAAAAAAAAAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTtgtg                                          <  2:809032/149‑1 (MQ=255)
                                                                                                           atgaGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAacac                                     >  2:647634/1‑149 (MQ=255)
                                                                                                              agagTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGt                                  >  1:941751/1‑149 (MQ=255)
                                                                                                                agTTTGGTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCa                                >  2:1101169/1‑149 (MQ=255)
                                                                                                                  tttGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCCAACAGAAAAAAGATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCata                              >  1:1158681/1‑149 (MQ=255)
                                                                                                                  tttGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCata                              >  1:248575/1‑149 (MQ=255)
                                                                                                                     ggTGTCTTTATGTGTCTGCAGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATAttt                           >  1:1170635/1‑149 (MQ=255)
                                                                                                                      gtgcccttATGTGTCTGGAGGGGACGGTCCACTCAATCTACCGTTTCCCGACCGTGAAAGTGCTAAAGAGATAAAGATTTGGAAAAAAAAAAAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTc                          <  2:335631/143‑1 (MQ=255)
                                                                                                                        gtCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACAGTCATATTTCCa                        >  1:966114/1‑149 (MQ=255)
                                                                                                                        gtCTTTATGTGTCAGCAGGCATCTTTCCATTAAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTAAAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCa                        <  2:1227370/149‑1 (MQ=255)
                                                                                                                         tCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGCTTTGGTTTCAAAAAGAAAAAAGATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCAt                       >  1:664963/1‑149 (MQ=255)
                                                                                                                           tttATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATag                     >  1:267640/1‑149 (MQ=255)
                                                                                                                                   tCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAGGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGtga             >  1:1024124/1‑149 (MQ=255)
                                                                                                                                   attGCAGGCAGCTTTCCAGTCAAAGTAACGTGTCGGTGCCGTGAAAGTGATAAAAAGATAAAGATTTAAATAAAAAAAAAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGtga             <  2:444280/147‑1 (MQ=255)
                                                                                                                                      gCAGGCCTCTTTCCATTCAAACTAACGTTTCGCCACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGCAAAAATATTTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGtgatga          <  1:292509/149‑1 (MQ=255)
                                                                                                                                       caagCATCTTTCCATTCAAACTAACGTTTCACAACCGTGATAGTGCTACAAAGATAAAGATTCGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTgatgat         <  1:63744/146‑1 (MQ=255)
                                                                                                                                        aGGCATCTTACCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGAACTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGata        >  2:450010/1‑149 (MQ=255)
                                                                                                                                             tCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTAc   >  1:1018091/1‑149 (MQ=255)
                                                                                                                                              cTTTCCATTCAAACTAACGTTTCGCTACCGTGAAATTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATACAGACGTGATGATATTTACa  <  1:898689/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGCCAGAATATTGCGCTTCATCATTATTTATTACCCTCATGGTTTTTTTTATGACACCTGCCACTGCCGTCAATAAGTTCTGTCAATAAAAATTTACGGAACTATTGATGAGAGTTTGGTGTCTTTATGTGTCTGCAGGCATCTTTCCATTCAAACTAACGTTTCGCTACCGTGAAAGTGCTACAAAGATAAAGATTTGGTTTCAAAAAGAAAAAATATGTAACCAAAAGTAAAATTTAAGGAACTTTGTGAACACCGTCATATTTCCATAGAGACGTGATGATATTTACA  >  CP009273/982415‑982705

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: ‑