Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 120,479 G→A 10.5% A358V (GCT→GTT)  aroP ← aromatic amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb120,4790GA10.5% 153.8 / 7.6 57A358V (GCT→GTT) aroParomatic amino acid transporter
Reads supporting (aligned to +/- strand):  ref base G (29/22);  new base A (2/4);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

CGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGACACCAGAATGG  >  W3110S.gb/120412‑120545
                                                                   |                                                                  
cgcCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCg                                                                  <  1:775876/70‑1 (MQ=255)
cgcCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCg                                                                  <  1:1676348/70‑1 (MQ=255)
 gcCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGa                                                                <  1:3562820/71‑1 (MQ=255)
         aaTCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTgggg                                                        >  1:2075254/1‑71 (MQ=255)
             atcGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAg                                                    >  1:3103113/1‑71 (MQ=255)
             atcGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAg                                                    >  1:2942604/1‑71 (MQ=255)
             atcGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAg                                                    >  1:2322084/1‑71 (MQ=255)
             atcGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAg                                                    >  1:1774488/1‑71 (MQ=255)
                 cccAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAg                                                <  1:2304229/71‑1 (MQ=255)
                    aGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGt                                                  <  1:464805/66‑1 (MQ=255)
                    aGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTa                                             <  1:1288954/71‑1 (MQ=255)
                     gTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGc                                                       <  1:2139925/60‑1 (MQ=255)
                        gATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAg                                                    >  1:1877899/1‑60 (MQ=255)
                         aTTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAg                                        <  1:1374351/71‑1 (MQ=255)
                         aTTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAg                                        <  1:3672495/71‑1 (MQ=255)
                            aCCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTAc                                     >  1:3614162/1‑71 (MQ=255)
                                   cAGATACCACCAGCGCCATTAACAGTCCGAAAGCg                                                                  >  1:3332038/1‑35 (MQ=255)
                                   cAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTg                                                           >  1:3634795/1‑42 (MQ=255)
                                   cAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTg                                                           >  1:3639269/1‑42 (MQ=255)
                                   cAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTg                                                           >  1:2713894/1‑42 (MQ=255)
                                   cAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAAc                              >  1:1708973/1‑71 (MQ=255)
                                   cAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAAc                              >  1:1072160/1‑71 (MQ=255)
                                   cAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAAc                              >  1:1803277/1‑71 (MQ=255)
                                     gATACCACCAGCGCCATTAACAGTCCGTAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGc                                   <  1:1303808/64‑1 (MQ=255)
                                     gATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCTAGGTAGTTAATCAGTACGc                                   <  1:70697/64‑1 (MQ=255)
                                     gATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGc                                   <  1:1525183/64‑1 (MQ=255)
                                      aTACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGcc                           <  1:1615065/71‑1 (MQ=255)
                                       tACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCa                                         >  1:2766084/1‑56 (MQ=255)
                                       tACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCa                                         >  1:794101/1‑56 (MQ=255)
                                       tACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCg                          >  1:597702/1‑71 (MQ=255)
                                        accaccAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACaa                               >  1:2932782/1‑65 (MQ=255)
                                        accacaAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACaa                               >  1:2416068/1‑65 (MQ=255)
                                            ccaGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTAcc                     <  1:3066508/71‑1 (MQ=255)
                                             caGCTCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcac                                 <  1:1888646/58‑1 (MQ=255)
                                             caGCGCCATTAACAGTCCGAAAACGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTAc                      >  1:289743/1‑69 (MQ=255)
                                             caGCGCCATTAACAGTCCGAAAACGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTAc                      >  1:3011736/1‑69 (MQ=255)
                                                cgcCATTAACAGTCCGAAAGCGGACTATGGGGCAAgg                                                   <  1:3777934/37‑1 (MQ=255)
                                                  ccATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcac                                 <  1:1360420/53‑1 (MQ=255)
                                                  ccATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcac                                 <  1:551948/53‑1 (MQ=255)
                                                  ccATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGcac                                 <  1:3656243/53‑1 (MQ=255)
                                                     ttAACAGTCCGAAAACGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGac            <  1:571710/71‑1 (MQ=255)
                                                     ttAACAGTCCGAAAACGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGac            <  1:177045/71‑1 (MQ=255)
                                                     ttAACAGTCCGAAAACGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGac            <  1:2562808/71‑1 (MQ=255)
                                                      tAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAg                                        <  1:3713120/42‑1 (MQ=255)
                                                      tAACAGTCCGAAAACGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGac            <  1:3611648/70‑1 (MQ=255)
                                                       aaCAGTCCGAAAGCGGACTCTGGGTCAAGGTAGTTAATCAGTACTCACAACGcc                           <  1:556724/54‑1 (MQ=255)
                                                       aaCAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCa                                         >  1:378599/1‑40 (MQ=255)
                                                       aaCAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGc                            >  1:3770419/1‑53 (MQ=255)
                                                       aaCAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGc                            >  1:1470897/1‑53 (MQ=255)
                                                       aaCAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGc                            >  1:883178/1‑53 (MQ=255)
                                                       aaCAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGacac          >  1:1389852/1‑71 (MQ=255)
                                                          aGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGa             >  1:648451/1‑65 (MQ=255)
                                                              cGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCa               >  1:3559696/1‑59 (MQ=255)
                                                              cGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCa               >  1:2809515/1‑59 (MQ=255)
                                                               gAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGACACCAGAATgg  <  1:3803904/71‑1 (MQ=255)
                                                                aaaGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGACAcc         >  1:3812380/1‑63 (MQ=255)
                                                                aaaGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGACACCAGAAt    >  1:1670899/1‑68 (MQ=255)
                                                                   |                                                                  
CGCCAGGCTAATCATCGCCCAGTTGATTACCAGTGCAGATACCACCAGCGCCATTAACAGTCCGAAAGCGGACTCTGGGGCAAGGTAGTTAATCAGTACGCACAACGCCGTTACCAGTGCAGACACCAGAATGG  >  W3110S.gb/120412‑120545

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