Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,982,860 C→A 7.1% intergenic (+49/‑98) proM → / → aslB tRNA‑Pro/putative AslA‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,982,8600CA 7.1% 191.8 / 3.1 56intergenic (+49/‑98)proM/aslBtRNA‑Pro/putative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (19/33);  new base A (2/2);  total (21/35)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGT  >  NC_000913/3982723‑3982992
                                                                                                                                         |                                                                                                                                    
accgCAGTCGTTCGGCGCGTTGCGCAGCGTGGGAGCGCATCTGGTTTGGGTCCAGTGGGTCGGAGGTGCGAATCGTCTCTCGCCGACCAATTTTGTACCCCGCGTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCAc                                                                                                                                    <  2:1227133/137‑1 (MQ=255)
   gCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCAccct                                                                                                                                 >  1:935985/1‑140 (MQ=255)
            cGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCa                                                                                                                        >  2:103283/1‑140 (MQ=255)
                     gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGcc                                                                                                               <  2:1117622/140‑1 (MQ=255)
                     gcgcAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACACCGCTTCGGCGGGGTTTTTTGTTTTATGTGCATTTCGTCACCCTCCCTTCGCAATAAACGcc                                                                                                               <  2:1158747/140‑1 (MQ=255)
                        cAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCGTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGAAATAAACGCCCGt                                                                                                            >  2:690828/1‑140 (MQ=255)
                         aGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGta                                                                                                           <  2:545857/140‑1 (MQ=255)
                           cTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGtaat                                                                                                         <  2:309830/140‑1 (MQ=255)
                           cTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGtaat                                                                                                         <  2:4866/140‑1 (MQ=255)
                                 aGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCa                                                                                                   <  1:1196439/140‑1 (MQ=255)
                                 aGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCa                                                                                                   <  1:593770/140‑1 (MQ=255)
                                 aGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCa                                                                                                   <  1:1195713/140‑1 (MQ=255)
                                     cAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGc                                                                                               >  1:490159/1‑140 (MQ=255)
                                     cAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGc                                                                                               >  1:110430/1‑140 (MQ=255)
                                       aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGccc                                                                                             <  1:1129318/140‑1 (MQ=255)
                                       aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGccc                                                                                             <  1:148999/140‑1 (MQ=255)
                                       aCTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGccc                                                                                             <  1:683608/140‑1 (MQ=255)
                                         tGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCa                                                                                           <  2:1121783/140‑1 (MQ=255)
                                         tGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCa                                                                                           <  2:1121781/140‑1 (MQ=255)
                                             ttGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGt                                                                                       >  1:282387/1‑140 (MQ=255)
                                                    cAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGAttt                                                                                <  1:103283/140‑1 (MQ=255)
                                                     aGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTc                                                                               >  1:576674/1‑140 (MQ=255)
                                                     aGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCAATGCCCCACGGTATGATTTc                                                                               >  2:734678/1‑140 (MQ=255)
                                                     aGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTAACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                                     >  2:16642/1‑118 (MQ=255)
                                                     aGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTAACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                                     >  2:15974/1‑118 (MQ=255)
                                                     aGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTAACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                                     <  1:16642/118‑1 (MQ=255)
                                                     aGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTAACCCTCCCTTCGCAATAAACGCCCGTAATAACt                                                                                                     <  1:15974/118‑1 (MQ=255)
                                                            cGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTa                                                                        >  1:336967/1‑140 (MQ=255)
                                                                    cgtaTCCTGTCTCGCCGTCCAATTTTGAACCCGGCTTGGGCGGGGTTTTTTGTTTTCTGTGCTTTTGGTCACCCTCCCTTGGCATTTTTCGCGCGTTATTAGTCATTGCCCGACGGTATGTTTTCGCCCTTTACGTTTTg                                                                <  2:277373/137‑1 (MQ=255)
                                                                    cGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTg                                                                >  1:977368/1‑140 (MQ=255)
                                                                    cGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTg                                                                <  2:277079/140‑1 (MQ=255)
                                                                    cGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTg                                                                <  2:277562/140‑1 (MQ=255)
                                                                         ccTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGa                                                           <  2:748746/140‑1 (MQ=255)
                                                                         ccTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGa                                                           <  2:276536/140‑1 (MQ=255)
                                                                          ctctctCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGAt                                                          <  2:841139/140‑1 (MQ=255)
                                                                          ctctctCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGAt                                                          <  1:719871/140‑1 (MQ=255)
                                                                           tctctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCAGGGTATGATTTCGCCCTTAACGTATTGAAGGATg                                                         >  2:786684/1‑140 (MQ=255)
                                                                           tctctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATg                                                         >  2:144194/1‑140 (MQ=255)
                                                                           tctctcGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATg                                                         >  2:454474/1‑140 (MQ=255)
                                                                               acGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGATTTCTGTGCATTACGTCACCCTCCCTTCCCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACtt                                                     <  1:254913/139‑1 (MQ=255)
                                                                                     aCCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCa                                               <  1:1006710/140‑1 (MQ=255)
                                                                                          ttttGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCATTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAg                                          <  1:417026/140‑1 (MQ=255)
                                                                                                     gCTTCGGCGGGGTATTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAtgctg                               <  2:631818/140‑1 (MQ=255)
                                                                                                      cTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCAtgctgc                              >  1:210999/1‑140 (MQ=255)
                                                                                                         cggcggGGATTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTACTAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAAc                           <  2:438878/140‑1 (MQ=255)
                                                                                                              gggTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGtt                      <  2:854632/140‑1 (MQ=255)
                                                                                                              gggTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGtt                      <  2:854644/140‑1 (MQ=255)
                                                                                                                      tGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGTCCATGCTGCAACAGGTTCCAAcgcg              <  2:922746/140‑1 (MQ=255)
                                                                                                                      tGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAcgcg              <  1:1029306/140‑1 (MQ=255)
                                                                                                                      tGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAAGATCTCATTGCCCCACGGTATGATTTCGCCCTTTCCGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAAcgcg              <  2:922758/140‑1 (MQ=255)
                                                                                                                           tCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCtt         >  2:509022/1‑140 (MQ=255)
                                                                                                                           tCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCtt         >  2:507445/1‑140 (MQ=255)
                                                                                                                           tCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCtt         >  1:949120/1‑140 (MQ=255)
                                                                                                                           tCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCtt         >  1:430597/1‑140 (MQ=255)
                                                                                                                              gtgCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTc      >  1:366783/1‑140 (MQ=255)
                                                                                                                                  aTTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCAtgt  <  1:189479/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
AAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGTTTGGGACCAGTGGGTCGGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGT  >  NC_000913/3982723‑3982992

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