Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 577,885 T→A 48.0% I280I (ATT→ATA macB → fused macrolide transporter subunits and ATP‑binding component and membrane component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE577,8850TA48.0% 0.0 / 57.4 50I280I (ATT→ATAmacBfused macrolide transporter subunits and ATP‑binding component and membrane component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (9/17);  new base A (7/17);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACA  >  minE/577819‑577938
                                                                  |                                                     
gACGATGGGATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATtattat                                                     <  1:129933/69‑1 (MQ=255)
gACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATtattat                                                     <  1:1751625/69‑1 (MQ=255)
         atggcGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGc                                            <  1:1323664/69‑1 (MQ=255)
               ggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggt                                      <  1:896416/69‑1 (MQ=255)
               ggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggt                                      <  1:2932724/69‑1 (MQ=255)
               ggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggt                                      <  1:2094643/69‑1 (MQ=255)
               ggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggt                                      <  1:2957719/69‑1 (MQ=255)
               ggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggt                                      <  1:3098138/69‑1 (MQ=255)
               ggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggt                                      <  1:3177726/69‑1 (MQ=255)
               ggcgCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggt                                      <  1:1601675/69‑1 (MQ=255)
                 cgcTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCggtgg                                    >  1:1442721/1‑69 (MQ=255)
                     ggCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTc                                <  1:642795/69‑1 (MQ=255)
                     ggCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTc                                <  1:808503/69‑1 (MQ=255)
                     ggCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTc                                <  1:261374/69‑1 (MQ=255)
                     ggCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTc                                <  1:2525838/69‑1 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTc                         >  1:2444572/1‑69 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTACATTGTc                         >  1:3337886/1‑69 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATAATCgg                                                  >  1:3004449/1‑44 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATAATCgg                                                  >  1:2110700/1‑44 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATAATCgg                                                  >  1:2192019/1‑44 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATAATCgg                                                  >  1:265051/1‑44 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATAATCgg                                                  >  1:2239839/1‑44 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACCATGCTGGGGATTATAATCgg                                                  >  1:2508178/1‑44 (MQ=255)
                            aaTAAAATGCGTACTTTACTGACAATGCTGGGGATTATAATCgg                                                  >  1:640216/1‑44 (MQ=38)
                             aTAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCg                        >  1:360717/1‑69 (MQ=255)
                             aTAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCg                        >  1:954744/1‑69 (MQ=255)
                             aTAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCg                        >  1:2414382/1‑69 (MQ=255)
                             aTAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCg                        >  1:222601/1‑69 (MQ=255)
                             aTAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCg                        >  1:1792908/1‑69 (MQ=255)
                               aaaaTGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCgtg                      <  1:3182731/69‑1 (MQ=255)
                               aaaaTGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCgtg                      <  1:1499048/69‑1 (MQ=255)
                               aaaaTGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCgtg                      <  1:2020100/69‑1 (MQ=255)
                                    gCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGtggtgg                  >  1:2989175/1‑68 (MQ=255)
                                                   cATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCcaaaca  <  1:210/69‑1 (MQ=255)
                                                    aTGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCcaaaca  <  1:1685828/68‑1 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:1412557/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:1053907/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:778569/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:733398/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:708562/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:1329413/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:1395814/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:468388/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:2620143/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:3225983/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:3156196/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:2683996/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:3003865/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGTCCGTGGTTTCCATTGTCGTggtgggc                <  1:2476200/52‑2 (MQ=255)
                                                      gCTGGGGATTATAATCGGTATTGCGGCGGTGGTTTCCATTGTCGTggtgggc                <  1:1005433/52‑2 (MQ=255)
                                                       cTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:1980843/51‑2 (MQ=255)
                                                       cTGGGGATTATAATCGGTATTGCGTCGGTGGTTTCCATTGTCGTggtgggc                <  1:1270155/51‑2 (MQ=255)
                                                                  |                                                     
GACGATGGCATGGCGGGCGCTGGCAGCGAATAAAATGCGTACTTTACTGACCATGCTGGGGATTATTATCGGTATTGCGTCGGTGGTTTCCATTGTCGTGGTGGGTGACGCCGCCAAACA  >  minE/577819‑577938

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: