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,9470GA98.0% 156.8 / NA 50*226* (TAG→TAAbioDdethiobiotin synthetase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/26);  total (23/27)
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

CGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAG  >  NC_000913/812811‑813079
                                                                                                                                        |                                                                                                                                    
cGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAAcc                                                                                                                                    >  2:1416645/1‑139 (MQ=255)
   gAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                                 <  2:1617259/139‑1 (MQ=255)
          tCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAtt                                                                                                                          <  1:961602/139‑1 (MQ=255)
              gCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGt                                                                                                                      <  2:1023362/139‑1 (MQ=255)
                    tatatGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAAtt                                                                                                                <  1:834776/139‑1 (MQ=255)
                            cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                                 <  1:1582709/114‑1 (MQ=255)
                            cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCAtt                                                                                                                                 >  2:1582709/1‑114 (MQ=255)
                            cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                                            <  2:860641/119‑1 (MQ=255)
                            cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                                            >  1:860641/1‑119 (MQ=255)
                            cACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgagcga                                                                                                        >  1:512959/1‑139 (MQ=255)
                              cGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCgaga                                                                                                      >  1:232557/1‑139 (MQ=255)
                                   aCCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                                            <  1:1536816/112‑1 (MQ=255)
                                   aCCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTa                                                                                                                            >  2:1536816/1‑112 (MQ=255)
                                    cccGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGt                                                                                                >  2:360074/1‑139 (MQ=255)
                                       gCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCtt                                                                                             <  1:1561216/139‑1 (MQ=255)
                                           gATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGAt                                                                                         <  2:130078/139‑1 (MQ=255)
                                             ttCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTg                                                                                       >  1:1376382/1‑139 (MQ=255)
                                               cccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGAc                                                                                     >  2:939555/1‑139 (MQ=255)
                                                ccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACt                                                                                    <  1:360074/139‑1 (MQ=255)
                                                ccGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTTTTGGTTATATTGCGAGCGAGATCGCGTCTTCGATTGACTg                                                                                   <  2:983648/139‑1 (MQ=255)
                                                  gcgcCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGc                                                                                  >  1:566195/1‑139 (MQ=255)
                                                        ctgctgGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTa                                                                            <  2:1376382/139‑1 (MQ=255)
                                                               gagagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAAtt                                                                     <  1:939555/139‑1 (MQ=255)
                                                                agagaTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTa                                                                    <  1:727582/139‑1 (MQ=255)
                                                                    aTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAtt                                                                                                                          >  1:1189410/1‑81 (MQ=255)
                                                                    aTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAtt                                                                                                                          <  2:1189410/81‑1 (MQ=255)
                                                                    aTCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAAtt                                                                <  1:180432/139‑1 (MQ=255)
                                                                     tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgag                                                                                                           <  2:1360987/95‑1 (MQ=255)
                                                                     tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTgcgag                                                                                                           >  1:1360987/1‑95 (MQ=255)
                                                                     tCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTc                                                               >  2:1119566/1‑139 (MQ=255)
                                                                              ttGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAt                                                                                                                           >  1:258288/1‑70 (MQ=255)
                                                                              ttGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTAt                                                                                                                           <  2:258288/70‑1 (MQ=255)
                                                                                 cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTa                                                                            >  1:961509/1‑114 (MQ=255)
                                                                                 cAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATGTa                                                                            <  2:961509/114‑1 (MQ=255)
                                                                                           aGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAAtt                                         >  2:1355176/1‑139 (MQ=255)
                                                                                           aGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAAtt                                         >  2:233121/1‑139 (MQ=255)
                                                                                                   cGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTcc                                 <  1:1572927/139‑1 (MQ=255)
                                                                                                           ggAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTGAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATg                         <  1:1119566/139‑1 (MQ=255)
                                                                                                                gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg                    <  1:1355176/139‑1 (MQ=255)
                                                                                                                gTACATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcggcg                    <  2:232557/139‑1 (MQ=255)
                                                                                                                   cATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGtt                                                                                                                     <  2:441227/39‑1 (MQ=255)
                                                                                                                   cATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGtt                                                                                                                     >  1:441227/1‑39 (MQ=255)
                                                                                                                   cATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAAtt                                                                                                                >  1:1251167/1‑44 (MQ=255)
                                                                                                                   cATAAACCTTGCCTTGTTGTAACCATTCTGTATTTGGTTAAAtt                                                                                                                <  2:1251167/44‑1 (MQ=255)
                                                                                                                            tgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcg                       <  2:1602146/124‑1 (MQ=255)
                                                                                                                            tgccttgTTGTAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATgcg                       >  1:1602146/1‑124 (MQ=255)
                                                                                                                                 tgttgtAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTa                                                                            >  2:1330756/1‑66 (MQ=255)
                                                                                                                                 tgttgtAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTa                                                                            <  1:1330756/66‑1 (MQ=255)
                                                                                                                                  gttgtAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAg  >  2:203648/1‑139 (MQ=255)
                                                                                                                                  gttgtAACCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAg  >  1:477684/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                    
CGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAAAGTACATAAACCTTGCCTTGTTGTAGCCATTCTGTATTTGGTTAAATTGCGAGCGAGATCGCGTCTTCGATTGACTGCAATTTAACCAATTAAATTCTAAAATAATCACGAAAAAAATTTTACTTCCGCCTCATGCGGCGAATGTGGGAATTGCCCAG  >  NC_000913/812811‑813079

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