Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,907,682 G→A 6.1% intergenic (‑91/+579) yobF ← / ← yebO DUF2527 family heat‑induced protein/putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,907,6820GA 6.1% 223.8 / 4.1 65intergenic (‑91/+579)yobF/yebODUF2527 family heat‑induced protein/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (25/36);  new base A (2/2);  total (27/38)
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.68e-01

ACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTT  >  NC_000913/1907545‑1907818
                                                                                                                                         |                                                                                                                                        
acgtcaacgaGAAAACACACGACTTCTTTACTGAAAATGCCACACATCAACAGAACTGTACCTCGTGTAACCCGAAATCTCTTATCACACACACCGTTAAATATGGCAAGCCATTTTTAAACACGTCTCGATCAGACGca                                                                                                                                        <  2:900224/140‑1 (MQ=255)
 cgtcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGcac                                                                                                                                       >  1:346411/1‑140 (MQ=255)
   tcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACaa                                                                                                                                     <  2:256110/140‑1 (MQ=255)
   tcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACaa                                                                                                                                     <  1:329350/140‑1 (MQ=255)
   tcaacgtGTTTACTCAGGACTTCTTTACTGAAAATGCAACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACaa                                                                                                                                     <  2:811996/140‑1 (MQ=255)
     aacgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAt                                                                                                                                   <  1:374921/140‑1 (MQ=255)
      acgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtc                                                                                                                                  >  1:517849/1‑140 (MQ=255)
      acgtGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAAtc                                                                                                                                  >  1:517883/1‑140 (MQ=255)
           attACTCAGGACTTCTTTACTGAAAATGCCACACAAAAACAGAACTGTACCTCGTTTAACCCGAAATCTCTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGAt                                                                                                                             <  1:150410/139‑1 (MQ=255)
            ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa                                                                                                                            >  2:588819/1‑140 (MQ=255)
            ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa                                                                                                                            >  2:394355/1‑140 (MQ=255)
            ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa                                                                                                                            >  2:589200/1‑140 (MQ=255)
            ttACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATa                                                                                                                            >  2:1158203/1‑140 (MQ=255)
              aCTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACg                                                                                                                          <  1:830902/140‑1 (MQ=255)
                         cttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTATACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttt                                                                                                               >  1:237646/1‑140 (MQ=255)
                         cttTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACttt                                                                                                               >  1:237644/1‑140 (MQ=255)
                              cTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTg                                                                                                          >  2:921241/1‑140 (MQ=255)
                                gAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGca                                                                                                        <  1:1013406/140‑1 (MQ=255)
                                 aaaaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCTGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGcac                                                                                                       >  1:1040599/1‑140 (MQ=255)
                                 aaaaTGCCACACATAAACAGAACTGTACCTCGTGTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGcac                                                                                                       >  2:200820/1‑140 (MQ=255)
                                   aaTGCCACACATAAACAGAACTGTACCTCGTTTAACGCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTTAACATGTCTCGTGCGGACGCACAAATCTCGATACGTCGTCGACTTTTTATGcagtg                                                                                                     >  2:895821/1‑137 (MQ=255)
                                   aaTGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCacac                                                                                                     >  1:495083/1‑140 (MQ=255)
                                         acacaTAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGc                                                                                               >  2:391209/1‑140 (MQ=255)
                                                aaCAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGaa                                                                                        <  2:346411/140‑1 (MQ=255)
                                                      aCTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATAcgcg                                                                                  <  1:664875/140‑1 (MQ=255)
                                                      aCTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATAcgcg                                                                                  <  1:601228/140‑1 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 >  2:931065/1‑140 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 <  1:1086527/140‑1 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 <  1:1086539/140‑1 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 <  1:590163/140‑1 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 <  1:1139790/140‑1 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 <  1:1163898/140‑1 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 <  1:895181/140‑1 (MQ=255)
                                                       cTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACATTTTTTGCACACTTATGCAACGGAATACgcgc                                                                                 <  1:895176/140‑1 (MQ=255)
                                                          tACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGa                                                                              >  1:336871/1‑140 (MQ=255)
                                                               cGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGACTACCACTATTTTTGCACACTTATGAAACGGAATACGCGCCGatatag                                                                         >  1:430490/1‑139 (MQ=255)
                                                                    aaCCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               <  1:499686/129‑1 (MQ=255)
                                                                    aaCCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               >  2:499686/1‑129 (MQ=255)
                                                                    aaCCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               >  2:499677/1‑129 (MQ=255)
                                                                    aaCCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               <  1:499677/129‑1 (MQ=255)
                                                                       ccGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGc                                                                 <  1:516271/140‑1 (MQ=255)
                                                                       ccGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGc                                                                 <  2:304913/140‑1 (MQ=255)
                                                                                 ttATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGc                                                       >  1:1059727/1‑140 (MQ=255)
                                                                                      acaTACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGtt                                                  <  2:561233/140‑1 (MQ=255)
                                                                                        aTACAACGATAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCa                                                <  2:1100957/140‑1 (MQ=255)
                                                                                         tACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAg                                               >  1:597248/1‑140 (MQ=255)
                                                                                           cAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGc                                             >  1:442036/1‑140 (MQ=255)
                                                                                            aaCGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                            <  1:887667/140‑1 (MQ=255)
                                                                                            aaCGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                            <  2:636656/140‑1 (MQ=255)
                                                                                            aaCGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                            <  1:940470/140‑1 (MQ=255)
                                                                                            aaCGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCa                                            <  1:887708/140‑1 (MQ=255)
                                                                                             aCGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAACTCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAg                                           <  1:553952/140‑1 (MQ=255)
                                                                                              cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACACACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               >  1:770338/1‑103 (MQ=255)
                                                                                              cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACACACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               >  1:770341/1‑103 (MQ=255)
                                                                                              cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACACACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               <  2:770338/103‑1 (MQ=255)
                                                                                              cGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACACACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCg                                                                               <  2:770341/103‑1 (MQ=255)
                                                                                                              ccATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTa                          <  1:1132811/140‑1 (MQ=255)
                                                                                                                        aCATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACa                <  1:356581/140‑1 (MQ=255)
                                                                                                                        aCATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACa                <  2:978398/140‑1 (MQ=255)
                                                                                                                         cATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCATATACTTAAACTGGCGTATTTTCTAACAt               <  1:328542/140‑1 (MQ=255)
                                                                                                                         cATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACAt               <  1:1183501/140‑1 (MQ=255)
                                                                                                                         cATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACAt               <  1:4744/140‑1 (MQ=255)
                                                                                                                         cATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACAt               <  1:388116/140‑1 (MQ=255)
                                                                                                                            gTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTAGTTTCTAACATAGt            >  1:1041108/1‑140 (MQ=255)
                                                                                                                                      gACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGtt  >  2:83181/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
ACGTCAACGTGTTTACTCAGGACTTCTTTACTGAAAATGCCACACATAAACAGAACTGTACCTCGTTTAACCCGAAATCTGTTATCACATACAACGTTAATTATGGCAAGCCATTTTTAAACATGTCTCGATCAGACGCACAAATCTCGATACGTCTTCCACTTTTTTTGCACACTTATGCAACGGAATACGCGCCGATATATCATTGTGCTTAACCTTGCCAGTTCAGGCAGATACTTAAACTGGCGTATTTTCTAACATAGTTCAATAAGTT  >  NC_000913/1907545‑1907818

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: