Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,472,714 T→C intergenic (‑181/+100) pyrL ← / ← yjgH pyrBI operon leader peptide/RutC family protein YjgH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,7140TC98.1% 197.5 / ‑3.9 53intergenic (‑181/+100)pyrL/yjgHpyrBI operon leader peptide/RutC family protein YjgH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/30);  minor base A (0/1);  total (22/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 = 9.99e-01

GCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGCGCTATCACTT  >  NC_000913/4472581‑4472842
                                                                                                                                     |                                                                                                                                
gcgcAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGctttt                                                                                                                            <  2:531949/140‑1 (MQ=255)
   cAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGctttttcc                                                                                                                         <  2:497259/140‑1 (MQ=255)
   cAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGctttttcc                                                                                                                         >  1:381542/1‑140 (MQ=255)
     aGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGctttttcctt                                                                                                                       >  2:85390/1‑140 (MQ=255)
                 aaTATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTtttcactttcac                                                                                                           <  1:671380/140‑1 (MQ=255)
                  atatTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACg                                                                                                          >  1:759559/1‑140 (MQ=255)
                    atTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCa                                                                                                        <  1:113335/140‑1 (MQ=255)
                         gTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAg                                                                                                   <  1:1037865/140‑1 (MQ=255)
                         gTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAg                                                                                                   <  1:468341/140‑1 (MQ=255)
                           cAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGat                                                                                                 >  1:345536/1‑140 (MQ=255)
                             gcgcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGatat                                                                                               <  2:644593/140‑1 (MQ=255)
                             gcgcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGatat                                                                                               <  2:218556/140‑1 (MQ=255)
                             gcgcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGatat                                                                                               >  1:502114/1‑140 (MQ=255)
                              cgcAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATg                                                                                              >  1:881357/1‑140 (MQ=255)
                                               aaTTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCAt                                                                             <  1:51363/140‑1 (MQ=255)
                                               aaTTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCAt                                                                             <  1:216182/140‑1 (MQ=255)
                                               aaTTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCAGCTTTTTCCTTTTTCTCTTTCACGCAATCAGATATGCATTTTATTCCATTCAt                                                                             <  2:586258/140‑1 (MQ=255)
                                                 ttAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATcc                                                                           <  1:99592/140‑1 (MQ=255)
                                                  tAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCg                                                                          <  1:172652/140‑1 (MQ=255)
                                                        tCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTtat                                                                    <  2:898440/140‑1 (MQ=255)
                                                            ggTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAggg                                                                <  1:116875/140‑1 (MQ=255)
                                                              tataTGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCg                                                              <  1:733958/140‑1 (MQ=255)
                                                              tataTGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCg                                                              <  2:827900/140‑1 (MQ=255)
                                                              tataTGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCg                                                              <  1:1027681/140‑1 (MQ=255)
                                                              tataTGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCg                                                              <  1:30974/140‑1 (MQ=255)
                                                                      tttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGa                                                      >  1:651130/1‑140 (MQ=255)
                                                                       ttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGaa                                                     >  2:618147/1‑140 (MQ=255)
                                                                       ttttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGaa                                                     >  2:608053/1‑140 (MQ=255)
                                                                        tttATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTTTAGGGCGGAGTGTGaaa                                                    >  2:364937/1‑140 (MQ=255)
                                                                          tATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGTGGAGTGTGAAAGc                                                  <  1:1434/140‑1 (MQ=255)
                                                                          tATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGc                                                  >  1:881868/1‑140 (MQ=255)
                                                                           aTTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCg                                                 <  2:502114/140‑1 (MQ=255)
                                                                            tttGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGa                                                <  1:334659/140‑1 (MQ=255)
                                                                            tttGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGa                                                <  2:46729/140‑1 (MQ=255)
                                                                                aTGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACgg                                            >  1:360434/1‑140 (MQ=255)
                                                                                 tgtgATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCGCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGc                                           >  1:238621/1‑140 (MQ=255)
                                                                                  gtgATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCt                                          <  2:273266/140‑1 (MQ=255)
                                                                                    gATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTaa                                        >  1:67765/1‑140 (MQ=255)
                                                                                    gATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTaa                                        >  1:113296/1‑140 (MQ=255)
                                                                                      ttCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAaca                                      <  1:85390/140‑1 (MQ=255)
                                                                                                gAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                            <  2:381542/140‑1 (MQ=255)
                                                                                                gAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                            <  2:113296/140‑1 (MQ=255)
                                                                                                gAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                            <  2:881868/140‑1 (MQ=255)
                                                                                                gAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                            <  2:601396/140‑1 (MQ=255)
                                                                                                gAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTa                            <  2:67765/140‑1 (MQ=255)
                                                                                                   cTGTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTg                         >  2:374005/1‑140 (MQ=255)
                                                                                                     gTCCTTCGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCt                       >  1:1038658/1‑140 (MQ=255)
                                                                                                           cGGTATCAATAATTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGa                 >  2:43334/1‑140 (MQ=255)
                                                                                                                  aataatTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATgcgcgc          >  2:558077/1‑140 (MQ=255)
                                                                                                                     aatTGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGCAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGCGCTAt       <  2:759559/140‑1 (MQ=255)
                                                                                                                       ttGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGCGCTATCa     >  2:861061/1‑140 (MQ=255)
                                                                                                                        tGCAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGCGCTATCAc    >  1:724908/1‑140 (MQ=255)
                                                                                                                          cAGACATTCCCCGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGCGCTATCACtt  >  1:782417/1‑140 (MQ=255)
                                                                                                                                     |                                                                                                                                
GCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGCATTTTATTCCATTCATCCGACTTATAGGGCGGAGTGTGAAAGCGAACGGCTAACACTATTGCTTACTGCTCAGGGATGCGCGCTATCACTT  >  NC_000913/4472581‑4472842

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