Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,282,421 A→C 6.7% intergenic (+439/+101) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,282,4210AC 6.7% 805.2 / 3.0 397intergenic (+439/+101)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (141/208);  new base C (28/6);  total (177/220)
Fisher's exact test for biased strand distribution p-value = 2.47e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 405 total aligned reads displayed.

GAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGA  >  CP009273/1282355‑1282488
                                                                  |                                                                   
  aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACg                                                                   <  1:12961628/67‑1 (MQ=35)
     aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCg                                                               >  1:22673304/1‑68 (MQ=35)
       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                             <  1:16417294/68‑1 (MQ=35)
       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                             <  1:26787807/68‑1 (MQ=35)
       tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                             <  1:2373204/68‑1 (MQ=35)
        aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                             >  1:27739820/1‑67 (MQ=35)
          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGGTGACCCTCGAGCCg                                                          >  1:5214937/1‑68 (MQ=255)
          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                          >  1:29693379/1‑68 (MQ=40)
          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                          >  1:2982609/1‑68 (MQ=40)
          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCg                                                          >  1:9464166/1‑68 (MQ=40)
          tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCGCGAGGCg                                                          >  1:3326832/1‑68 (MQ=255)
           cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                         <  1:28310560/68‑1 (MQ=40)
            aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                         >  1:13531886/1‑67 (MQ=40)
             gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGaa                                                        >  1:21839112/1‑67 (MQ=25)
             gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                        >  1:7597872/1‑67 (MQ=40)
             gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGaa                                                        >  1:4656959/1‑67 (MQ=25)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTTCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCt                                                     >  1:4367583/1‑68 (MQ=255)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAAccg                                                     >  1:6592830/1‑67 (MQ=255)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                        <  1:20705009/65‑1 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                        <  1:1423756/65‑1 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                      <  1:22041170/67‑1 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                       <  1:616291/66‑1 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                       <  1:15680441/66‑1 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAc                                                       <  1:7375243/66‑1 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                     >  1:26285092/1‑68 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                     >  1:16891073/1‑68 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                     >  1:18836472/1‑68 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                     >  1:21495274/1‑68 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                     >  1:10547598/1‑68 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                     >  1:13283898/1‑68 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCt                                                     >  1:1318765/1‑68 (MQ=40)
               gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCGCGAGCCGAAccg                                                     >  1:20182732/1‑67 (MQ=255)
                 tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                   <  1:6449760/68‑1 (MQ=255)
                  gcgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAg                                                  <  1:4058016/68‑1 (MQ=255)
                    gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGCGCCGAACCTTAGTc                                                >  1:25257026/1‑68 (MQ=255)
                    gcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTc                                                >  1:17110065/1‑68 (MQ=255)
                         tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGaa                                             >  1:17675310/1‑66 (MQ=255)
                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                          <  1:8990472/68‑1 (MQ=255)
                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                          <  1:4959817/68‑1 (MQ=255)
                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                          <  1:15409832/68‑1 (MQ=255)
                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                          <  1:8139678/68‑1 (MQ=255)
                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                          <  1:19514515/68‑1 (MQ=255)
                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                          <  1:19726292/68‑1 (MQ=255)
                          cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                          <  1:136830/68‑1 (MQ=255)
                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                     >  1:9571956/1‑67 (MQ=255)
                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCa                                     >  1:17527187/1‑67 (MQ=255)
                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                    <  1:16667428/68‑1 (MQ=255)
                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                    <  1:27372627/68‑1 (MQ=255)
                                tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAt                                    <  1:9251027/68‑1 (MQ=255)
                                   ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCt                                 >  1:14766225/1‑68 (MQ=255)
                                   ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCt                                 >  1:6824825/1‑68 (MQ=255)
                                     tCGTTGCCTGCGTCTACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                <  1:12364767/67‑1 (MQ=255)
                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGCGCCGCCCCTTAGTCGAAGCTTCTCAGCCtt                                >  1:12501616/1‑67 (MQ=14)
                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                <  1:3025461/67‑1 (MQ=255)
                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                <  1:23238459/67‑1 (MQ=255)
                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                <  1:10074550/67‑1 (MQ=255)
                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                <  1:22510764/67‑1 (MQ=255)
                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCtt                                <  1:16195027/67‑1 (MQ=255)
                                     tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCtt                                >  1:23915413/1‑67 (MQ=255)
                                      cGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                               <  1:28108687/67‑1 (MQ=255)
                                         tGCCTGCGGCAACGCTCTCTCGCTTACGCTCTAGGCGAACCTTAGTCGAAGCTTCTCATCCTTcccct                           <  1:7989980/68‑2 (MQ=14)
                                          gCCTGCGGCAACGCTCTCTCGCTGACGCTCGGGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                          >  1:18188628/1‑68 (MQ=255)
                                          gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGCACCTTAGTCGAAGCTTCTCATCCTTCCCCGc                          >  1:5919137/1‑68 (MQ=255)
                                          gCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAGCCTTCCCCga                          >  1:24986680/1‑67 (MQ=255)
                                           ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                         >  1:5081653/1‑68 (MQ=255)
                                              gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                      <  1:15134459/68‑1 (MQ=255)
                                              gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCCTAGTCGAACCTTCTCATCCTTCCCCGCATgg                      <  1:20101168/68‑1 (MQ=255)
                                                ggCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                     <  1:24603001/67‑1 (MQ=255)
                                                ggCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                     <  1:10241787/67‑1 (MQ=255)
                                                ggCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                     <  1:3826962/67‑1 (MQ=255)
                                                 gCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGAAggggga                   >  1:28907363/1‑66 (MQ=14)
                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                 <  1:25699142/68‑1 (MQ=255)
                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                 <  1:17936207/68‑1 (MQ=255)
                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                 <  1:11691384/68‑1 (MQ=255)
                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                 <  1:28739963/68‑1 (MQ=255)
                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                 <  1:5679171/68‑1 (MQ=255)
                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGa                 <  1:24158599/68‑1 (MQ=255)
                                                   aaCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGa                 >  1:7716368/1‑68 (MQ=255)
                                                     cGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:12775982/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:5522399/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:25817399/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:8347015/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:6989524/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:8654641/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:10668034/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:10608993/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:27018375/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:21446585/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:17313518/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:29225739/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:18568160/68‑1 (MQ=255)
                                                     cGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAat               <  1:20927860/68‑1 (MQ=255)
                                                      gCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAata              >  1:20573410/1‑68 (MQ=255)
                                                        tctctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat             <  1:14016070/67‑1 (MQ=255)
                                                           ctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtttta         >  1:14627230/1‑66 (MQ=255)
                                                           ctcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAtttaa         >  1:11470426/1‑66 (MQ=255)
                                                            tcGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt        <  1:11286456/68‑1 (MQ=255)
                                                            tcGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAt        <  1:2928863/68‑1 (MQ=255)
                                                              gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAtt       >  1:11912557/1‑67 (MQ=255)
                                                                  aCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCCTCCTTCCGAGAATGGGCTCAATAAATGATTGCggt  >  1:10879884/1‑67 (MQ=12)
                                                                  |                                                                   
GAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGA  >  CP009273/1282355‑1282488

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