Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 808,205 A→C 22.1% R351R (CGT→CGG)  bioA ← 7,8‑diaminopelargonic acid synthase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913808,2050AC22.1% 67.1 / 3.2 45R351R (CGT→CGG) bioA7,8‑diaminopelargonic acid synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (13/21);  new base C (5/6);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATCAGGTAAATCAGTTTGCCAAAAGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCA  >  NC_000913/808057‑808351
                                                                                                                                                    |                                                                                                                                                  
gcATCAGGTAAATCAGTTTGCCAAAAGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCa                                                                                                                                                    >  2:405224/1‑149 (MQ=255)
  aTCAGGTAAATCAGTTTGCCAAAAGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCAGATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCCTCCcg                                                                                                                                                  >  1:484962/1‑149 (MQ=255)
        tAAATCAGTTTGCCAAAAGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCAcgggcggg                                                                                                                                            <  1:379904/149‑1 (MQ=255)
        tAAATCAGTTTGCCAAAAGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCAcgggcggg                                                                                                                                            <  1:461052/149‑1 (MQ=255)
           atcatgttGCCAAAAGGCGGGATCCAGACACCCTGTTCGACATAGTATTTTTGCAGCGCCGCAATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACGATTTCGGCATCACGGGCGGGGGc                                                                                                                                         <  2:1001756/143‑1 (MQ=255)
              aGTTTGCCAAAAGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAg                                                                                                                                      <  2:245858/149‑1 (MQ=255)
                      aaaaGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCCGCGCCGCCATATTCCCCGGATGAGTGGTTTCGACCACGCCACTGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCCCGGGCGGGAGCAAGTTCCTctc                                                                                                                              >  2:711870/1‑147 (MQ=255)
                          gggcGGCTCCAGACACCCTGAGCGACTAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTGCGACCACGCCCATGGCCCCCAGTACGCGCACATCGTCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAg                                                                                                                          <  1:915676/146‑1 (MQ=255)
                            ccGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCTGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGCCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCt                                                                                                                        <  2:544843/149‑1 (MQ=255)
                            ccGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCt                                                                                                                        <  1:653115/149‑1 (MQ=255)
                              ggATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCACCCATTTCGGCAGCCCGGGCGGGGGCAAGTTGCTCGCGCAGCTGt                                                                                                                      >  1:111834/1‑149 (MQ=255)
                                aTCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTAc                                                                                                                    <  1:132693/149‑1 (MQ=255)
                                    aGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCa                                                                                                                <  1:970592/149‑1 (MQ=255)
                                      acacCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCACCCATTTCGGCATCACGGGCGGGGTCAAGTTGCTCGCGCAGCTGTACTTCAat                                                                                                              >  1:610732/1‑149 (MQ=255)
                                      acacCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAat                                                                                                              <  1:496040/149‑1 (MQ=255)
                                        acCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTGACCGGATGAGTGGTTTCGACCACGCAAATGGCCCGCAGTACGCGCACATCGGCCACCATATCGGCAGCCCGGGGGGGGGCAAGTTGCTCGGGCTGCGGGAATTCAagat                                                                                                            >  2:284918/1‑149 (MQ=255)
                                         cccTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCAATTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATc                                                                                                           <  1:405224/149‑1 (MQ=255)
                                              ttCGACAAAGAATTTTTGCAGCGCCGCCATCTTCACCGGATGAGGGGGTTCGACCACGCCAATGGCCCCCAGGACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCTAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCa                                                                                                      <  1:966524/149‑1 (MQ=255)
                                                 gACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCt                                                                                                   >  1:528937/1‑149 (MQ=255)
                                                    aaaGAATTTTTGCAGCGCCGCCATATTCACCGGATGCGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCCCGGGCGGGGGCAAGTTGCTCGCGCAGCGGTACTTCAATATCCGCCCCCTGtt                                                                                                >  1:719706/1‑149 (MQ=255)
                                                    aaaGAATTTTTGCAGCGCCGCCAGATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACCTCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGGGCAGCTGTACTTCAATATCCGCCACCTGtt                                                                                                >  1:716014/1‑149 (MQ=255)
                                                                 aGCGCCGCCATATTCACCGGATGAGTGGGTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGcc                                                                                   >  2:737643/1‑149 (MQ=255)
                                                                        ccATATTCACCGGATGAGTGGAGTCGACCACGCCAATCGCCCCCAGCACGCGCACATCCGCAACCATTTCGGCATCCCGGGCGGGGGCAAGGAGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCg                                                                            <  1:584290/149‑1 (MQ=255)
                                                                              tCACCGGATGCGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCAACTGTTGCTGCCAGTCGCCAGATTCGAGAATc                                                                      <  2:1013428/149‑1 (MQ=255)
                                                                                aCCGGATGAGTGGGTTCGACCACGCCAATCGCCCCCCGTACGCGCACATCGGCAACCATTTCGGCAGCCCGGCCGGGGGCAAGTTGCTCGCGCAGCTGGACTTCAATATCCGCCCCCTGTGGCTGCCAGTCGCCAGATTCGAGAATCGc                                                                    <  2:328678/149‑1 (MQ=255)
                                                                                            gTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGtt                                                        <  1:870446/149‑1 (MQ=255)
                                                                                            gTTTCGACCACGCCAATGGCCCCAAGTACGCGCACATCGGAAACCATTTCGGCAGAACGGGCGGGGGCAAGTAGCTCGCGCAGATCTACGTCAATATCCGCCAACTGTAGCTGCGAGTCGCGAGGATGGAGCAGCGACCGGCTAGCGtt                                                        >  1:278084/1‑149 (MQ=255)
                                                                                             tttCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGttt                                                       <  1:657047/149‑1 (MQ=255)
                                                                                                cGACCCCGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGGACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATGCGAGAATCGCCAGGCTGGCGTTtgct                                                    >  2:220303/1‑149 (MQ=255)
                                                                                                  aCCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTtgctgc                                                  >  1:45985/1‑149 (MQ=255)
                                                                                                   ccACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGcc                                                 >  1:941410/1‑149 (MQ=255)
                                                                                                   ccACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGAAGCTGTACTTCAATATCCGCCCCCTGTTGATGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGcc                                                 >  2:572569/1‑149 (MQ=255)
                                                                                                    cACGCCAATGGCCCCCAGTACGCGCACATCCGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCATTATCCTCCACCTGTTGCTGCCAGTCGCCAGATTCGAGATTAGCCTGGCTGGCGTTtgctgcgg                                                >  1:961525/1‑147 (MQ=255)
                                                                                                          aagggCCCCCCGTACTCGCACAGCGGCAACCATATCGGCATCACGGCCGGGGGCAAGTTGCTCGCGCAGCGGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCGGGCGTTTGCTGCCGCGCAgg                                          <  1:171560/146‑1 (MQ=255)
                                                                                                          aaTGGCCCCCAGTACGCGCACATAGGCATCCATGGCGGCATCCCGGGCGGGGGCAAGGGGCTCGCGCAGCTGTACTTCAATATGCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAgg                                          <  1:584107/149‑1 (MQ=255)
                                                                                                                  ccAGTACGCGCACCTCGCCAACCATTTCGGCATCCCGGGCGGGGGCAAGTTGCTCGCGCAGCGGTACTTCAATATCCGCCACCTGTTGCTGCCAGGCGCGAGATTCGAGAATCGCCAGGCTGGCGGTTGCTGCCGCGCAGGCCAGCGGa                                  <  1:131143/149‑1 (MQ=255)
                                                                                                                      tACGCGCACAGCGGCAACCATTGCGGCATCACGGGCGGGGGCAAGTGGCGCGCGCAGCTGAACTTCAATATCCGCCCCCTGTGGCTGCCAGTCGCCACATTAGCGACTCGCCAGGCGGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGc                              <  1:711870/149‑1 (MQ=255)
                                                                                                                        cgcgCACATCGGCAACCATTTCGGCATCCCGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGccc                            <  1:513711/149‑1 (MQ=255)
                                                                                                                           gcACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAg                                           >  2:198067/1‑131 (MQ=255)
                                                                                                                           gcACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAg                                           <  1:198067/131‑1 (MQ=255)
                                                                                                                               aTCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCAAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGACCATaatat                     >  2:956677/1‑146 (MQ=255)
                                                                                                                                  ggCAACCATATCGGCATACCGGGCGGGGGCCAGTTGCGCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGGCAGTCGCCAGATTCGAGAAGCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTg                  <  1:1024737/149‑1 (MQ=255)
                                                                                                                                         aactCGGCAGCACCGGCGCGGGCAAGTTGCTCGCGCAGCTGTACTTCAAAAACCGCCACCTGGAGCTGCCAGTCGACAGATTCGACAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCATCGGATTGCCCATAAAAGTTGGCCCATg           <  2:484962/146‑1 (MQ=255)
                                                                                                                                               gCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTCCTTCAATATCCGCCCCCTGTTGCCGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGACAGCGGATTGCCAATAAAAGTTGGCCCATGCATaaa     <  1:757521/149‑1 (MQ=255)
                                                                                                                                               gCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATaaa     <  1:481795/149‑1 (MQ=255)
                                                                                                                                                  tCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCa  >  1:899409/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCATCAGGTAAATCAGTTTGCCAAAAGGCCGGATCCAGACACCCTGTTCGACAAAGAATTTTTGCAGCGCCGCCATATTCACCGGATGAGTGGTTTCGACCACGCCAATGGCCCCCAGTACGCGCACATCGGCAACCATTTCGGCATCACGGGCGGGGGCAAGTTGCTCGCGCAGCTGTACTTCAATATCCGCCACCTGTTGCTGCCAGTCGCCAGATTCGAGAATCGCCAGGCTGGCGTTTGCTGCCGCGCAGGCCAGCGGATTGCCCATAAAAGTTGGCCCATGCATAAAGCA  >  NC_000913/808057‑808351

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