Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 437,107 G→A *173* (TAG→TAA pgpA → phosphatidylglycerophosphatase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913437,1070GA98.0% 157.9 / ‑4.4 50*173* (TAG→TAApgpAphosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/27);  minor base T (0/1);  total (22/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCA  >  NC_000913/436973‑437240
                                                                                                                                      |                                                                                                                                     
tGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATtgt                                                                                                                                   <  1:561442/139‑1 (MQ=255)
   ccGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCa                                                                                                                                <  1:1480288/139‑1 (MQ=255)
      aTCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCg                                                                                                                             >  1:178061/1‑139 (MQ=255)
        ccGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGAt                                                                                                                           <  2:1705027/139‑1 (MQ=255)
          gCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGc                                                                                                                         <  2:1384748/139‑1 (MQ=255)
                                      cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTATTTGTGCAcc                                                                                                                              <  2:1230528/106‑1 (MQ=255)
                                      cATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCAcc                                                                                                                              >  1:1230528/1‑106 (MQ=255)
                                        tGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGAt                                                                                                                           <  1:1477322/107‑1 (MQ=255)
                                        tGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGAt                                                                                                                           >  2:1477322/1‑107 (MQ=255)
                                         ggggATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGc                                                                                          <  1:345018/139‑1 (MQ=255)
                                         ggggATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGc                                                                                          <  1:329574/139‑1 (MQ=255)
                                                      gacgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTAtt                                                                             >  2:1733804/1‑139 (MQ=255)
                                                       acgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  1:817224/1‑120 (MQ=255)
                                                       acgaTATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  2:817224/120‑1 (MQ=255)
                                                          atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGGCACCGATGCCTGATGTGACGCTTGTc                                                                                                        <  2:895147/107‑1 (MQ=255)
                                                          atatTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGGCACCGATGCCTGATGTGACGCTTGTc                                                                                                        >  1:895147/1‑107 (MQ=255)
                                                            atTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATc                                                                       <  2:178061/139‑1 (MQ=255)
                                                              tGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         <  1:737840/103‑1 (MQ=255)
                                                              tGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGt                                                                                                         >  2:737840/1‑103 (MQ=255)
                                                               gTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTa                                                                    >  1:1518293/1‑139 (MQ=255)
                                                                  gCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTCATTGTGCACCGATGc                                                                                                                         <  1:1225609/83‑1 (MQ=255)
                                                                  gCCGGGGTGATTTCCGCAGGCATCCTGTAGTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGTACCGATGc                                                                                                                         >  2:1225609/1‑83 (MQ=255)
                                                                   ccGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGc                                                                                                                         <  1:576486/82‑1 (MQ=255)
                                                                   ccGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGc                                                                                                                         >  2:576486/1‑82 (MQ=255)
                                                                       ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTa                                                                           >  1:181871/1‑124 (MQ=255)
                                                                       ggTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTa                                                                           <  2:181871/124‑1 (MQ=255)
                                                                        gTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGAt                                                           <  1:1733804/139‑1 (MQ=255)
                                                                             ttCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGt                                                      <  1:1428953/139‑1 (MQ=255)
                                                                                        tCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGt                                           >  1:1581082/1‑139 (MQ=255)
                                                                                                  tATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCGTCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCg                                 <  2:1459367/139‑1 (MQ=255)
                                                                                                   aTCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCgg                                >  1:1367873/1‑139 (MQ=255)
                                                                                                    tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                    >  1:271914/1‑70 (MQ=255)
                                                                                                    tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                    <  2:271914/70‑1 (MQ=255)
                                                                                                    tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGc                                                         <  2:1539522/113‑1 (MQ=255)
                                                                                                    tCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGc                                                         >  1:1539522/1‑113 (MQ=255)
                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCTTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  2:1226023/72‑1 (MQ=255)
                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  1:904264/1‑72 (MQ=255)
                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               >  1:1226023/1‑72 (MQ=255)
                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCtcat                                                                                               <  2:904264/72‑1 (MQ=255)
                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAct                                                                                   >  2:317246/1‑84 (MQ=255)
                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGAct                                                                                   <  1:317246/84‑1 (MQ=255)
                                                                                                       gTCATCACTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAAt                            >  1:1002545/1‑139 (MQ=255)
                                                                                                             aCTGGCCGCTGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTc                      <  2:1581082/139‑1 (MQ=255)
                                                                                                                      tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCaa             >  1:468443/1‑139 (MQ=255)
                                                                                                                      tGGGTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCaa             >  1:14761/1‑139 (MQ=255)
                                                                                                                       gggTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCaaa            <  2:1002545/139‑1 (MQ=255)
                                                                                                                        ggTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                    >  1:174329/1‑50 (MQ=255)
                                                                                                                        ggTATTCTGTCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGt                                                                                                    <  2:174329/50‑1 (MQ=255)
                                                                                                                                 tCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCa  <  2:468443/139‑1 (MQ=255)
                                                                                                                                 tCGTAATTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCa  <  1:608369/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                     
TGGCCGATCCGCTGGTTTGATCGCAATGTGCATGGCGGCATGGGGATCATGATCGACGATATTGTCGCCGGGGTGATTTCCGCAGGCATCCTGTATTTTATCGGTCATCACTGGCCGCTGGGTATTCTGTCGTAGTTGTGCACCGATGCCTGATGTGACGCTTGTCACGTCTCATCAGGCCTGGACTCTTATTTAAATCCTACGACAGGATGCGGTTTATACGGCGTTTCCAGTTCGGCAATCTGTTCCGGCTTCAAAGTGATATCCA  >  NC_000913/436973‑437240

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