Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,271,941 G→A 7.8% D48D (GAC→GAT) ‡ PP_2869 ← FMN‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,271,9410GA 7.8% 148.0 / 2.3 52D48D (GAC→GAT) ‡PP_2869FMN‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (22/25);  new base A (2/2);  total (25/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

CCGTGCGCATGGCAGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAGTGGGAAGGCGGTCTGCATCAGGC  >  NC_002947/3271795‑3272089
                                                                                                                                                   |                                                                                                                                                    
ccGTGCGCATGGCAGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTc                                                                                                                                                    <  7:288690/149‑1 (MQ=255)
      gcATGGCAGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTg                                                                                                                                              <  1:181455/149‑1 (MQ=255)
           gCAGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAgg                                                                                                                                         >  4:38001/1‑149 (MQ=255)
           gCAGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGAGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGCTCAGGCCGACGCCGCCGGCCGCGCGCTCGCGCTGGGAgg                                                                                                                                         >  8:74505/1‑149 (MQ=255)
             aGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGTACGATCAGGCCGACGCCGCCTACCAGCGCGACCGCGCTgccggcg                                                                                                                                        <  7:174004/149‑7 (MQ=255)
             aGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTAAGTGGCCATGAGCACATGCGAGGTGTAGCGGACACAGTCGTGCACCCCGCCGACCTGCAGCACGATCAGGCCGAAGCCGCCGGCAGCCCGGTCCCCCTGGTAGGc                                                                                                                                        >  7:111139/1‑148 (MQ=255)
              ggCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGa                                                                                                                                      <  1:181651/149‑1 (MQ=255)
                 cTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCa                                                                                                                                   >  8:268029/1‑149 (MQ=255)
                    ggCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAg                                                                                                                                  <  8:61794/147‑1 (MQ=255)
                    ggCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAg                                                                                                                                  <  8:204234/147‑1 (MQ=255)
                        aGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTc                                                                                                                            <  3:38001/149‑1 (MQ=255)
                        aGCTGGCGGTAGCCGTCGATGCAGGCATCGGCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTc                                                                                                                            <  7:74505/149‑1 (MQ=255)
                             gCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGtt                                                                                                                         <  4:27207/147‑1 (MQ=255)
                             gCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGa                                                                                                                       <  8:109275/149‑1 (MQ=255)
                                       tCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTc                                                                                                             <  2:75045/149‑1 (MQ=255)
                                            gCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTgg                                                                                                        >  6:228544/1‑149 (MQ=255)
                                               ggCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGTCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCtt                                                                                                                               <  8:47356/123‑1 (MQ=255)
                                               ggCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCtt                                                                                                                               >  7:47356/1‑123 (MQ=255)
                                                     gTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGc                                                                                                                                 >  3:123945/1‑115 (MQ=255)
                                                     gTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGc                                                                                                                                 <  4:123945/115‑1 (MQ=255)
                                                       cAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTg                                                                                                         <  2:158278/137‑1 (MQ=255)
                                                       cAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTg                                                                                                         <  2:158276/137‑1 (MQ=255)
                                                       cAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTg                                                                                                         >  1:158276/1‑137 (MQ=255)
                                                       cAGAGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTg                                                                                                         >  1:158278/1‑137 (MQ=255)
                                                        aGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGgtgt                                                                                            >  2:112197/1‑149 (MQ=255)
                                                           ggCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGt                                                                                         <  5:228544/149‑1 (MQ=255)
                                                               aTTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGAc                                                                                      <  3:118620/148‑1 (MQ=255)
                                                               aTTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGAc                                                                                      >  4:118620/1‑148 (MQ=255)
                                                               aTTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGACCGCGCTgccg                                                                                                                                           >  4:112008/1‑92 (MQ=255)
                                                               aTTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGACCGCGCTgccg                                                                                                                                           <  3:112008/95‑4 (MQ=255)
                                                                    cacaTGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTgg                                                                                <  7:268029/149‑1 (MQ=255)
                                                                     acaTGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGa                                                                               >  1:148998/1‑149 (MQ=255)
                                                                     acaTGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGa                                                                               <  5:121610/149‑1 (MQ=255)
                                                                       aTGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACa                                                                             >  8:265506/1‑149 (MQ=255)
                                                                                 tAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGt                                                                    <  1:112197/148‑1 (MQ=255)
                                                                                             tCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGttt                                                       <  6:199902/149‑1 (MQ=255)
                                                                                                 gCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCcg                                                   >  6:95902/1‑149 (MQ=255)
                                                                                                  cACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCcg                                                   <  5:95902/148‑1 (MQ=255)
                                                                                                   aCCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGGTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGGTTTTCAAGCGTTTGCcgcg                                                 >  5:204039/1‑149 (MQ=255)
                                                                                                     cccGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTTTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGttttt                                                                <  2:244863/132‑1 (MQ=255)
                                                                                                     cccGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGttttt                                                                >  1:244863/1‑132 (MQ=255)
                                                                                                                  agcaCGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTc                                  <  2:148998/149‑1 (MQ=255)
                                                                                                                       gATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGaa                               >  4:219183/1‑147 (MQ=255)
                                                                                                                        aTCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAg                            <  5:44819/149‑1 (MQ=255)
                                                                                                                           aGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAg                         >  8:221147/1‑149 (MQ=255)
                                                                                                                             gCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAGTg                       >  8:182717/1‑149 (MQ=255)
                                                                                                                               cGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAGTgg                      >  4:285277/1‑148 (MQ=255)
                                                                                                                                   gccgccGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGc                                                           >  2:247138/1‑107 (MQ=255)
                                                                                                                                   gccgccGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGc                                                           <  1:247138/107‑1 (MQ=255)
                                                                                                                                   gccgccGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAGTGGGAAgg                 >  4:9133/1‑149 (MQ=255)
                                                                                                                                     cgccgGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAGTGGGAAGGCg               >  2:112973/1‑149 (MQ=255)
                                                                                                                                       cctaccAGCGCGACCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCGGGAAGGGATCGAACAGGAGTCTGATGGCgg              >  8:296640/6‑149 (MQ=255)
                                                                                                                                                   gTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAGTGGGAAGGCGGTCTGCATCAggc  <  7:265506/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCGTGCGCATGGCAGGCCTCGGCAAGCTGGCGGTAGCCGTCGATGCAGGCATCGTCAGTGGCCATTAGCACATGCGAGGTGTAGCGGGCACTGTCGTGCACCCCGGCGACCTGCAGCACGATCAGGCCGACGCCGCCGGCAGCGCGGTCGCGCTGGTAGGCGATCAGCTTGTCGTTGACCAGGTTGTCGGTGGGCAGGCAGGTGTCGTGACCGGTGGACATGATGCGGTTTTTCAAGCGTTTGCCGCGGATCTGCAAGGGCTCGAACAGGAGTGGGAAGGCGGTCTGCATCAGGC  >  NC_002947/3271795‑3272089

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