Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,532,298 A→T T15T (ACA→ACT)  murN → FemAB family peptidoglycan biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,532,2980AT97.8% 142.0 / ‑3.6 46T15T (ACA→ACT) murNFemAB family peptidoglycan biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (19/26);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

ATGTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532160‑2532433
                                                                                                                                           |                                                                                                                                       
aTGTGTATTTGTAATCAGTTTTTGATAAATTAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTtactac                                                                                                                                       >  2:53115/1‑141 (MQ=255)
 tgtgTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTactaca                                                                                                                                      <  2:153968/141‑1 (MQ=255)
    gtATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAg                                                                                                                                   >  2:442177/1‑141 (MQ=255)
     tATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGc                                                                                                                                  <  1:129488/141‑1 (MQ=255)
     tATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATAAATAAATT‑AATAATCATTTGATAAATTCAATGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGc                                                                                                                                  <  1:266952/141‑1 (MQ=255)
           tAAGCAGTTTGTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCTGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCAttttt                                                                                                                            <  2:168684/141‑1 (MQ=255)
            aaTCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTc                                                                                                                           >  1:440788/1‑141 (MQ=255)
                aGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACAt                                                                                                                       <  1:468966/141‑1 (MQ=255)
                    tttGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTaca                                                                                                                   >  2:101952/1‑141 (MQ=255)
                       gATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTacacac                                                                                                                <  1:510024/141‑1 (MQ=255)
                              taaaaAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGc                                                                                                         <  1:101952/141‑1 (MQ=255)
                                aaaaaaTAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTa                                                                                                       >  1:371364/1‑141 (MQ=255)
                                 aaaaaTAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTAt                                                                                                      >  1:291839/1‑141 (MQ=255)
                                  aaaaTAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTAtt                                                                                                     >  2:74413/1‑141 (MQ=255)
                                      tAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCAtt                                                                                                 <  1:424563/141‑1 (MQ=255)
                                        gCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                               <  2:371364/141‑1 (MQ=255)
                                        gCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTAc                                                                                               <  1:224536/141‑1 (MQ=255)
                                                  taaATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATagag                                                                                     <  2:440788/141‑1 (MQ=255)
                                                   aaaTTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGt                                                                                    <  2:103495/141‑1 (MQ=255)
                                                    aaTTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGtt                                                                                   <  1:194372/141‑1 (MQ=255)
                                                       tCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAt                                                                                <  1:180079/141‑1 (MQ=255)
                                                        cAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTc                                                                                                    >  1:498595/1‑120 (MQ=255)
                                                        cAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTc                                                                                                    <  2:498595/120‑1 (MQ=255)
                                                         aaaaaTTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGAttt                                                                              <  1:538721/141‑1 (MQ=255)
                                                          aaaaTTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTa                                                                             <  2:383667/141‑1 (MQ=255)
                                                            aaTTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGAGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTAGATTTaaa                                                                           >  2:395491/1‑141 (MQ=255)
                                                               ta‑taGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAgg                                                                        <  2:426640/141‑1 (MQ=255)
                                                                a‑taGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGc                                                                       <  1:101917/141‑1 (MQ=255)
                                                                  taGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCg                                                                      >  1:95032/1‑141 (MQ=255)
                                                                  taGTTAAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCg                                                                      <  2:210876/141‑1 (MQ=255)
                                                                    gTTGAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATCACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGAt                                                                    >  2:514784/1‑141 (MQ=255)
                                                                                            aaaaTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGAt                                            <  1:301005/141‑1 (MQ=255)
                                                                                               aTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                         <  1:91191/141‑1 (MQ=255)
                                                                                               aTGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGAc                                         <  2:374602/141‑1 (MQ=255)
                                                                                                     tttGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGt                                   <  2:509671/141‑1 (MQ=255)
                                                                                                                  cTTTGGATGAGTTTAATGTGGTTACTACATAACATTTTTTACATTAGACATAAATAGCTATTCATTATAATCATAGAGTTGATTTAAAATGCGATGTGCATCTTGTAGGGGATAAAGATGACAAGGGGGAAGGGATTGCAg                      >  2:333502/1‑141 (MQ=255)
                                                                                                                  cTTCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAg                      >  2:402653/1‑141 (MQ=255)
                                                                                                                   ttCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAgg                     >  2:220103/1‑141 (MQ=255)
                                                                                                                   ttCGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAgg                     <  1:152873/141‑1 (MQ=255)
                                                                                                                     cGGATGAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGAt                   <  1:39084/141‑1 (MQ=255)
                                                                                                                          gAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                            >  1:53504/1‑95 (MQ=255)
                                                                                                                          gAGTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCt                                                            <  2:53504/95‑1 (MQ=255)
                                                                                                                            gTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTAtt            >  1:209834/1‑141 (MQ=255)
                                                                                                                            gTTCAATGCGTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTAtt            >  2:496490/1‑141 (MQ=255)
                                                                                                                                aaTGCGTTTACTACAAAGCATTTTTCACATTAcacaca                                                                                                               >  1:393574/1‑38 (MQ=255)
                                                                                                                                aaTGCGTTTACTACAAAGCATTTTTCACATTAcacaca                                                                                                               <  2:393574/38‑1 (MQ=255)
                                                                                                                                     gTTTACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgc   >  1:182821/1‑141 (MQ=255)
                                                                                                                                      tttACTACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAgca  >  1:216694/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                       
ATGTGTATTTGTAATCAGTTTTTGATAAAATAAAAAAATAGCATTATTAATAAATTCAAAAATTA‑TAGTTCAATTCAGAGAGGGATGGGATAAAATGAACTTTGTAACGTTGACTTCGGATGAGTTCAATGCGTTTACAACAAAGCATTTTTCACATTACACACAATCAGCTATTCATTACAATCATAGAGTTGATTTAAAAGGCGATGTGCATCTTGTAGGGGTTAAAGATGACAATGGTCAAGTGATTGCAGGATGCTTATTGACAGAAGCA  >  CP000730/2532160‑2532433

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