Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,626,268 A→G pseudogene (359/1368 nt) ydfJ ← MHS family MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,626,2680AG94.7% 121.3 / ‑0.6 38pseudogene (359/1368 nt)ydfJMHS family MFS transporter
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/25);  total (12/26)
Fisher's exact test for biased strand distribution p-value = 5.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

GCATATTCCGCCAGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGAT  >  NZ_CP009273/1626187‑1626347
                                                                                 |                                                                               
gcATATTCCGCCAGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGt                                                                         <  2:537163/90‑1 (MQ=255)
   tatTCCGCCAGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCa                                                                      <  2:148614/90‑1 (MQ=255)
     ttCCGCCAGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCaaa                                                                    <  2:397925/90‑1 (MQ=255)
     ttCCGCCAGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCaaa                                                                    <  1:684852/90‑1 (MQ=255)
           cAGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACAcc                                                              <  1:562136/90‑1 (MQ=255)
            aGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCAAGATGGGTGCAAAAACACCa                                                             <  2:6270/90‑1 (MQ=255)
               aTCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATc                                                          <  1:324029/90‑1 (MQ=255)
               aTCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATc                                                          <  1:757757/90‑1 (MQ=255)
                         cACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAgg                                                <  2:312856/90‑1 (MQ=255)
                           ccGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTc                                              <  1:743147/90‑1 (MQ=255)
                            cGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCg                                             >  2:768723/1‑90 (MQ=255)
                             ggAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGAGCATAGGTCgg                                            >  1:40357/1‑90 (MQ=255)
                                aaTTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTaa                                         <  2:417317/90‑1 (MQ=255)
                                 aTTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAac                                        <  2:126006/90‑1 (MQ=255)
                                 aTTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAac                                        <  1:889792/90‑1 (MQ=255)
                                    tCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAacaca                                     <  2:668575/90‑1 (MQ=255)
                                          ccTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAAtt                               <  2:907901/90‑1 (MQ=255)
                                            tGCACCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTaa                             >  1:629841/1‑90 (MQ=255)
                                               aCCCAACCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAggt                          <  1:369084/90‑1 (MQ=255)
                                                    aCCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTAc                     >  1:666107/1‑90 (MQ=255)
                                                    aCCCCTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTAc                     >  2:430641/1‑90 (MQ=255)
                                                    aCCCCTGAATAATACGCAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTAc                     >  1:355373/1‑90 (MQ=255)
                                                       ccTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAAtt                               <  1:288757/77‑1 (MQ=255)
                                                       ccTGAATAATACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAAtt                               >  2:288757/1‑77 (MQ=255)
                                                           aataatACGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATccc              <  2:372861/90‑1 (MQ=255)
                                                                taCGCAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCcatca         >  1:60171/1‑90 (MQ=255)
                                                                    cAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTAc                     <  2:48004/74‑1 (MQ=255)
                                                                    cAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTAc                     >  1:48004/1‑74 (MQ=255)
                                                                    cAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGt     >  1:38127/1‑90 (MQ=255)
                                                                    cAACGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGt     >  2:71878/1‑90 (MQ=255)
                                                                     aaCGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTg    <  2:675721/90‑1 (MQ=255)
                                                                     aaCGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTg    <  1:309632/90‑1 (MQ=255)
                                                                      aCGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAacaca                                     >  1:73455/1‑56 (MQ=255)
                                                                      aCGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAacaca                                     <  2:73455/56‑1 (MQ=255)
                                                                       cGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGAt  <  1:858851/90‑1 (MQ=255)
                                                                       cGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGAt  <  1:720977/90‑1 (MQ=255)
                                                                       cGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGAt  <  1:602583/90‑1 (MQ=255)
                                                                       cGTCACCAGCGAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGAt  <  1:135454/90‑1 (MQ=255)
                                                                                 |                                                                               
GCATATTCCGCCAGCATCGTACCGGCACCGGAAATTTCCGCACCTGCACCCAACCCCTGAATAATACGCAACGTCACCAGCAAGATGGGTGCAAAAACACCAATCTGTGCATAGGTCGGTAACACACCAATTAAGGTGGTACAGATCCCCATCATGGTGAT  >  NZ_CP009273/1626187‑1626347

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