Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 812,947 G→A *226* (TAG→TAA bioD → dethiobiotin synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913812,9470GA97.1% 101.4 / ‑6.2 34*226* (TAG→TAAbioDdethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/16);  total (17/17)
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

GGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812812‑813072
                                                                                                                                       |                                                                                                                             
gggAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGTGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCa                                                                                                                            <  2:1074130/139‑1 (MQ=255)
   aaaaCGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTc                                                                                                                         <  2:370151/139‑1 (MQ=255)
                           gaCGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGAACCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAATCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTAGGTTAAATTgctagcga                                                                                                 >  2:750027/2‑139 (MQ=255)
                           cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagcga                                                                                                 >  2:717456/1‑139 (MQ=255)
                           cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagcga                                                                                                 >  1:844596/1‑139 (MQ=255)
                             cGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCgaga                                                                                               >  2:869832/1‑139 (MQ=255)
                                              cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:844596/139‑1 (MQ=255)
                                              cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              <  2:220617/139‑1 (MQ=255)
                                              cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                              >  2:839377/1‑139 (MQ=255)
                                                               agagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                             <  1:839377/139‑1 (MQ=255)
                                                                                                 gCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTc                           >  2:337299/1‑139 (MQ=255)
                                                                                                  cGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgc                                                                                                       <  2:573092/62‑1 (MQ=255)
                                                                                                  cGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgc                                                                                                       >  1:573092/1‑62 (MQ=255)
                                                                                                        ccGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCa                    <  1:717456/139‑1 (MQ=255)
                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             >  1:604991/1‑64 (MQ=255)
                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             >  1:358041/1‑64 (MQ=255)
                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             <  2:604991/64‑1 (MQ=255)
                                                                                                          ggAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATc                                                                                             <  2:358041/64‑1 (MQ=255)
                                                                                                           gAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           >  1:370162/1‑97 (MQ=255)
                                                                                                           gAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           <  2:370162/97‑1 (MQ=255)
                                                                                                               gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg             <  1:337299/139‑1 (MQ=255)
                                                                                                               gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg             <  2:937585/139‑1 (MQ=255)
                                                                                                                   aTAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAAtg         >  2:618638/1‑139 (MQ=255)
                                                                                                                       aCCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGc                        <  1:237350/120‑1 (MQ=255)
                                                                                                                       aCCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGc                        >  2:237350/1‑120 (MQ=255)
                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           >  1:865309/1‑82 (MQ=255)
                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTaaa                                                           <  2:865309/82‑1 (MQ=255)
                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  >  2:275175/1‑139 (MQ=255)
                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  >  1:262657/1‑139 (MQ=255)
                                                                                                                          ttgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAt  <  1:243843/139‑1 (MQ=255)
                                                                                                                            gccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATcgc                                                                                           <  2:1094479/48‑1 (MQ=255)
                                                                                                                            gccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATcgc                                                                                           >  1:1094479/1‑48 (MQ=255)
                                                                                                                               ttgttgTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 >  1:966359/1‑55 (MQ=255)
                                                                                                                               ttgttgTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAtt                                                                                 <  2:966359/55‑1 (MQ=255)
                                                                                                                                       |                                                                                                                             
GGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAAT  >  NC_000913/812812‑813072

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