Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,703,303 G→C 6.9% E710Q (GAA→CAA)  rsxC → SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,703,3030GC 6.9% 913.9 / 17.1 299E710Q (GAA→CAA) rsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base G (147/128);  new base C (20/4);  total (167/132)
Fisher's exact test for biased strand distribution p-value = 4.80e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 330 total aligned reads displayed.

CGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCG  >  CP009273/1703236‑1703369
                                                                   |                                                                  
 gAAGCGGCTATTGCACGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCagaaga                                                                   <  1:30850911/68‑1 (MQ=255)
   aGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAAc                                                                 >  1:20344090/1‑68 (MQ=255)
      tgCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:3947995/66‑1 (MQ=255)
      ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:28412228/67‑1 (MQ=255)
      ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:10814126/67‑1 (MQ=255)
      ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:15953765/67‑1 (MQ=255)
      ggCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:7019917/67‑1 (MQ=255)
      cgCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:13398690/66‑1 (MQ=255)
      cgCTATCGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:10348462/66‑1 (MQ=255)
      agCTATTGCACGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAg                                                               <  1:5663546/66‑1 (MQ=255)
              cACGCGCCAACGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATcc                                                       >  1:16743625/1‑67 (MQ=255)
              cACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATcc                                                       >  1:3492587/1‑67 (MQ=255)
              cACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATcc                                                       >  1:28642350/1‑67 (MQ=255)
              cACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATcc                                                       >  1:9458355/1‑67 (MQ=255)
                    ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCaaa                                                >  1:30179776/1‑68 (MQ=255)
                     cAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCaaa                                                <  1:13072556/67‑1 (MQ=25)
                     cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCaaa                                                <  1:5162634/67‑1 (MQ=255)
                      aaaGCTCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGc                                              <  1:29641179/68‑1 (MQ=255)
                        aGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCg                                             >  1:12739218/1‑67 (MQ=255)
                         gCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGc                                           >  1:12815183/1‑68 (MQ=255)
                         gCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGc                                           >  1:15111685/1‑68 (MQ=255)
                         gCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGc                                           >  1:11808843/1‑68 (MQ=255)
                           aCGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAg                                         >  1:19579317/1‑68 (MQ=255)
                              cAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGta                                       >  1:6659959/1‑66 (MQ=11)
                                aGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGcc                                    >  1:27906452/1‑68 (MQ=255)
                                 gCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGcc                                    >  1:27113583/1‑67 (MQ=255)
                                   tGGACCAGCAACAGGCTAATGAGGTACAAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgcg                                 <  1:23880398/68‑1 (MQ=18)
                                   tGGAACAGCAACAGGCTATTGCGGTAACAGATCAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgcg                                 <  1:7633330/68‑1 (MQ=18)
                                   tGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgcg                                 <  1:10550408/68‑1 (MQ=255)
                                   tGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgcg                                 <  1:25469785/68‑1 (MQ=255)
                                   gggAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgcg                                 <  1:21983497/67‑1 (MQ=255)
                                    ctaaCAGAAACGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:8302974/66‑1 (MQ=18)
                                     gAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCcgcg                                 >  1:27563832/1‑66 (MQ=255)
                                     caaCAGCAACCGGCTAATGCGGAACCAGAACAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:310558/67‑1 (MQ=17)
                                      aaCCGCAACCGGCTAAAGCGCAACCAGAAGAACAGGTCGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:336175/67‑1 (MQ=14)
                                      aaCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:1916571/67‑1 (MQ=255)
                                      aaCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:11184527/67‑1 (MQ=255)
                                      aaCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:22351950/67‑1 (MQ=255)
                                      aaCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:12174132/67‑1 (MQ=255)
                                      aaCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGc                               <  1:31734829/67‑1 (MQ=255)
                                       aCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTa                             >  1:18149909/1‑68 (MQ=255)
                                       aCAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTa                             >  1:4268690/1‑68 (MQ=255)
                                       aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGt                                        >  1:32089916/1‑57 (MQ=11)
                                       aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTTGCCGCGGCTa                             >  1:8402679/1‑68 (MQ=17)
                                       aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTAGCAGCGGCTa                             >  1:33919442/1‑68 (MQ=14)
                                       aCAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTAGACGCGGCTa                             >  1:20497635/1‑68 (MQ=14)
                                       aCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTTGATCCGCGCAAAGCGGCAGTAGACGCGGCTa                             >  1:1125831/1‑68 (MQ=14)
                                       aCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCGGCAGTTGCCGCGGCTa                             >  1:26115487/1‑68 (MQ=14)
                                       aCAGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCCGCAGTTGCCGCGGCTa                             >  1:33352974/1‑68 (MQ=14)
                                       aCAGCAACAGGCTAATGCGGAACCAAAACAACAGGTAGATCCGCGGAAAGCGGCAGTAGCCGCCGCTa                             >  1:14568449/1‑68 (MQ=11)
                                         aGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAtt                           >  1:13674821/1‑68 (MQ=255)
                                         aGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCGGCCGTAGCCGCGGCTAtt                           >  1:5983627/1‑68 (MQ=11)
                                         aGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAtt                           >  1:7685258/1‑68 (MQ=17)
                                         aGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCGCGCAAAGCGGCAGTAGCCGCGGCTAtt                           >  1:22851576/1‑68 (MQ=14)
                                            aaCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGcc                        <  1:17779694/68‑1 (MQ=255)
                                            aaCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGcc                        <  1:6468419/68‑1 (MQ=255)
                                            aaCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGcc                        <  1:24763417/68‑1 (MQ=255)
                                            aaCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGcc                        <  1:26349739/68‑1 (MQ=255)
                                            aaCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGcc                        <  1:10064314/68‑1 (MQ=255)
                                            aaCAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGc                         <  1:32874533/67‑1 (MQ=255)
                                            aaCAGGCTAATGCGGTACCAGAAGAACACGATGATCAGCGCAAAGCGGCAGTTGCCGAGGCTATAGc                         >  1:1377014/1‑67 (MQ=21)
                                            aaCAGGCTAATGCGGAACCAGAAGAACAGGTAGATCCGCGCAAAGCGGCAGTAGCCGCGGCTATTGc                         >  1:4466427/1‑67 (MQ=17)
                                                 gCTAATGCGGTACCAGAAGAACAGGTTGTTCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgc                     >  1:16074075/1‑66 (MQ=255)
                                                 gCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgc                     >  1:30535580/1‑66 (MQ=255)
                                                 gCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgc                     >  1:14029382/1‑66 (MQ=255)
                                                 gCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGTGGCCGCGGCTATTGCCcgcg                    >  1:789541/1‑67 (MQ=11)
                                                  cTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCgcgc                   >  1:9331566/1‑67 (MQ=255)
                                                  cTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCgcgc                   <  1:15400727/67‑1 (MQ=255)
                                                  cTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCgcgc                   <  1:31449003/67‑1 (MQ=255)
                                                  cTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCgcgc                   >  1:11592953/1‑67 (MQ=255)
                                                  cTAATGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCcgcg                    <  1:21650876/66‑1 (MQ=255)
                                                  atAATGCGGTACAAGCAGAACAGGTTGATCCGCGCAAAGCTGCAGTTGCCGCGGCTATTGCCCgcgc                   <  1:33758710/66‑1 (MQ=18)
                                                   tAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCt                  <  1:14236570/67‑1 (MQ=255)
                                                    aaTGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCa                <  1:32371860/68‑1 (MQ=255)
                                                      tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:21132202/68‑1 (MQ=255)
                                                      tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:30022111/68‑1 (MQ=255)
                                                      tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:17963699/68‑1 (MQ=255)
                                                      tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:28163534/68‑1 (MQ=255)
                                                      tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:17071110/68‑1 (MQ=255)
                                                      tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:26685635/68‑1 (MQ=255)
                                                      tGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:26220060/68‑1 (MQ=255)
                                                      tGCGGAACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              <  1:7174644/68‑1 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              >  1:25867050/1‑67 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              >  1:2486841/1‑67 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              >  1:24495761/1‑67 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              >  1:23006725/1‑67 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              >  1:29148974/1‑67 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              >  1:29349966/1‑67 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGTTGATCCGAGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAgg              >  1:25386686/1‑67 (MQ=255)
                                                       gCGGTACCAGAAGAACAGGATGATCCGCGCAAAGCGGCAGTTGCCGAGGCAATAGCCCTCGCTCacg              >  1:9619465/1‑65 (MQ=25)
                                                            aCCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaa        >  1:15691439/1‑68 (MQ=255)
                                                              cAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCaaaaaa      >  1:18761157/1‑68 (MQ=255)
                                                                 aagaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc   >  1:19896971/1‑68 (MQ=255)
                                                                 aagaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc   >  1:1927084/1‑68 (MQ=255)
                                                                 aagaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc   >  1:4776135/1‑68 (MQ=255)
                                                                 aagaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc   >  1:452416/1‑68 (MQ=255)
                                                                 aagaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTAATGCCCGCGCACAGGCCAAAAAAgcc   >  1:26827447/1‑68 (MQ=255)
                                                                 aagaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCACGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc   >  1:10995650/1‑68 (MQ=255)
                                                                  agaaCAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAgcc   <  1:3094178/67‑1 (MQ=255)
                                                                  agaaCAGGTCGATCCGCGCAAAGCGGCAGTAGCCGCGGCTATTGCCCGCGCGAAGGCCAAAAAAgcc   >  1:14748851/1‑67 (MQ=25)
                                                                   |                                                                  
CGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGTACCAGAAGAACAGGTTGATCCGCGCAAAGCGGCAGTTGCCGCGGCTATTGCCCGCGCTCAGGCCAAAAAAGCCG  >  CP009273/1703236‑1703369

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