Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,708,586 C→T 8.6% P222P (CCG→CCA cysI ← sulfite reductase hemoprotein subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,708,5860CT 8.6% 160.4 / 5.4 59P222P (CCG→CCAcysIsulfite reductase hemoprotein subunit beta
Reads supporting (aligned to +/- strand):  ref base C (27/26);  new base T (3/2);  total (30/29)
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.22e-01

ATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGCG  >  NC_002947/2708440‑2708734
                                                                                                                                                   |                                                                                                                                                     
aTTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGc                                                                                                                                                    >  2:214694/1‑150 (MQ=255)
 tttCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCa                                                                                                                                                   <  1:214694/150‑1 (MQ=255)
   tCGGCCTCGACCTTTTCAGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAgg                                                                                                                                                 <  1:123430/150‑1 (MQ=255)
     ggCCTCGACCTTTTCGGCGAACACTTCGAGTTTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGac                                                                                                                                               <  1:240/150‑2 (MQ=255)
          cGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGt                                                                                                                                          >  2:1424218/1‑150 (MQ=255)
                    ggCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGa                                                                                                                               <  1:931490/150‑1 (MQ=255)
                       gtacacTTCCGGGGGAAGGGCCTTGACCAGGATCTTGATCCGCTCACTGTAGTTTTTGTCACGGCGACCGTACCGGTTGTACACACGCAGGATGTCGTCGAGGTAGCTGATTAGGTCCTTCTACGGCAGGAACTCGTTGATGAActaacc                                                                                                                            >  2:1388436/3‑145 (MQ=255)
                          cacTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCccac                                                                                                                         <  2:27447/150‑1 (MQ=255)
                             ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCcaccac                                                                                                                      <  1:1068991/150‑1 (MQ=255)
                             ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACACAACACGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCcaccac                                                                                                                      >  1:1206415/1‑150 (MQ=255)
                               ccGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCg                                                                                                                    <  2:998507/150‑1 (MQ=255)
                                      aaGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACg                                                                                                             <  2:395281/150‑1 (MQ=255)
                                        gggCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGc                                                                                                           >  2:117457/1‑150 (MQ=255)
                                        gggCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGc                                                                                                           >  1:533504/1‑150 (MQ=255)
                                        gggCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACCCGCAGGATGGCGGCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCTTGAAGCAACGTGCTCACCACCGGAGGACGGc                                                                                                           >  2:616976/1‑150 (MQ=255)
                                           ccTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCa                                                                                                        >  1:326107/1‑150 (MQ=255)
                                            cTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAg                                                                                                       <  1:364651/150‑1 (MQ=255)
                                                   aGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacca                                                                                                <  2:449596/150‑1 (MQ=255)
                                                        cTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGCCAGCAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTAc                                                                                           <  1:1170757/150‑1 (MQ=255)
                                                         ttGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTAcc                                                                                          >  2:411303/1‑150 (MQ=255)
                                                                       gTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGcc                                                                            >  1:1382489/1‑150 (MQ=255)
                                                                          cTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAg                                                                         >  1:1377499/1‑150 (MQ=255)
                                                                           ttgttgTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCTCCGGAGTACTGCCCAGGCCACCACCTACGAGCACACGGAAGCCCAGt                                                                        <  2:521307/150‑1 (MQ=255)
                                                                           ttgttgTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGt                                                                        >  2:1151315/1‑150 (MQ=255)
                                                                                 tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGcc                                                                                                    >  1:1074700/1‑116 (MQ=255)
                                                                                 tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGcc                                                                                                    <  2:1074700/116‑1 (MQ=255)
                                                                                 tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGcc                                                                   <  2:412184/149‑1 (MQ=255)
                                                                                 tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGcc                                                                   >  1:412184/1‑149 (MQ=255)
                                                                                 tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCg                                                                  <  1:1083387/150‑1 (MQ=255)
                                                                                 tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCg                                                                  >  1:1215853/1‑150 (MQ=255)
                                                                                         aCCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGtt                                                          <  2:638079/150‑1 (MQ=255)
                                                                                          ccGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCa                                                                          <  2:1005117/133‑1 (MQ=255)
                                                                                          ccGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCa                                                                          >  1:1005117/1‑133 (MQ=255)
                                                                                              tAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTgcgc                                                      >  1:816996/1‑150 (MQ=255)
                                                                                               aaCGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCa                                                     <  2:1215853/150‑1 (MQ=255)
                                                                                                aCGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCAc                                                    <  1:507636/150‑1 (MQ=255)
                                                                                                  ggTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCg                                                  >  1:680955/1‑150 (MQ=255)
                                                                                                    ttGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGac                                                                                                                                               >  1:1089524/1‑55 (MQ=255)
                                                                                                    ttGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGac                                                                                                                                               <  2:1089524/56‑2 (MQ=255)
                                                                                                              gcAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCaa                                      >  2:803839/1‑150 (MQ=255)
                                                                                                               cAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAAt                                     <  1:411303/150‑1 (MQ=255)
                                                                                                                        gTCGAGGTATCTGATCAGGTCCTGCCCAGGCCGAAACTCGTTGATGAACGAGCCCACCACCGGGGTAAGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCAc                            <  2:1366763/150‑1 (MQ=255)
                                                                                                                             ggTAGCTGATCATGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCAAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTCGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGa                       <  1:616976/150‑1 (MQ=255)
                                                                                                                             ggTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGa                       <  1:1424218/150‑1 (MQ=255)
                                                                                                                                  cTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATggcg                  <  2:326107/150‑1 (MQ=255)
                                                                                                                                   tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCGCACCACCGGGGTACGTCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGGTggcgg                 >  1:522669/1‑150 (MQ=255)
                                                                                                                                   tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGCAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGc                                                >  2:629757/1‑119 (MQ=255)
                                                                                                                                   tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGc                                                <  1:629757/119‑1 (MQ=255)
                                                                                                                                      cgagttCCTGCCGTGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGgcgc              >  1:994434/6‑150 (MQ=255)
                                                                                                                                       gagttCCTGCCGTGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGcgcg             >  2:285055/5‑150 (MQ=255)
                                                                                                                                        aGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAg                                                                              >  1:859575/1‑84 (MQ=255)
                                                                                                                                        aGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAg                                                                              <  2:859575/84‑1 (MQ=255)
                                                                                                                                          gTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTNCGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGc                                                >  1:873778/1‑112 (MQ=255)
                                                                                                                                          gTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGc                                                <  2:873778/112‑1 (MQ=255)
                                                                                                                                          gTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTc          >  1:204400/1‑150 (MQ=255)
                                                                                                                                             cTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTc       >  2:228257/1‑150 (MQ=255)
                                                                                                                                                ccACGGCAGGAACTCGTTGATGAACGAGCCCACCGCCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTGCGATGGCGGCGCGGTCTTCCTg    <  2:616792/150‑1 (MQ=255)
                                                                                                                                                ccACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGACCCACCTTGAGCGACCAGCGGTGCGCACCGGCTCAAGCACCATGACGACCACCACGGGGGCGGCGGGGTCTTCCTg    >  1:26031/1‑150 (MQ=255)
                                                                                                                                                  aCGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGCg  >  2:1367672/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGCG  >  NC_002947/2708440‑2708734

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