Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 44,301 (T)9→8 intergenic (+132/‑359) oprD_1 → / → czcR_2 Porin D/Transcriptional activator protein CzcR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB12244,2930T.94.6% 224.1 / 9.6 59intergenic (+124/‑367)oprD_1/czcR_2Porin D/Transcriptional activator protein CzcR
Reads supporting (aligned to +/- strand):  ref base T (3/0);  major base . (22/31);  minor base G (3/0);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 8.30e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGAAGCC  >  ACB122/44173‑44421
                                                                                                                        |                                                                                                                                
gCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCGTTTTTTTGCCTaaa                                                                                                                   >  1:66885/1‑135 (MQ=255)
  tAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCGTTTTTTTGCCTtaaa                                                                                                                  >  5:247212/1‑134 (MQ=255)
     cTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAAtt                                                                                                                >  7:48953/1‑133 (MQ=255)
     cTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCa                                                                                                              <  6:191690/135‑1 (MQ=255)
     cTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCa                                                                                                              <  3:153375/135‑1 (MQ=255)
          ccAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCGCCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAg                                                                                                         <  2:299327/135‑1 (MQ=255)
          ccAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAg                                                                                                         <  8:48953/135‑1 (MQ=255)
          ccAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAg                                                                                                         <  7:244326/135‑1 (MQ=255)
           cAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAgg                                                                                                        >  4:227998/1‑135 (MQ=255)
            aGTTATGCATTGGCGGGTGCTGATAGTAAAACGAGCATTACGCATTCAGCTTGAGTGGCTTAAGTGCTGCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                                       <  4:47166/135‑1 (MQ=255)
            aGTTATCGATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                                       <  2:274401/135‑1 (MQ=255)
            aGTTATCCCATGGCCGGTCGTGTTAGAAAAACGAGCATTCCGGATACCCCTTGAGTGGGGTATGTGATCCTTTGGTAAAAGGAGGATATTAGGCGCCCATCGCGCGCCGTTTTTTTGCCTAAAATTCATGTAGGc                                                                                                       <  4:174557/135‑1 (MQ=255)
            aGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                                       <  1:201416/135‑1 (MQ=255)
                      tGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCTAAAATTAAAGTAGGCGTTCGATTGt                                                                                             >  3:338689/1‑135 (MQ=255)
                         cGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCGTTTTTTTTCCTAAAATTAAAGTAGGCGTTCGATTGTAcg                                                                                          >  2:127725/1‑135 (MQ=255)
                         cGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGttc                                                                                           >  8:33061/1‑132 (MQ=255)
                           ggTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTCCTAAAATTAAAGTTGGAGTTGGATTTTAagca                                                                                        >  2:368255/1‑135 (MQ=255)
                           ggTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgca                                                                                        >  8:136231/1‑135 (MQ=255)
                           ggTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCGTTTTTTTGCCTAAAATTCAAGTAGGGGTTCGATTGTAcgca                                                                                        >  5:203554/1‑135 (MQ=255)
                             tCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatc                                                                                      <  7:147921/135‑1 (MQ=255)
                                tGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTAcgcatcgca                                                                                   >  4:319618/1‑135 (MQ=255)
                                      tAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                              <  7:147009/134‑1 (MQ=255)
                                      tAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                              <  1:127725/134‑1 (MQ=255)
                                          aCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                                       >  2:318466/1‑105 (MQ=255)
                                          aCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGc                                                                                                       <  1:318466/105‑1 (MQ=255)
                                          aCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                         <  5:10086/135‑1 (MQ=255)
                                          aCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                         >  6:386874/1‑135 (MQ=255)
                                          aCGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGaa                                                                         <  3:319618/135‑1 (MQ=255)
                                            gAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGt                                                                       >  2:312928/1‑135 (MQ=255)
                                             aGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCGGTTTTTTTGGCTTAAATTGAAGTAGGCGTTGCATTGTACGGCTCGCCTGCCGGTGAAGt                                                                       >  6:303790/1‑134 (MQ=255)
                                                 ttACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATg                                                                  <  7:140190/135‑1 (MQ=255)
                                                  tACGCATTCACCTTGAGTTGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGa                                                                 <  1:56726/135‑1 (MQ=255)
                                                       aTTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTaca                                                            <  7:136231/135‑1 (MQ=255)
                                                            ccTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCaa                                                        >  7:101455/1‑134 (MQ=255)
                                                            ccTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAAt                                                       <  7:272305/135‑1 (MQ=255)
                                                               tGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTAAACCAATGAc                                                    >  4:266137/1‑135 (MQ=255)
                                                                 aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                    <  8:313405/133‑1 (MQ=255)
                                                                 aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                    <  7:68529/133‑1 (MQ=255)
                                                                 aGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                    <  5:163863/133‑1 (MQ=255)
                                                                  gTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAg                                                 >  8:127962/1‑135 (MQ=255)
                                                                    ggCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTGGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGtc                                               >  7:172931/1‑135 (MQ=255)
                                                                    ggCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGt                                                <  8:101455/134‑1 (MQ=255)
                                                                         aaGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGCGAAGTTCATGACTACACCAATGAc                                                    >  4:36903/1‑125 (MQ=255)
                                                                         aaGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                    <  5:389039/125‑1 (MQ=255)
                                                                         aaGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGAc                                                    >  6:389039/1‑125 (MQ=255)
                                                                         aaGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGCGAAGTTCATGACTACACCAATGAc                                                    <  3:36903/125‑1 (MQ=255)
                                                                              cTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTAAAGTAGTCGTTCGATTTTACTCATCGCATGCCGGTGAATTTCATGACTTCACCAAATACAACTCTCGTCaaat                                      >  6:905/1‑132 (MQ=255)
                                                                                  tttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGGCTAAAATTAAAGTAGGCGTTCGATTGTACGCATCGCAAGCCGGTGAAGTTCATGACTACACCAATTAACAGAATCTTCAATGTGAGc                                 >  6:161353/1‑135 (MQ=255)
                                                                                  tttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGtt                                                                                                    >  1:280762/1‑68 (MQ=255)
                                                                                  tttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGtt                                                                                                    <  2:280762/68‑1 (MQ=255)
                                                                                   ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGGCTAAAATTGAAGTAGGCGTTAGATTGTACGCATGGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTATCTTCAATTTGAGCg                                >  1:245226/1‑135 (MQ=255)
                                                                                   ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTGGATTGTACGCATCGCATGCCGGTGAAGTTCATGAATACACAAATGACCAGTCTCGTCAATGTGAGCg                                >  6:56395/1‑135 (MQ=255)
                                                                                   ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                              <  4:370045/89‑1 (MQ=255)
                                                                                   ttGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCg                                                                              >  3:370045/1‑89 (MQ=255)
                                                                                        aaaaGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTc                           <  5:167507/135‑1 (MQ=255)
                                                                                                         cgcCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCgg           <  1:312928/134‑1 (MQ=255)
                                                                                                         cgcCCATTGGGCGCCTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCgg           <  2:66885/134‑1 (MQ=255)
                                                                                                             cATTGGGCGGCGTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGa      <  1:362108/135‑1 (MQ=255)
                                                                                                                 gggCGCCGTTTTTTTTGCGAAAATTCAAGTAGGCGTTGGATTGTACGCATCGCATGCCGGTGTAGTTGATGACTAAACCAATGACCAGTCTCGTCAATGTTAGCGCCCTCGCCGTCATCTGAGCGGTACGAAGcc  >  2:221893/1‑135 (MQ=255)
                                                                                                                        |                                                                                                                                
GCTAACTGGGCCAGTTATCCATTGGCGGGTCCTGTTAGTAAAACGAGCATTACGCATTCACCTTGAGTGGCTTAAGTGCTCCTTTGGTAAAAGGTGGATATTTGGCGCCCATTGGGCGCCTTTTTTTTTGCCTAAAATTCAAGTAGGCGTTCGATTGTACGCATCGCATGCCGGTGAAGTTCATGACTACACCAATGACCAGTCTCTTCAATGTGAGCGACCTCGCCGTCATCTGAGCGGAACGAAGCC  >  ACB122/44173‑44421

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