Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,332,615 T→A 11.3% intergenic (‑25/+102) yhbE ← / ← rpmA conserved inner membrane protein/50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,332,6150TA11.3% 188.8 / 15.7 71intergenic (‑25/+102)yhbE/rpmAconserved inner membrane protein/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (4/4);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

CTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCG  >  AP009048/3332545‑3332684
                                                                      |                                                                     
cTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCAttt                                                                       <  1:3015660/71‑1 (MQ=255)
   tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTg                                                                      <  1:1820230/69‑1 (MQ=255)
   tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGc                                                                     <  1:1205504/70‑1 (MQ=255)
     gtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTAGCGGa                                                                  >  1:1501544/1‑71 (MQ=255)
       gAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGAt                                                                 >  1:1304204/1‑70 (MQ=255)
        aGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGAtt                                                                <  1:489202/70‑1 (MQ=255)
         gcgcCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTaa                                                              <  1:4057308/71‑1 (MQ=255)
          cgcCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAAc                                                             <  1:3861070/71‑1 (MQ=255)
             cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTc                                                          <  1:538794/71‑1 (MQ=255)
             cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTc                                                          >  1:2765363/1‑71 (MQ=255)
              aaaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCg                                                         >  1:4433435/1‑71 (MQ=255)
              aaaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCg                                                         >  1:1768590/1‑71 (MQ=255)
               aaaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCg                                                         >  1:4079586/1‑70 (MQ=255)
                aaGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCa                                                       >  1:1452629/1‑71 (MQ=255)
                  gAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCAcc                                                     <  1:1049844/71‑1 (MQ=255)
                  gAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCAc                                                      >  1:3166635/1‑70 (MQ=255)
                       ccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCAt                                                   >  1:1328817/1‑68 (MQ=255)
                       ccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATgcg                                                <  1:4108968/71‑1 (MQ=255)
                        caatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATgcg                                                >  1:4195488/1‑70 (MQ=255)
                          atgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTaa                                                              >  1:4196276/1‑54 (MQ=255)
                          atgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCa                                                    >  1:2436133/1‑64 (MQ=255)
                          atgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAt                                             <  1:2761564/71‑1 (MQ=255)
                          atgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAt                                             <  1:3495770/71‑1 (MQ=255)
                          atgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAt                                             >  1:282979/1‑71 (MQ=255)
                          atgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCAt                                             <  1:2835740/71‑1 (MQ=255)
                           tgcctgccTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtt                                            <  1:2680553/71‑1 (MQ=255)
                           tgcctgccTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtt                                            <  1:3056681/71‑1 (MQ=255)
                           tgcctgccTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtt                                            <  1:3075669/71‑1 (MQ=255)
                             cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtttt                                          <  1:3472586/71‑1 (MQ=255)
                             cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtttt                                          >  1:164298/1‑71 (MQ=255)
                             cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtttt                                          <  1:544296/71‑1 (MQ=255)
                              ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCAc                                                      <  1:2762809/58‑1 (MQ=255)
                              ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATgcgc                                               <  1:1243248/65‑1 (MQ=255)
                              ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATgcgc                                               <  1:4387989/65‑1 (MQ=255)
                              ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtttt                                          >  1:4035572/1‑70 (MQ=255)
                              ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAt                                             <  1:2712001/67‑1 (MQ=255)
                              ctgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCAttttt                                         <  1:2105963/71‑1 (MQ=255)
                               tgcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTg                                                                      >  1:1419660/1‑41 (MQ=255)
                                 cctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTAc                                       <  1:2134198/70‑1 (MQ=255)
                                  ctgctgCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAtttt                                          >  1:535564/1‑66 (MQ=255)
                                   tgctgcTTCATACCGTATTCTCTGCAAAGCGCATTTGCGGAt                                                                 <  1:3883794/42‑1 (MQ=255)
                                   tgctgcTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGAt                                                                 <  1:2501047/42‑1 (MQ=255)
                                    gctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAt                                             <  1:1959243/61‑1 (MQ=255)
                                    gctgctTCATACCGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCATTTTTACCGGt                                   >  1:2472270/1‑71 (MQ=255)
                                     ctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCgg                                                                   <  1:1492044/38‑1 (MQ=255)
                                     ctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGtt                                  <  1:4240906/71‑1 (MQ=255)
                                     ctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGtt                                  <  1:3520368/71‑1 (MQ=255)
                                     ctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGtt                                  <  1:3515199/71‑1 (MQ=255)
                                       gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTAt                                >  1:2016403/1‑71 (MQ=255)
                                        ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAt                                             <  1:3551537/57‑1 (MQ=255)
                                        ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTa                                        >  1:2183502/1‑62 (MQ=255)
                                          tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCAt                                                   >  1:1090082/1‑49 (MQ=255)
                                          tCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTAt                                >  1:1279057/1‑68 (MQ=255)
                                           cATACCGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGaa                            >  1:154926/1‑71 (MQ=255)
                                             tACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCAttttt                                         <  1:1380479/56‑1 (MQ=255)
                                              aCCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCa                                                    >  1:1221017/1‑44 (MQ=255)
                                               ccGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGt                         >  1:2738651/1‑70 (MQ=255)
                                               ccGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTa                        <  1:4175929/71‑1 (MQ=255)
                                                cGTATTCCCTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTaa                       <  1:1995016/71‑1 (MQ=255)
                                                cGTATTCACTGCAAAGCGCATTAGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTa                        >  1:1702118/1‑70 (MQ=255)
                                                  tATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTaaaa                     <  1:3192366/71‑1 (MQ=255)
                                                    ttCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTAcc                                      >  1:2932634/1‑52 (MQ=255)
                                                    ttCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAg                   <  1:46863/71‑1 (MQ=255)
                                                      cccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGcc                 >  1:1663341/1‑71 (MQ=255)
                                                        cTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTaaaa                     >  1:2340012/1‑65 (MQ=255)
                                                          gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCg              >  1:2019688/1‑70 (MQ=255)
                                                             aaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAAc          <  1:1726543/71‑1 (MQ=255)
                                                              aGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAAca         <  1:551219/71‑2 (MQ=255)
                                                               gcgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATc                               >  1:3632204/1‑48 (MQ=255)
                                                               gcgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATc                               >  1:3127551/1‑48 (MQ=255)
                                                               gcgcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACgt        <  1:1029988/71‑1 (MQ=255)
                                                                   aTTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTg    <  1:289853/71‑1 (MQ=255)
                                                                    tttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGcccc               <  1:1098720/59‑1 (MQ=255)
                                                                     ttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAg                   >  1:4553732/1‑54 (MQ=255)
                                                                      tGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCg  <  1:482126/70‑1 (MQ=255)
                                                                      |                                                                     
CTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCG  >  AP009048/3332545‑3332684

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