Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,265,166 G→A 7.0% A288V (GCA→GTA)  tdcB ← L‑threonine dehydratase, catabolic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,265,1660GA 7.0% 188.5 / 4.4 57A288V (GCA→GTA) tdcBL‑threonine dehydratase, catabolic
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

CGAATTAATTAAGCGTCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGGCG  >  NC_000913/3265031‑3265305
                                                                                                                                       |                                                                                                                                           
cGAATTAATTAAGCGTCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCag                                                                                                                                         <  1:813371/140‑1 (MQ=255)
cGAATTAATTAAGCGTCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCag                                                                                                                                         <  1:812411/140‑1 (MQ=255)
cGAATTAATTAAGCGTCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCag                                                                                                                                         <  1:813367/140‑1 (MQ=255)
             cGTCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAgcgcctgcgcct                                                                                                                            >  1:1612913/1‑140 (MQ=255)
               tCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTc                                                                                                                          >  1:337917/1‑140 (MQ=255)
               tCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGAGCCTGCGCCTTc                                                                                                                          >  1:1186758/1‑140 (MQ=255)
                 aaCGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCgg                                                                                                                        >  1:1219795/1‑140 (MQ=255)
                    gAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTga                                                                                                                     <  2:993870/140‑1 (MQ=255)
                     aaaCCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTgac                                                                                                                    <  1:1157335/140‑1 (MQ=255)
                     aaaCCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTgac                                                                                                                    <  1:414032/140‑1 (MQ=255)
                       aCCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTgacga                                                                                                                  >  2:1155732/1‑140 (MQ=255)
                       aCCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTgacga                                                                                                                  <  1:547257/140‑1 (MQ=255)
                                 tGAGAGACGCGAGAAAGATCTATATTGCCGCCGGTATTTATATTGATGGTTTTTTTGTTTTGAATATATTGGTCTAATTTACCGCTTTATATTGCAGCTGATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTgcg                                                                                                        <  2:374125/140‑1 (MQ=255)
                                 tGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAACTAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTgcg                                                                                                        <  1:1072692/140‑1 (MQ=255)
                                 tGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTgcg                                                                                                        <  1:1620433/140‑1 (MQ=255)
                                 tGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTgcg                                                                                                        <  1:486120/140‑1 (MQ=255)
                                 tGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTgcg                                                                                                        <  2:231004/140‑1 (MQ=255)
                                 tGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAAAAATACTCACGGAATTTCTGAATTGAATATATTGCTCTAATTTACCGCTTAATAAAGCAGCACATGCCAGAGCGCCTGCGCCATCGGTGACGACATTATTgcg                                                                                                        <  2:232951/140‑1 (MQ=255)
                                  gagagaCGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTgcgc                                                                                                       <  2:1096126/140‑1 (MQ=255)
                                         gctataaaGATCGATAATGACTCCTTAAATAAAACATAAGGTTTTACTGTTTTGAATATATATTTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTTCGCCTTCTGTGACGACTTTATTGCGCTGAATTa                                                                                                >  2:1469617/6‑140 (MQ=255)
                                          cgAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTaa                                                                                               <  1:946/140‑1 (MQ=255)
                                          cgAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTaa                                                                                               <  2:357816/140‑1 (MQ=255)
                                          cgAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTaa                                                                                               <  2:95613/140‑1 (MQ=255)
                                            agatagatCGTTATTGCCGCCGGAAATAATACTGAGGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAgg                                                                                             <  1:812059/140‑1 (MQ=255)
                                            agaAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAgg                                                                                             <  1:284185/140‑1 (MQ=255)
                                            agaAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAgg                                                                                             <  1:812040/140‑1 (MQ=255)
                                             gaAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGc                                                                                            <  2:237788/140‑1 (MQ=255)
                                             gaAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGc                                                                                            <  2:236166/140‑1 (MQ=255)
                                                    cGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATa                                                                                     <  1:1166058/140‑1 (MQ=255)
                                                     gATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATAc                                                                                    >  1:395230/1‑140 (MQ=255)
                                                            ccgccgGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATACCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCt                                                                             >  1:1013985/1‑140 (MQ=255)
                                                            ccgccgGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATACCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCt                                                                             >  1:1013341/1‑140 (MQ=255)
                                                                 ggAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGAttt                                                                        >  2:364351/1‑140 (MQ=255)
                                                                 ggAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGAttt                                                                        >  2:364355/1‑140 (MQ=255)
                                                                    aataatACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGt                                                                     <  1:191149/140‑1 (MQ=255)
                                                                         taCTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCg                                                                <  1:1013769/140‑1 (MQ=255)
                                                                              aCGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGAc                                                           >  1:320711/1‑140 (MQ=255)
                                                                              aCGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGAc                                                           <  1:75136/140‑1 (MQ=255)
                                                                                 gTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAg                                                        >  2:97360/1‑140 (MQ=255)
                                                                                    tttCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCAc                                                     >  1:143087/1‑140 (MQ=255)
                                                                                    ttacTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCACTCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCAc                                                     >  1:142972/4‑140 (MQ=255)
                                                                                            ttGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCAt                                             >  1:1043045/1‑140 (MQ=255)
                                                                                                  atatTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTa                                       >  1:37294/1‑140 (MQ=255)
                                                                                                         tCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCacg                                >  2:1278684/1‑140 (MQ=255)
                                                                                                         tCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCacg                                >  2:1278712/1‑140 (MQ=255)
                                                                                                         tCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCacg                                >  2:1278706/1‑140 (MQ=255)
                                                                                                                   cGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCg                      >  1:884814/1‑140 (MQ=255)
                                                                                                                     cTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTa                    >  2:1316587/1‑140 (MQ=255)
                                                                                                                      ttAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTaa                   <  2:344728/140‑1 (MQ=255)
                                                                                                                      ttAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTaa                   <  2:437043/140‑1 (MQ=255)
                                                                                                                      ttAATAATGCAGCACATACCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTaa                   <  2:1013341/140‑1 (MQ=255)
                                                                                                                      ttAATAATGCAGCACATACCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTaa                   <  2:1013985/140‑1 (MQ=255)
                                                                                                                       taataaTGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAg                  >  2:441291/1‑140 (MQ=255)
                                                                                                                             tGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTTAGTTAAAt            >  1:567586/1‑140 (MQ=255)
                                                                                                                               cagcaCATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTa          >  2:1016905/1‑140 (MQ=255)
                                                                                                                                  cacaTGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCacgaa                              >  1:941769/1‑117 (MQ=255)
                                                                                                                                  cacaTGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCacgaa                              <  2:941769/117‑1 (MQ=255)
                                                                                                                                       gCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCcgggcg  >  2:416374/1‑140 (MQ=255)
                                                                                                                                       gCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCcgggcg  >  2:417775/1‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                           
CGAATTAATTAAGCGTCAACGAAACCGGTGATTTGAGAGACGCGAGAAAGATCGATATTGCCGCCGGAAATAATACTGACGGTTTTTCTGTTTTGAATATATTGGTCTAATTTACCGCTTAATAATGCAGCACATGCCAGAGCGCCTGCGCCTTCGGTGACGACTTTATTGCGCTGAATTAAGGCAATCATACTGTTTCTGATTTCGTCTTCGCTGACCAGCACGATGTCATCGACTAATTCACGAACGATTTCGTAAGTTAAATTACCCGGGCG  >  NC_000913/3265031‑3265305

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: