Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,244,145 G→T 13.7% F750L (TTC→TTA plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,244,1450GT13.7% 126.0 / 3.9 45F750L (TTC→TTAplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base T (3/3);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

ACACAACACATCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCATCCGGCAATTCTCTCTGATTACCCTTCGCCCTGCGTCGCACTCTCAATCGTCAAACGCACGTCTGATGTAATCAACTCCGCCAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCCAGCAGCTGTAGCGTGCGAGAATGCGCCGGGTTGATATGCAACTCATCATCTTGCAGGGTAATCAGCCCCTGACGTTGCATCTCATTTGCCAG  >  CP009273/4243887‑4244390
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
acacaacacaTCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCATCCGGCAATTCTCTCTGATTACCCTTCGCCCTGCGTCGCACTCTCAATCNTCAAACGCACGGCTGATGTAATCAACTCCGCCAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTg                                                                                                                                                                                                                                        <  1:598737/274‑1 (MQ=255)
 cacaacacaTCGCAATTTATTTAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCATCCGGCAATTCTCTCTGATTACCCTTCGCCCTGCGTCGCACTCTCAATCGTCAAACGCACGTCTGATGTAATCAACTCCGCCAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGACGCTATCGCTGATATACACTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTCGTCCAACAACTCCGGCGCGTTGAt                                                                                                                                                                                                                                      >  2:609644/1‑275 (MQ=255)
                      taaTTTGCAGATTATGGAAGGCCTGATAAGGCGTTTTCGCCGCATCCGGCAATTCTCTCTGATTACCCTTCGCCCTGCGTCGCACTCTCAATCGTCAAACGCACGTCTGATGTAATNAACTCCGCCAGTAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACg                                                                                                                                                                                                                 <  1:609655/274‑1 (MQ=255)
                                                                                                                     cAAACGCACGTCTGATGTAATCAACTCCGCCAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACGCAACAACCAGAAGGTGATAGCATAAAGTTGCa                                                                                                                     >  1:331024/1‑272 (MQ=255)
                                                                                                                           cACGTCTGATGTAATCAACTCCGCCAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAATAACTCCGGCGCGTTGATGCCGTGCAGCACg                                                                                                                                                                                                                         >  1:343506/1‑166 (MQ=255)
                                                                                                                           cACGTCTGATGTAATCAACTCCGCCAGAAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATTACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCAAAGAACTCCGGCGCGTTGATGc                                                                                                                                                                                                                                    >  1:132181/1‑155 (MQ=255)
                                                                                                                               tCTGATGTAATCAACTCCGCCAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAATAACTCCGGCGCGTTGATGCCGTGCAGCACGgaga                                                                                                                                                                                                                     <  2:343501/166‑1 (MQ=255)
                                                                                                                               tCTGATGTAATCAACTCCGCCAGCAACTGAAAAACCATCATCGTTTCAGCCGGTTCTGCATCGCCGCTATCGCTGATATACCCGTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTg                                                                                                                                                                                                                                <  2:132180/155‑1 (MQ=255)
                                                                                                                                                  ccAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTgc                                                                                                                                                                                                             >  2:67043/1‑155 (MQ=255)
                                                                                                                                                  ccAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCAACAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCCCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGAAGGGTTCCGGGCCCGCCTGGGGCCaacagccgg                                                                                     >  1:310909/1‑272 (MQ=255)
                                                                                                                                                    acacacggctAAGCCTTGATCATTGCCGCCGTTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCCGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACGCCGGCGCGTTGAGGGCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCCAGCAGCTGTa                                                                                    <  2:124473/265‑1 (MQ=255)
                                                                                                                                                      cAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGa                                                                                                                                                                                                         <  1:67043/155‑1 (MQ=255)
                                                                                                                                                                           ttCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTg                                                                                                                                                                                                    >  1:166951/1‑139 (MQ=255)
                                                                                                                                                                               gCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGc                                                                                                                                                                                                <  2:166948/139‑1 (MQ=255)
                                                                                                                                                                                              ccgcTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAATAACTCCGGCGCGTTGATGCCGTGCAGCACTGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTg                                                                                                                                                          >  2:135074/1‑162 (MQ=255)
                                                                                                                                                                                                  tATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAATAACTCCGGCGCGTTGATGCCGTGCAGCACTGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCAc                                                                                                                                                      <  1:135075/162‑1 (MQ=255)
                                                                                                                                                                                                              aCCCTTCATCACGCAGTTGCAGCACCAGCGAACTTAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCCAGCAGCTGTAGCGTGCGAGAATGCGCCGGGTTGATATGCAACTCATCATCTTGCAGGGTAATCAGcccc                         <  1:225759/275‑1 (MQ=255)
                                                                                                                                                                                                                 cTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGc                                                                                                                                                        >  2:231608/1‑145 (MQ=255)
                                                                                                                                                                                                                  ttCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGc                                                                                                                                                                                                                                    >  1:128557/1‑68 (MQ=255)
                                                                                                                                                                                                                  ttCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGAc                                                                                                                                                                >  2:320463/1‑136 (MQ=255)
                                                                                                                                                                                                                    catcaCGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgc                                                                                                                                                                                                                               >  1:328651/1‑71 (MQ=255)
                                                                                                                                                                                                                    catcaCGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCCACAACCAGAAGGGGCTGGCATAACGTTGCAGCGTTTCCCGCGCGCCGGc                                                                                                  >  2:621093/1‑196 (MQ=255)
                                                                                                                                                                                                                     atcaCGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTc                                                                                                                                                    <  1:231612/145‑1 (MQ=255)
                                                                                                                                                                                                                      tcaCGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTg                                                                                                                                                                                                                                <  2:128556/68‑1 (MQ=255)
                                                                                                                                                                                                                      tcaCGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGt                                                                                                                                                            <  1:320468/136‑1 (MQ=255)
                                                                                                                                                                                                                        tcGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGTACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACAAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGcc                                                                                              <  1:621104/195‑1 (MQ=255)
                                                                                                                                                                                                                        aCGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgcagca                                                                                                                                                                                                                           <  2:328646/71‑1 (MQ=255)
                                                                                                                                                                                                                         cGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTgc                                                                                                                                                                                                             >  2:617174/1‑84 (MQ=255)
                                                                                                                                                                                                                             gtgtCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGa                                                                                                                                                                                                         <  1:617185/84‑1 (MQ=255)
                                                                                                                                                                                                                             gtgtCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCCAGCAGCTGTAGCGTGCGAGAATGCGCCGGGTTg                                                             <  1:307057/224‑1 (MQ=255)
                                                                                                                                                                                                                              tgtCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGtt                                                                                                                                                                                                                                         >  2:486137/1‑51 (MQ=255)
                                                                                                                                                                                                                              tgtCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTg                                                                                                                                                                                                    >  2:603882/1‑88 (MQ=255)
                                                                                                                                                                                                                              tgtCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGAt                                                                                                                                   >  2:326144/1‑153 (MQ=255)
                                                                                                                                                                                                                                  aGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATg                                                                                                                                                                                                                                     <  1:486145/51‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGc                                                                                                                                                                                                <  1:603893/88‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCa                                                                                                                               <  1:326149/153‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaagaACGGAANCCCGCCTTGGAGAAGGACCCCGGCGCGTTGATGCCGTGCAGGCCGGGGAGACGTTTCGCGACGGTGGGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCAGAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCCAGCAGCTGTAGCGTGCGAGAATGCGCCGGGTTGATATGCAACTCATCATCTTGCAGGGTAATCAGCCCCTGACGTTGCATCTCATTTGCCAg  <  2:501065/270‑1 (MQ=255)
                                                                                                                                                                                                                                        agagAACTGAACACCGCCTTGTCGAAGATCTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCCAGCNGCTGTAGCGTGCGAGCATGCGCCGGGTTGATATGCAACTCATCATCTTGCAGGGTAATCAGCCCCTGACGTTGCCTCTCATTTg      >  1:501074/1‑268 (MQ=255)
                                                                                                                                                                                                                                        agagAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCaacaac                                                                                                                                              >  2:397390/1‑132 (MQ=255)
                                                                                                                                                                                                                                        agagAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCAt                                                                                                                              >  2:597410/1‑148 (MQ=255)
                                                                                                                                                                                                                                          agaACTGAACACCGCCTTGTCGAATAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTtc                                                                                                                                                                                        >  1:92510/1‑88 (MQ=255)
                                                                                                                                                                                                                                            aaCTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGa                                                                                                                                          <  1:397395/132‑1 (MQ=255)
                                                                                                                                                                                                                                            aaCTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACg                                                                                                                          <  1:597421/148‑1 (MQ=255)
                                                                                                                                                                                                                                              cTGAACACCGCCTTGTCGAATAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCa                                                                                                                                                                                    <  2:92510/88‑1 (MQ=255)
                                                                                                                                                                                                                                                           tGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCgg                                                                                                >  1:31763/1‑159 (MQ=255)
                                                                                                                                                                                                                                                               gaagaaCTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCcag                                                                                            <  2:31763/159‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gaaCTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGAc                                                                                                                                                                >  1:178291/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     
ACACAACACATCGCAATTTATTGAATTTGCAGATTATGGAAGGCCGGATAAGGCGTTTTCGCCGCATCCGGCAATTCTCTCTGATTACCCTTCGCCCTGCGTCGCACTCTCAATCGTCAAACGCACGTCTGATGTAATCAACTCCGCCAGCAACTGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCGCGGTTGATCGACGGGTTGGCACTCAACAACCAGAAGGTGATGGCATAACGTTGCAGCGTTTCGCGCGCGCCTGCGGCCAGCAGCTGTAGCGTGCGAGAATGCGCCGGGTTGATATGCAACTCATCATCTTGCAGGGTAATCAGCCCCTGACGTTGCATCTCATTTGCCAG  >  CP009273/4243887‑4244390

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