Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,090,860 G→A A223A (GCG→GCA hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,8600GA92.9% 124.1 / 0.5 42A223A (GCG→GCAhisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (21/18);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGC  >  NC_000913/2090771‑2090954
                                                                                         |                                                                                              
tgcctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCaa                                                                                      >  1:704673/1‑100 (MQ=255)
  cctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCaaaa                                                                                    <  1:1100461/100‑1 (MQ=255)
  cctgATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCaaaa                                                                                    <  1:518362/100‑1 (MQ=255)
      aTTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAACCAAAATAca                                                                                >  2:392115/1‑100 (MQ=255)
              cgATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCatgatg                                                                        >  2:1052303/1‑100 (MQ=255)
                aTGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGAGTGATCCAGGCACGCGAATCAAAATACATCATGATgca                                                                      >  2:1103481/1‑100 (MQ=255)
                      ctaTGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCacc                                                                <  2:908106/98‑1 (MQ=255)
                       aaTGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccg                                                               <  1:135914/100‑1 (MQ=255)
                        aTGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccag                                                              >  1:908114/1‑98 (MQ=255)
                         tGGAAGCATCCAAACACCAACTGATCGACAAACTGCTGACCCGTAATCACGGTGTGACACACACACCCGACCCAAAATACATCATGACGCACGCaccga                                                              >  2:997080/1‑99 (MQ=255)
                         tGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCaccgac                                                             >  2:513449/1‑100 (MQ=255)
                            aagaaTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCAccgaccga                                                          <  1:431301/100‑1 (MQ=255)
                              gaaTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAAc                                                        >  1:693485/1‑100 (MQ=255)
                                aTCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACg                                                       >  1:424371/1‑99 (MQ=255)
                                 tCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTc                                                     <  1:992259/100‑1 (MQ=255)
                                 tCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTc                                                     <  1:315772/100‑1 (MQ=255)
                                   cAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTg                                                   <  2:704669/100‑1 (MQ=255)
                                    aaaCAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTgg                                                  >  2:783130/1‑100 (MQ=255)
                                    aaaCAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATTAAAATACATCATGATGCACGCACCGACCGAACGt                                                      >  2:1038209/1‑96 (MQ=255)
                                          cAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAg                                            >  1:793346/1‑100 (MQ=255)
                                          cAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAg                                            <  1:138485/100‑1 (MQ=255)
                                              tGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCAt                                        >  1:708542/1‑100 (MQ=255)
                                                 tCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGc                                     <  1:706897/100‑1 (MQ=255)
                                                     cAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctg                                 >  2:281268/1‑100 (MQ=255)
                                                       aaCTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctgct                               >  1:896441/1‑100 (MQ=255)
                                                        aCTGCTGACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctgctg                              <  1:539999/100‑1 (MQ=255)
                                                           gctgACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCa                           >  1:438031/1‑100 (MQ=255)
                                                              gACCCGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGt                        >  1:947046/1‑100 (MQ=255)
                                                                cccGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGc                      >  1:384068/1‑100 (MQ=255)
                                                                 ccGTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGcc                     <  1:466951/100‑1 (MQ=255)
                                                                   gTATTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGa                   <  2:1193555/100‑1 (MQ=255)
                                                                     aTTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAAc                 >  2:161170/1‑100 (MQ=255)
                                                                     aTTCAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACAGAACGGCTGGATGGAGTCATCGCCCTGCTGCCAGGTGCCGAAc                 >  2:586030/1‑100 (MQ=255)
                                                                        cAGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGcc              <  1:970507/100‑1 (MQ=255)
                                                                         aGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGccc             >  2:497084/1‑100 (MQ=255)
                                                                         aGGGTGTGATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGccc             <  2:845222/100‑1 (MQ=255)
                                                                            gtgtGATCCAGGCACGCGAATCTCAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAAc          <  1:1103491/100‑1 (MQ=255)
                                                                               tgATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTAt       >  1:718788/1‑100 (MQ=255)
                                                                                gATCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTAtt      >  1:847815/1‑100 (MQ=255)
                                                                                 aTCCAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTc     <  1:189045/100‑1 (MQ=255)
                                                                                   ccAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTg   <  1:721011/100‑1 (MQ=255)
                                                                                    cAGGCACGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGc  >  1:62667/1‑100 (MQ=255)
                                                                                         |                                                                                              
TGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGC  >  NC_000913/2090771‑2090954

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