Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,029,236 T→C intergenic (+512/‑80) ompC → / → hchA porin OmpC/protein deglycase HchA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,029,2360TC95.7% 144.7 / ‑2.2 46intergenic (+512/‑80)ompC/hchAporin OmpC/protein deglycase HchA
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (21/23);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTTCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAATACACT  >  NZ_CP009273/2029150‑2029315
                                                                                      |                                                                               
gCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCa                                                                              >  2:141293/1‑90 (MQ=255)
gCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCa                                                                              >  1:185530/1‑90 (MQ=255)
    gAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCt                                                                          <  2:621626/90‑1 (MQ=255)
              cAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAAt                                                                >  2:753851/1‑90 (MQ=255)
              cAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAAt                                                                >  2:204092/1‑90 (MQ=255)
              cAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAAt                                                                >  1:27836/1‑90 (MQ=255)
              cAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAAt                                                                >  2:611552/1‑90 (MQ=255)
                     aGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGAcc                                                         >  2:627226/1‑90 (MQ=255)
                     aGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGAcc                                                         <  2:195629/90‑1 (MQ=255)
                     aGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCGGGTTAAAAATAGCGAcc                                                         >  2:110086/1‑90 (MQ=255)
                      gTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCa                                                        >  1:167037/1‑90 (MQ=255)
                              tcatcatTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTc                                                <  2:149907/90‑1 (MQ=255)
                                atcatTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGc                                              >  2:580969/1‑90 (MQ=255)
                                atcatTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGc                                              >  1:383973/1‑90 (MQ=255)
                                atcatTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGc                                              <  1:627226/90‑1 (MQ=255)
                                  catTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCaa                                            <  1:634848/90‑1 (MQ=255)
                                     tACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAata                                         >  2:105122/1‑90 (MQ=255)
                                       cATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAAtata                                       <  1:510554/90‑1 (MQ=255)
                                       cATACTTACAACCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAAtata                                       <  1:505112/90‑1 (MQ=255)
                                        aTACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAg                                      >  1:300905/1‑90 (MQ=255)
                                           cTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGa                                   <  2:300905/90‑1 (MQ=255)
                                           cTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGa                                   <  1:539473/90‑1 (MQ=255)
                                           cTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGa                                   <  1:355032/90‑1 (MQ=255)
                                           cTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGa                                   >  2:157668/1‑90 (MQ=255)
                                            ttACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGAc                                  <  2:257558/90‑1 (MQ=255)
                                            ttACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGAc                                  >  2:166615/1‑90 (MQ=255)
                                                aGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTAcc                              >  1:622622/1‑90 (MQ=255)
                                                   cgccgACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCcta                           >  1:765481/1‑90 (MQ=255)
                                                   cgccgACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCcta                           >  1:131520/1‑90 (MQ=255)
                                                    gccgACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCctaa                          >  2:528788/1‑90 (MQ=255)
                                                     ccgACTACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCctaac                         <  1:675507/90‑1 (MQ=255)
                                                          tACCGCTTTCGATATATTGACTCAATTTTCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGc                    <  1:379676/90‑1 (MQ=255)
                                                          tACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGc                    <  2:44790/90‑1 (MQ=255)
                                                          tACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGc                    <  1:250453/90‑1 (MQ=255)
                                                          tACCGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGc                    <  1:611552/90‑1 (MQ=255)
                                                             cGCTTTCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGcaac                 <  2:27836/90‑1 (MQ=255)
                                                                 ttCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATa             <  2:383973/90‑1 (MQ=255)
                                                                 ttCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATa             <  1:186651/90‑1 (MQ=255)
                                                                 ttCGATATATTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATa             <  1:166615/90‑1 (MQ=255)
                                                                      tatatTGACTCAATTTCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAAt       <  1:553718/90‑1 (MQ=255)
                                                                       atatTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAAt       <  2:131520/90‑1 (MQ=255)
                                                                       atatTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAAt       <  1:576282/90‑1 (MQ=255)
                                                                       atatTGACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAAt       <  1:204092/90‑1 (MQ=255)
                                                                            gACTCAATTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAATACACt  >  2:567825/1‑90 (MQ=255)
                                                                                 aaTTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAAt       <  1:473960/80‑1 (MQ=255)
                                                                                 aaTTTCCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAAt       >  2:473960/1‑80 (MQ=255)
                                                                                      |                                                                               
GCGTGAATAATATTCAACACCAGTGCAGCATCATCATTACATACTTACAGCCGCCGACTACCGCTTTCGATATATTGACTCAATTTTCCAGCCTGGTTAAAAATAGCGACCAGCTCAGTCGCAAATATAGTGACTACCCTAACTAAGCAACAATAAGGAATACACT  >  NZ_CP009273/2029150‑2029315

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