Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,245,344 G→C 17.2% R64G (CGC→GGC)  pyrC ← dihydroorotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,245,3440GC17.2% 120.8 / 18.5 48R64G (CGC→GGC) pyrCdihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (21/18);  new base C (5/3);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGATGGATGT  >  NC_002947/1245195‑1245493
                                                                                                                                                     |                                                                                                                                                     
gCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCcgggcg                                                                                                                                                       >  2:894415/1‑150 (MQ=255)
      aTCGGAGTTGGTGGTGGCGCCGGCCGGGTATAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGc                                                                                                                                                 >  1:523855/1‑150 (MQ=255)
         ggAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAg                                                                                                                                              >  1:734842/1‑150 (MQ=255)
               ggtggtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATgcg                                                                                                                                        <  2:731243/150‑1 (MQ=255)
                   gtggCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCa                                                                                                                                    >  1:704235/1‑150 (MQ=255)
                              cggGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGcc                                                                                                                         >  2:179264/1‑150 (MQ=255)
                                        ttGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCgg                                                                                                                 >  2:956118/1‑148 (MQ=255)
                                        ttGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCgg                                                                                                                 <  1:956118/148‑1 (MQ=255)
                                              gcgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTg                                                                                                         >  1:299957/1‑150 (MQ=255)
                                               cgTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                                        <  2:704235/150‑1 (MQ=255)
                                                                gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                                            <  1:934539/129‑1 (MQ=255)
                                                                gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACa                                                                                                            >  2:934539/1‑129 (MQ=255)
                                                                gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCAGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                                        <  2:1009740/133‑1 (MQ=255)
                                                                gCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCAGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTgc                                                                                                        >  1:1009740/1‑133 (MQ=255)
                                                                     ggCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTggg                                                                                  <  1:179264/150‑1 (MQ=255)
                                                                               aCGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatgatgatg                                                                        <  1:36530/150‑1 (MQ=255)
                                                                                    cTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACg                                                                   >  2:687831/1‑150 (MQ=255)
                                                                                     tCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACgg                                                                  >  1:268408/1‑150 (MQ=255)
                                                                                     tCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACgg                                                                  >  1:637685/1‑150 (MQ=255)
                                                                                          gCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCAGCACGAGGTTGGGCATGATGATGGCACGGGCGaa                                                             <  2:179724/150‑1 (MQ=255)
                                                                                                 cggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                                      <  1:894415/150‑1 (MQ=255)
                                                                                                 cggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGCCAGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                                      <  2:938328/150‑1 (MQ=255)
                                                                                                  ggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                                      <  1:394069/149‑1 (MQ=255)
                                                                                                  ggtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTAcgc                                                      >  2:394069/1‑149 (MQ=255)
                                                                                                   gtcggtGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACgcgc                                                    >  1:152842/1‑150 (MQ=255)
                                                                                                         gAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTc                                              <  1:128155/150‑1 (MQ=255)
                                                                                                           ggTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatga                                                                             >  1:263824/1‑117 (MQ=255)
                                                                                                           ggTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatga                                                                             <  2:263824/117‑1 (MQ=255)
                                                                                                            gTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGAGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGcc                                           >  1:860776/1‑150 (MQ=255)
                                                                                                                aGCACCATCAGCGGCTCGAAGCAGGGGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACg                                       >  1:766204/1‑150 (MQ=255)
                                                                                                                  cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTa                                                         <  2:1060468/130‑1 (MQ=255)
                                                                                                                  cACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTa                                                         >  1:1060468/1‑130 (MQ=255)
                                                                                                                     cATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCATGGGCGGCACGAGGTTGGGCATCATGATCATATCGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                                  <  1:834672/150‑1 (MQ=255)
                                                                                                                     cATCAGCGGCTCGAAGCCGCGGGTGGCCGGGAGGGCCGCCAGGATGCGCTGACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGTGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGtgtg                                  <  2:196899/150‑1 (MQ=255)
                                                                                                                      aTCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTgg                                 >  2:447632/1‑150 (MQ=255)
                                                                                                                       tCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTggg                                <  2:152842/150‑1 (MQ=255)
                                                                                                                       tCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTggg                                >  1:277753/1‑150 (MQ=255)
                                                                                                                          gCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGaa                                                             <  2:30732/118‑1 (MQ=255)
                                                                                                                          gCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTAACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGaa                                                             >  1:30732/1‑118 (MQ=255)
                                                                                                                          gCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGATGCGGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCa                             >  1:812754/1‑150 (MQ=255)
                                                                                                                                 attctggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGc                      >  1:751775/6‑150 (MQ=255)
                                                                                                                                     tggCTGCCCGCCCGGCCGGCAGCCGGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatga                                                                             >  2:416682/2‑91 (MQ=38)
                                                                                                                                     tggCTGCCCGCCCGGCCGGCAGCCGGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCatga                                                                             <  1:416682/90‑1 (MQ=38)
                                                                                                                                      ggCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCAt                 >  1:1061644/1‑150 (MQ=255)
                                                                                                                                       gCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGc                                            >  1:1019752/1‑122 (MQ=255)
                                                                                                                                       gCTGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGc                                            <  2:1019752/122‑1 (MQ=255)
                                                                                                                                         tGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGACACCATcgc              >  1:1067738/1‑150 (MQ=255)
                                                                                                                                         tGCCCGCCCGGCCGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGACACCATcgc              <  2:1067738/150‑1 (MQ=255)
                                                                                                                                              gccggGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAga         <  2:637685/150‑1 (MQ=255)
                                                                                                                                                cgggcgggCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtgt                                   >  2:447196/1‑122 (MQ=255)
                                                                                                                                                cgggcgggCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGgtgt                                   <  1:447196/122‑1 (MQ=255)
                                                                                                                                                     gggcAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGATGGAtgt  >  1:428397/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGGAATCGGAGTTGGTGGTGGCGCCGGCCGGGTACAACTTGGCGGCGTACACGAAACCGCTGGCCTTGGCCGCGCGGACGTCCTCGGGGCTGGTGCGGTCGGTGAGGTACAGCACCATCAGCGGCTCGAAGCGGCTGCCGGCCGGGCGGGCAGCCAGAATGCGCTCACGGTAGGCGCCGGCCTCGGTGACATTGCGCACGGGCGGCACGAGGTTGGGCATGATGATGGCACGGGCGAAGGTACGCGCCACGTCGCCAACGGTGTGGGGCAGGACGGCACCATCGCGCAGATGGATGT  >  NC_002947/1245195‑1245493

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