Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,093,040 G→A 100% A215T (GCT→ACT)  hisC → histidinol‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,093,0400GA98.6% 116.2 / ‑6.3 44A215T (GCT→ACT) hisChistidinol‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (20/22);  minor base C (1/0);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTTCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGG  >  NC_000913/2092898‑2093178
                                                                                                                                              |                                                                                                                                          
ttttCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTg                                                                                                                                      >  1:310387/1‑149 (MQ=255)
ttttCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTg                                                                                                                                      <  1:719953/149‑1 (MQ=255)
tgttCGCACCCTGCTGGAGTGAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATCTCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAACACTTTTg                                                                                                                                      <  2:135799/147‑1 (MQ=255)
             ctgGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCt                                                                                                                         >  1:782500/1‑149 (MQ=255)
             ctgGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCt                                                                                                                         >  1:546968/1‑149 (MQ=255)
                gAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCg                                                                                                                      <  1:284072/149‑1 (MQ=255)
                gAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCg                                                                                                                      <  1:520211/149‑1 (MQ=255)
                  gTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGtt                                                                                                                    <  1:109937/149‑1 (MQ=255)
                  gTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGtt                                                                                                                    <  1:72841/149‑1 (MQ=255)
                  gTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGCCGAAAACTTTTGCTCTGGCGGGGCTTCGtt                                                                                                                    >  1:531044/1‑149 (MQ=255)
                  gTTAACCCGCGGTAAGGCGATTGTGGTTGACGATGAAGCCTATATCGAGTTTTACCCGCAGACATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCGGGCGGGGCTTCGtt                                                                                                                    >  1:478179/1‑149 (MQ=255)
                        ccGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGGCCACGGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGAt                                                                                                              >  2:74791/1‑149 (MQ=255)
                        ccGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGAt                                                                                                              <  1:505790/149‑1 (MQ=255)
                            ggTAAGGCGAGTGGGGTTGCCGATGCAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCCATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTTCGGATTTAc                                                                                                          <  2:599657/149‑1 (MQ=255)
                             gTAAGGCGCTTGTGGTTGCCGATGACGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTCTGCTCTGGCTGGGCTTCGTTGCGGATTTAcc                                                                                                         <  2:592455/149‑2 (MQ=255)
                                            ttGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaagaa                                                                                          >  1:126222/1‑149 (MQ=255)
                                              gCCGATGTAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGt                                                                                        <  1:399856/149‑1 (MQ=255)
                                              gCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGt                                                                                        <  1:274525/149‑1 (MQ=255)
                                                cGATGAAGCCTATATCGAGTTTTGCACGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGgca                                                                                      >  1:391728/1‑149 (MQ=255)
                                                 gATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGtcat                                                                                     >  2:698746/1‑149 (MQ=255)
                                                       gCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTCCGCACACTGTCGAAAACCTTTGCTCTGGCGGGGCGTCGTTGCGGATTTACGCTGGAAAACGAGGAAGTCATCAACct                                                                               >  2:777390/1‑149 (MQ=255)
                                                         cTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATGTACGCTGGCAAACGAAGAAGTCATCAACctga                                                                             >  2:439231/1‑148 (MQ=255)
                                                                   ttttGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCCTCAACCTGCTGATGCAAGt                                                                   >  2:672279/1‑149 (MQ=255)
                                                                     ttGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAACGTGa                                                                 >  2:600766/1‑149 (MQ=255)
                                                                     ttGCCCCCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGACGTCATCAACCTGCTGATGAAAGTGa                                                                 >  2:259541/1‑149 (MQ=255)
                                                                      tGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCAACTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCCGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGAt                                                                <  1:71987/149‑1 (MQ=255)
                                                                       gCCCTCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATc                                                               <  2:734612/149‑1 (MQ=255)
                                                                       gCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATc                                                               <  1:605197/149‑1 (MQ=255)
                                                                        cccGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCg                                                              >  2:346137/1‑149 (MQ=255)
                                                                         ccGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTCTTCTACGCAAAATGTCGAAAACTTTTGCTATGGCGGGGCTTCGCTGCGGATTTACTCTGGCGAACGATGAAGTCAGCAACCCGCTGATGACAGGGGTCGc                                                             >  2:264607/1‑149 (MQ=255)
                                                                          cGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCCAACGAAGAAGTCATCAACCTGATGATGAAAGTGATCAcc                                                            >  1:672000/1‑149 (MQ=255)
                                                                                          tggctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGa                                            <  2:83820/149‑1 (MQ=255)
                                                                                               ggctggcGGAATCTCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGTTTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCg                                       <  2:10097/149‑1 (MQ=255)
                                                                                               ggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCg                                       <  1:629422/149‑1 (MQ=255)
                                                                                                  tggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGGCATCAACCTGCTGATGAAAGTGATCGCGCCCTACCCGCTCTCAACGCCGGtt                                    >  2:48649/1‑149 (MQ=255)
                                                                                                     cGGCATATCCGCACCGGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCGGCTCTCGACGCCGGTTGcc                                 <  2:531524/149‑1 (MQ=255)
                                                                                                              cGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGAACGCCCCCTACCCGCTCTCGGCGCCGGTTGCCGACATTGCg                        >  2:199742/1‑149 (MQ=255)
                                                                                                                cACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGc                      <  2:723233/149‑1 (MQ=255)
                                                                                                                   ctgtcgaTTTCACGCACACTGTCGAAGTCTTTTGCTCTGGCGGGGCTTCGTTGCCGATTTACGCTGGCAAACGAACAAGTCATCAACCTGCTGATGCAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCa                   <  2:148562/143‑1 (MQ=255)
                                                                                                                     ggacaTTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGTCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGGTCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAgg                 <  1:339389/145‑1 (MQ=255)
                                                                                                                      gCTATTTTACGCACACTGTCGAAACCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACCTTGCGGCCCAGGc                >  2:125019/1‑149 (MQ=255)
                                                                                                                      gCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATATCCGCTGGCAAACGACGAAGTCATCAAGCTGCTGATGAAAGTGATCGCCCCCTACACGCTCTCGCCGCCGGTAGCCGCCATTGCGGCCCAGGc                >  2:58687/1‑149 (MQ=255)
                                                                                                                            ttACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAg          <  1:305952/149‑1 (MQ=255)
                                                                                                                                    acTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAgg  <  1:698746/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                          
TTTTCGCACCCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGG  >  NC_000913/2092898‑2093178

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