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.5% 101.0 / ‑6.6 37A215T (GCT→ACT) hisChistidinol‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/18);  minor base C (0/2);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

CCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGAAT  >  NC_000913/2092907‑2093182
                                                                                                                                     |                                                                                                                                              
ccTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACCTTTGCTCTGGCgg                                                                                                                                 >  1:668818/1‑149 (MQ=255)
 ctgctgGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCggg                                                                                                                                >  1:159203/1‑149 (MQ=255)
  tgctgGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCgggg                                                                                                                               <  1:546883/149‑1 (MQ=255)
             aCCCGCGGTAAGGCGATTGGTGTTGCCGATGACGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATGCTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCgg                                                                                                                    <  1:5655/149‑1 (MQ=255)
                cgcgtcaaGGCGACTGTGGTAGCCGATGAGGCCTTTAACGGGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGAtt                                                                                                                 <  2:218857/143‑1 (MQ=255)
                     tAAGGCGATTGTGGTTGCCGATGAAGCCTACATCGCGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTCTTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGc                                                                                                            <  1:25489/149‑1 (MQ=255)
                       aGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTg                                                                                                          <  1:552488/149‑1 (MQ=255)
                            aTTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGcgcg                                                                                                     >  2:729458/1‑146 (MQ=255)
                                tgGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaa                                                                                                 <  1:425886/149‑1 (MQ=255)
                                 ggTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaag                                                                                                <  1:138127/149‑1 (MQ=255)
                                 ggTTGCCGATGAAGCCTAGATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTCTTTTACGCACAATGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaag                                                                                                <  2:429992/149‑1 (MQ=255)
                                   ttGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTCACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACgaagaa                                                                                              <  1:339699/149‑1 (MQ=255)
                                    tGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACCCTGGCCAACGgagaag                                                                                             >  2:70893/1‑149 (MQ=255)
                                                       gAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTCCGCTGGCAAACGAAGACGTCATCAACCTGCTGATGaa                                                                          >  1:151566/1‑149 (MQ=255)
                                                       gAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGCTGaa                                                                          >  1:76390/1‑149 (MQ=255)
                                                       gAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGAGGAATATCCGCACCTGGCTATGTTACGCACACTGTCGAAAACTTTTGGTGTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGCAGTCATACACCTGCTGGTGaa                                                                          >  2:661409/1‑149 (MQ=255)
                                                                   cAGGCATCGCTGGCTGGTTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGcccc                                                              <  1:165937/149‑1 (MQ=255)
                                                                       cATCGCTGGCTGGCTGGCTGGCGGAATATCCGCCTGTCGTTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCGTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTAc                                                          <  1:202134/149‑1 (MQ=255)
                                                                             tggctggctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGctc                                                    <  2:78073/149‑1 (MQ=255)
                                                                               gcgagctggctagcGTAAGAACGGCACCAGGCTATTTTACGCACACGGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGACGTACGCTGGCAAACGAAGAAGGCATCAACCGGCTGATGAAAGTGATCGCCCCCTACCCGctccc                                                  >  2:378735/5‑147 (MQ=255)
                                                                                   gctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGAGCGCCCCCTACCCGCTCTCGACg                                              >  1:189078/1‑149 (MQ=255)
                                                                                    ctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGc                                             <  2:346003/149‑1 (MQ=255)
                                                                                    ctggctggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACCAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGc                                             <  1:324855/149‑1 (MQ=255)
                                                                                         tggcGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTCACGGGGCTTCGTTGCGGATTTACGCGAGGGAAGGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCCACGACGGtt                                        >  1:328415/1‑149 (MQ=255)
                                                                                                     cGCACCTGGCTATTTTACGCACACTGTCGAAACCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCg                            <  2:668818/149‑1 (MQ=255)
                                                                                                     cGCACCTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCg                            >  2:655315/1‑149 (MQ=255)
                                                                                                      gacccTGGCTATTTTCCGCACACTGTCGAAACCTTTGGCTCGGGCGGGGCTTCATTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCgg                           <  2:318007/146‑1 (MQ=255)
                                                                                                          cTGGCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCa                       >  1:671446/1‑149 (MQ=255)
                                                                                                           tGGCTATGTTACGCACACAGTCGAACACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAg                      <  2:282317/149‑1 (MQ=255)
                                                                                                            ggCTATTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAgg                     >  1:459646/1‑149 (MQ=255)
                                                                                                            ggCTACTTTACGGACACTGTCGCAGTGTGTGTCTCTGGCGGGGCTTCGTTCCGGATTCACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAgg                     <  2:663207/149‑1 (MQ=255)
                                                                                                                aTTTTACGCACACTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGTATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGtt                 <  2:77458/149‑1 (MQ=255)
                                                                                                                   ttACGCACCCTGTCGAAAACTTTGGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAg              <  1:214747/149‑1 (MQ=255)
                                                                                                                            cTGTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAGGCCCCCAggg     >  1:512252/1‑149 (MQ=255)
                                                                                                                              gTCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCTAACGAAGGAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACAGTGCGGACCAGGCGTTAAGCCCACAGGGaa   >  2:460800/1‑149 (MQ=255)
                                                                                                                              gTCGAAAACTCTTGCAATCGCGGGGCTGCGTTGCGGTCTTCCGCTGGCAAACGAAGAAGTCATCAACCCGCTGATGAAAGGGATCGCCCCCTACCCGCGCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGaa   <  2:660502/149‑1 (MQ=255)
                                                                                                                               tCGAAAACTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGCAt  >  1:396661/1‑149 (MQ=255)
                                                                                                                                     |                                                                                                                                              
CCTGCTGGAGTTAACCCGCGGTAAGGCGATTGTGGTTGCCGATGAAGCCTATATCGAGTTTTGCCCGCAGGCATCGCTGGCTGGCTGGCTGGCGGAATATCCGCACCTGGCTATTTTACGCACACTGTCGAAAGCTTTTGCTCTGGCGGGGCTTCGTTGCGGATTTACGCTGGCAAACGAAGAAGTCATCAACCTGCTGATGAAAGTGATCGCCCCCTACCCGCTCTCGACGCCGGTTGCCGACATTGCGGCCCAGGCGTTAAGCCCACAGGGAAT  >  NC_000913/2092907‑2093182

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