Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,328,417 G→A R188C (CGC→TGC)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,328,4170GA98.1% 191.2 / ‑3.9 54R188C (CGC→TGC) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (18/34);  total (18/36)
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.85e-01

TCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATC  >  REL606/3328344‑3328486
                                                                         |                                                                     
tCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGcaca                                                                     >  2:142266/1‑76 (MQ=255)
 cACCGATCAGAGAACCGTAATTACGACGCACATAGTTGCTGATAGCTTCGTCGAAATGGTCACCACCAATGCTcag                                                                    >  2:1986696/1‑75 (MQ=255)
   ccGATCAGAGAACCGTAATTACGACGCAGATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACag                                                                  <  1:306423/76‑1 (MQ=255)
   ccGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACag                                                                  <  1:909823/76‑1 (MQ=255)
    cGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACaga                                                                 <  1:680273/76‑1 (MQ=255)
         agagAACTGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAgagg                                                            <  2:821498/76‑1 (MQ=255)
         agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAgagg                                                            <  1:1071013/76‑1 (MQ=255)
         agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAgagg                                                            >  2:423653/1‑76 (MQ=255)
            gaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGt                                                         <  2:1516035/76‑1 (MQ=255)
             aaCCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTa                                                        >  2:1379402/1‑76 (MQ=255)
               ccGTAATTACGNCGNACATAGTTGATCATAGCTTCGTCGAAACGGTCACCACCAATGCTCACAGAAGAGGAGTaaa                                                      <  2:306970/76‑1 (MQ=255)
               ccGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTaaa                                                      <  1:1378567/76‑1 (MQ=255)
               ccGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTaaa                                                      >  2:1567492/1‑76 (MQ=255)
                 gTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAAcc                                                    >  1:1379880/1‑76 (MQ=255)
                 gTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAAcc                                                    <  1:634158/76‑1 (MQ=255)
                  tAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACca                                                   <  2:1400054/76‑1 (MQ=255)
                   aaTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACcac                                                  <  2:1854378/76‑1 (MQ=255)
                   aaTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACcac                                                  <  1:489151/76‑1 (MQ=255)
                    aTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACcaca                                                 <  1:945512/76‑1 (MQ=255)
                      tACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACAcc                                               >  1:1336626/1‑76 (MQ=255)
                          acgCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTc                                           >  1:376541/1‑76 (MQ=255)
                            gCACATAGTTGATGNTAGCTTCGTCGAAACGGTCACCACCAATGCACTCAGAAGAGGAGTAAACCAAACCNTTCaa                                         <  2:2043244/76‑1 (MQ=255)
                              acaTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAgg                                       >  1:765540/1‑76 (MQ=255)
                                aTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGgag                                     <  1:773418/76‑1 (MQ=255)
                                 tAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGgaga                                    >  1:576525/1‑76 (MQ=255)
                                  aGTTGATGCTAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGAt                                   >  1:865994/1‑76 (MQ=255)
                                  aGTTGATGATAGNTTNGTCGAAACGGTCACCACCNATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGAt                                   <  1:1668368/76‑1 (MQ=255)
                                  aGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACNCAGAAGAGGAGTAAACCANACNGTTCAAGGAGAt                                   >  1:1658290/1‑76 (MQ=255)
                                  aGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGAt                                   <  1:904758/76‑1 (MQ=255)
                                   gTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGANATa                                  <  1:1823185/76‑1 (MQ=255)
                                   gTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATa                                  >  1:875508/1‑76 (MQ=255)
                                   gTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATa                                  >  1:562527/1‑76 (MQ=255)
                                    ttGATGATAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATaa                                 <  1:1263478/76‑1 (MQ=255)
                                        tgatAGCTTCGGCGAAACGGTCACCACCAATGCACACTGAAGAGGAGTAAACCACACCGTTCAAGTAGATAATAGc                             >  1:53298/1‑76 (MQ=255)
                                         gatAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCa                            <  1:1567493/76‑1 (MQ=255)
                                           tAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAAc                          >  2:1297560/1‑76 (MQ=255)
                                           tAGCTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAAc                          >  2:1586534/1‑76 (MQ=255)
                                              cTTCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTc                       <  2:1555972/76‑1 (MQ=255)
                                               ttCGTCGAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCa                      <  1:280848/76‑1 (MQ=255)
                                                 cgtcgAAACGGTCANCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGACAACAGCAACTTCAGt                    <  1:1748405/76‑1 (MQ=255)
                                                  gtcgAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTg                   <  1:1263724/76‑1 (MQ=255)
                                                  gtcgAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTg                   >  1:1014468/1‑76 (MQ=255)
                                                    cgAAACGGTCACCACCAATGCACACAGNAGAGGAGTAAAACACACCGTTCAAGGAGATAACAGCAACTTCAGTGGt                 >  1:1784116/1‑76 (MQ=255)
                                                     gAAACGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTa                <  2:1265877/76‑1 (MQ=255)
                                                      aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac               <  2:425661/76‑1 (MQ=255)
                                                      aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac               >  2:551135/1‑76 (MQ=255)
                                                      aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac               <  1:1140040/76‑1 (MQ=255)
                                                      aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac               <  1:80503/76‑1 (MQ=255)
                                                      aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac               <  1:442121/76‑1 (MQ=255)
                                                      aaaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac               <  2:1337320/76‑1 (MQ=255)
                                                       aaCGGTCACCACCAATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTacc              <  1:973218/76‑1 (MQ=255)
                                                        aCGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTacca             <  1:686626/76‑1 (MQ=255)
                                                               ccaccaATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGa      <  2:652510/76‑1 (MQ=255)
                                                                 accaATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGata    <  2:1482079/76‑1 (MQ=255)
                                                                  ccaATGCACACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGatat   <  2:905720/76‑1 (MQ=255)
                                                                   caATGCACACAGACGAGGAGNAAACCACACCGTTCANGGAGATAACAGCAGCGTCAGTCGTACCACCACCGATATc  <  2:2018476/76‑1 (MQ=37)
                                                                         |                                                                     
TCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATC  >  REL606/3328344‑3328486

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