Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,412,459 +A pseudogene (59/66 nt) yhdW → pseudogene, amino‑acid transporter homology;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,412,4561.A93.8% 144.3 / 4.2 49pseudogene (56/66 nt)yhdWpseudogene, amino‑acid transporter homology;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base A (27/18);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAGGTATTGAT  >  CP009273/3412311‑3412582
                                                                                                                                                  |                                                                                                                               
tactTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTGGc                                                                                                                                  >  1:399077/1‑146 (MQ=255)
   ttttGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATAATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTCGCCGTTGCacaagc                                                                                                                           >  2:303994/1‑144 (MQ=255)
    tttGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCaaaatca                                                                                                                          >  2:538353/1‑145 (MQ=255)
          tataAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATGGCCACGCTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCgct                                                                                                                    >  2:298013/1‑149 (MQ=255)
             aaGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCgctgct                                                                                                                 <  2:869840/149‑1 (MQ=255)
                     ggCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGCCCCGAc                                                                                                         >  2:848621/1‑149 (MQ=255)
                     ggCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAAGTTTCGCCTCCCAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGAc                                                                                                         >  2:729208/1‑149 (MQ=255)
                          cagacagGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGc                                                                                                       >  1:651386/1‑146 (MQ=255)
                               agGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCAAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGAt                                                                                                  >  1:51848/1‑146 (MQ=255)
                                 gTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGc                                                                                                >  1:969494/1‑146 (MQ=255)
                                    aaaGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAGTCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCa                                                                                          >  2:642004/1‑149 (MQ=255)
                                     aaGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTTGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGtt                                                                                            <  1:512511/146‑1 (MQ=255)
                                      aGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTc                                                                                           >  1:740717/1‑146 (MQ=255)
                                           tttCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGaaaaaa                                                                                   <  2:739624/149‑1 (MQ=255)
                                               ccAAACGTCCTATAACGATAATTTTCACCTCTCTGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATTCAAGCGCTGCTGGCGCCACGCTTGATGCCGGTCTGAAAATAg                                                                                  >  1:578647/1‑146 (MQ=255)
                                                      tCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAACTCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTtgt                                                                           <  1:501296/146‑1 (MQ=255)
                                                            aaCGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATg                                                                     >  1:1031901/1‑146 (MQ=255)
                                                             aCGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCggg                                                                 >  2:735589/1‑149 (MQ=255)
                                                             aCGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCggg                                                                 >  2:533035/1‑149 (MQ=255)
                                                               gATAATTTTAGCCTCACAGGAAGAATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGGGCTGCTGGCGCGACGGTTGATGCCGTTCAGAAAGAAGGTTTTGGGCAATGCgg                                                                  >  1:554610/1‑146 (MQ=255)
                                                                aTAATTTTCGCCTCACAGGAGGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATc                                                              >  2:222237/1‑149 (MQ=255)
                                                                aTAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCggg                                                                 >  1:860874/1‑146 (MQ=255)
                                                                  aaTTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGAt                                                               >  1:154221/1‑146 (MQ=255)
                                                                    ttttCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTg                                                          <  2:636597/149‑1 (MQ=255)
                                                                      ttCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGAt                                                        <  2:622965/149‑1 (MQ=255)
                                                                       tCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGACGTCATGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATg                                                       <  2:561474/149‑1 (MQ=255)
                                                                         gCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGa                                                     <  2:589327/149‑1 (MQ=255)
                                                                          ccgcacaGGAAGCATTATGAAAAAGATGATGGTGGCCACACTGGCTGCGGCCCACGTGCGGCGTGCCGTTGCAAAATCGGGCGCTGCTGGCGCGACGCTTGATGCCGGTCAGAAAAAAAGTTTTGTGCAATGCGGGATCAGTGATg                                                       <  1:604645/143‑1 (MQ=255)
                                                                          ccTCACGGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATg                                                       <  1:497819/146‑1 (MQ=255)
                                                                          ccTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATg                                                       <  1:1020238/146‑1 (MQ=255)
                                                                          ccTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATg                                                       >  1:151486/1‑146 (MQ=255)
                                                                               caGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTa                                                  >  1:669803/1‑146 (MQ=255)
                                                                                aGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTg                                              <  2:102708/149‑1 (MQ=255)
                                                                                    ggCATTATGAAAAAGGTGGTGATAGCCACCCTGGTGGCCGCAAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCGGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAGTGCGGGATCAGTGATGGATTACCTgg                                             <  1:131478/145‑1 (MQ=255)
                                                                                     gCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTggg                                            >  1:988069/1‑146 (MQ=255)
                                                                                        ttATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTctct                                      >  2:908956/1‑149 (MQ=255)
                                                                                            gAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCGAGATCGGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATg                                  >  2:509459/1‑149 (MQ=255)
                                                                                             aaaaaGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCGGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGAAGGATTACCTGGGTTCTCTTa                                    >  1:422890/1‑146 (MQ=255)
                                                                                             aaaaaGATGATGATAGCCACACTGGCGGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGc                                 <  2:711193/149‑1 (MQ=255)
                                                                                             aaaaaGAGGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGc                                 <  2:568667/149‑1 (MQ=255)
                                                                                                aaGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGc                                 <  1:441545/146‑1 (MQ=255)
                                                                                                     gatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGGGATGGATTACCTGGGTTCTCTTATGCCGATGCta                         >  2:292845/1‑148 (MQ=255)
                                                                                                     gatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGGGATGGAATACCTGGGTTCTCTTATGCCGATg                            >  1:317566/1‑146 (MQ=255)
                                                                                                             ccACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACggg                    >  1:333242/1‑145 (MQ=255)
                                                                                                             ccACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGt                    >  1:943166/1‑146 (MQ=255)
                                                                                                             ccACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAg                 >  2:579054/1‑149 (MQ=255)
                                                                                                             ccACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTCCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGt                    >  1:52702/1‑146 (MQ=255)
                                                                                                                  cTGGCTGCCGCCAGCTTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGtt               <  1:791637/146‑1 (MQ=255)
                                                                                                                   tGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCa           <  2:602235/149‑1 (MQ=255)
                                                                                                                     gCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCATTCAGAAAAAAGGTTTTGTGCAATGCGGGAGCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAgg         >  2:607876/1‑149 (MQ=255)
                                                                                                                          cgccAGTGTGCTGCTTGCCGTCGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAGGTATTg    <  2:696025/149‑1 (MQ=255)
                                                                                                                            ccAGCGTGCTGCTTGCCGTTGCAAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAGGTATTGAt  >  2:204747/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                               
TACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAGGTATTGAT  >  CP009273/3412311‑3412582

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